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.

Implication of scaling hierarchy associated with nonequilibrium: field and particulate

The advent of nanotechnology has necessitated a better understanding of how material microstructure changes at the atomic level would affect the macroscopic properties that control the performance. Such a challenge has uncovered many phenomena that were not previously understood and taken for granted. Among them are the basic foundation of dislocation theories which are now known to be inadequate. Simplifying assumptions invoked at the macroscale may not be applicable at the micro- and/or nanoscale. There are implications of scaling hierrachy associated with in-homegeneity and nonequilibrium. of physical systems. What is taken to be homogeneous and equilibrium at the macroscale may not be so when the physical size of the material is reduced to microns. These fundamental issues cannot be dispensed at will for the sake of convenience because they could alter the outcome of predictions. Even more unsatisfying is the lack of consistency in modeling physical systems. This could translate to the inability for identifying the relevant manufacturing parameters and rendering the end product unpractical because of high cost. Advanced composite and ceramic materials are cases in point. Discussed are potential pitfalls for applying models at both the atomic and continuum levels. No encouragement is made to unravel the truth of nature. Let it be partiuclates, a smooth continuum or a combination of both. The present trend of development in scaling tends to seek for different characteristic lengths of material microstructures with or without the influence of time effects. Much will be learned from atomistic simulation models to show how results could differ as boundary conditions and scales are changed. Quantum mechanics, continuum and cosmological models provide evidence that no general approach is in sight. Of immediate interest is perhaps the establishment of greater precision in terminology so as to better communicate results involving multiscale physical events.

Mesofracture mechanics: a necessary link

Classical fracture mechanics is based on the premise that small scale features could be averaged to give a larger scale property such that the assumption of material homogeneity would hold. Involvement of the material microstructure, however, necessitates different characteristic lengths for describing different geometric features. Macroscopic parameters could not be freely exchanged with those at the microscopic scale level. Such a practice could cause misinterpretation of test data. Ambiguities arising from the lack of a more precise range of limitations for the definitions of physical parameters are discussed in connection with material length scales. Physical events overlooked between the macroscopic and microscopic scale could be the link that is needed to bridge the gap. The classical models for the creation of free surface for a liquid and solid are oversimplified. They consider only the translational motion of individual atoms. Movements of groups or clusters of molecules deserve attention. Multiscale cracking behavior also requires the distinction of material damage involving at least two different scales in a single simulation. In this connection, special attention should be given to the use of asymptotic solution in contrast to the full field solution when applying fracture criteria. The former may leave out detail features that would have otherwise been included by the latter. Illustrations are provided for predicting the crack initiation sites of piezoceramics. No definite conclusions can be drawn from the atomistic simulation models such as those used in molecular dynamics until the non-equilibrium boundary conditions can be better understood. The specification of strain rates and temperatures should be synchronized as the specimen size is reduced to microns. Many of the results obtained at the atomic scale should be first identified with those at the mesoscale before they are assumed to be connected with macroscopic observations. Hopefully, "mesofracture mechanics" could serve as the link to bring macrofracture mechanics closer to microfracture mechanics.

智能交通仿真与车辆Agent决策策略的研究

文章介绍了Agent技术及特点,根据Agent的原理及其仿真的优势,着重叙述了Agent技术在交通仿真中的应用,详细探讨了基于Agent的智能仿真系统中车辆Agent的决策特点及模糊决策方法,同时,还分析了驾驶员因素对驾驶行为的影响。

基于时间Petri网的制造过程模型研究

分析了时间Petri网的激发规则、托肯可用时间和抑制弧等特性,以及制造过程中随机故障的特征。提出不同的时间关联方式对应的多种建模方法,考虑不同的故障发现模式、不同的作业处理策略,建立相应的单机制造过程模型。在此基础上采用模块化和层次化方法可以构建复杂制造过程的时间着色Petri网模型,并可以转换成仿真模型,进一步分析随机机器故障对制造过程性能的影响。

Effects of Subgrid-Scale Modeling on Time Correlations in Large Eddy Simulation

The effects of the unresolved subgrid-scale (SGS) motions on the energy balance of the resolved scales in large eddy simulation (LES) have been investigated actively because modeling the energy transfer between the resolved and unresolved scales is crucial to constructing accurate SGS models. But the subgrid scales not only modify the energy balance, they also contribute to temporal decorrelation of the resolved scales. The importance of this effect in applications including the predictability problem and the evaluation of sound radiation by turbulent flows motivates the present study of the effect of SGS modeling on turbulent time correlations. This paper compares the two-point, two-time Eulerian velocity correlation in isotropic homogeneous turbulence evaluated by direct numerical simulation (DNS) with the correlations evaluated by LES using a standard spectral eddy viscosity. It proves convenient to express the two-point correlations in terms of spatial Fourier decomposition of the velocity field. The LES fields are more coherent than the DNS fields: their time correlations decay more slowly at all resolved scales of motion and both their integral scales and microscales are larger than those of the DNS field. Filtering alone is not responsible for this effect: in the Fourier representation, the time correlations of the filtered DNS field are identical to those of the DNS field itself. The possibility of modeling the decorrelating effects of the unresolved scales of motion by including a random force in the model is briefly discussed. The results could have applications to the problem of computing sound sources in isotropic homogeneous turbulence by LES

Large Eddy Simulation of Complex Turbulent Flows: Physical Aspects and Research Trends

In the current paper, we have primarily addressed one powerful simulation tool developed during the last decades-Large Eddy Simulation (LES), which is most suitable for unsteady three-dimensional complex turbulent flows in industry and natural environment. The main point in LES is that the large-scale motion is resolved while the small-scale motion is modeled or, in geophysical terminology, parameterized. With a view to devising a subgrid-scale(SGS) model of high quality, we have highlighted analyzing physical aspects in scale interaction and-energy transfer such as dissipation, backscatter, local and non-local interaction, anisotropy and resolution requirement. They are the factors responsible for where the advantages and disadvantages in existing SGS models come from. A case study on LES of turbulence in vegetative canopy is presented to illustrate that LES model is more based on physical arguments. Then, varieties of challenging complex turbulent flows in both industry and geophysical fields in the near future-are presented. In conclusion; we may say with confidence that new century shall see the flourish in the research of turbulence with the aid of LES combined with other approaches.

Modeling and Simulation of Liquid Pulsed Particulate Fluidized Beds

Pulsed fluidization is of considerable interest in process engineering for improving fluidization quality. Quantitative understanding of the pulsed two-phase flow behaviors is very important for proper design and optimum operation of such contactors. The

Simulation of Coal Combustion by AUSM Turbulence-Chemistry Char Combustion Model and a Full Two-Fluid Model

An algebraic unified second-order moment (AUSM) turbulence-chemistry model of char combustion is introduced in this paper, to calculate the effect of particle temperature fluctuation on char combustion. The AUSM model is used to simulate gas-particle flows, in coal combustion in a pulverized coal combustor, together with a full two-fluid model for reacting gas-particle flows and coal combustion, including the sub-models as the k-epsilon-k(p) two-phase turbulence niodel, the EBU-Arrhenius volatile and CO combustion model, and the six-flux radiation model. A new method for calculating particle mass flow rate is also used in this model to correct particle outflow rate and mass flow rate for inside sections, which can obey the principle of mass conservation for the particle phase and can also speed up the iterating convergence of the computation procedure effectively. The simulation results indicate that, the AUSM char combustion model is more preferable to the old char combustion model, since the later totally eliminate the influence of particle temperature fluctuation on char combustion rate.

Effects Of Subgrid-Scale Modeling On Lagrangian Statistics In Large-Eddy Simulation

The application of large-eddy simulation (LES) to turbulent transport processes requires accurate prediction of the Lagrangian statistics of flow fields. However, in most existing SGS models, no explicit consideration is given to Lagrangian statistics. In this paper, we focus on the effects of SGS modeling on Lagrangian statistics in LES ranging from statistics determining single-particle dispersion to those of pair dispersion and multiparticle dispersion. Lagrangian statistics in homogeneous isotropic turbulence are extracted from direct numerical simulation (DNS) and the LES with a spectral eddy-viscosity model. For the case of longtime single-particle dispersion, it is shown that, compared to DNS, LES overpredicts the time scale of the Lagrangian velocity correlation but underpredicts the Lagrangian velocity fluctuation. These two effects tend to cancel one another leading to an accurate prediction of the longtime turbulent dispersion coefficient. Unlike the single-particle dispersion, LES tends to underestimate significantly the rate of relative dispersion of particle pairs and multiple-particles, when initial separation distances are less than the minimum resolved scale due to the lack of subgrid fluctuations. The overprediction of LES on the time scale of the Lagrangian velocity correlation is further confirmed by a theoretical analysis using a turbulence closure theory.

On the computation of space-time correlations by large-eddy simulation

The effect of subgrid-scale (SGS) modeling on velocity (space-) time correlations is investigated in decaying isotropic turbulence. The performance of several SGS models is evaluated, which shows superiority of the dynamic Smagorinsky model used in conjunction with the multiscale large-eddy simulation (LES) procedure. Compared to the results of direct numerical simulation, LES is shown to underpredict the (un-normalized) correlation magnitude and slightly overpredict the decorrelation time scales. This can lead to inaccurate solutions in applications such as aeroacoustics. The underprediction of correlation functions is particularly severe for higher wavenumber modes which are swept by the most energetic modes. The classic sweeping hypothesis for stationary turbulence is generalized for decaying turbulence and used to analyze the observed discrepancies. Based on this analysis, the time correlations are determined by the wavenumber energy spectra and the sweeping velocity, which is the square root of the total energy. Hence, an accurate prediction of the instantaneous energy spectra is most critical to the accurate computation of time correlations. (C) 2004 American Institute of Physics.

Modeling, numerical methods, and simulation for particle-fluid two-phase flow problems

In this paper, we study the issues of modeling, numerical methods, and simulation with comparison to experimental data for the particle-fluid two-phase flow problem involving a solid-liquid mixed medium. The physical situation being considered is a pulsed liquid fluidized bed. The mathematical model is based on the assumption of one-dimensional flows, incompressible in both particle and fluid phases, equal particle diameters, and the wall friction force on both phases being ignored. The model consists of a set of coupled differential equations describing the conservation of mass and momentum in both phases with coupling and interaction between the two phases. We demonstrate conditions under which the system is either mathematically well posed or ill posed. We consider the general model with additional physical viscosities and/or additional virtual mass forces, both of which stabilize the system. Two numerical methods, one of them is first-order accurate and the other fifth-order accurate, are used to solve the models. A change of variable technique effectively handles the changing domain and boundary conditions. The numerical methods are demonstrated to be stable and convergent through careful numerical experiments. Simulation results for realistic pulsed liquid fluidized bed are provided and compared with experimental data. (C) 2004 Elsevier Ltd. All rights reserved.

Chemical equilibrium analysis and CFD simulation of coal combustion and NOx formation

A full two-fluid model of reacting gas-particle flows with an algebraic unified second-order moment (AUSM) turbulence-chemistry model is used to simulate Beijing coal combustion and NOx formation. The sub-models are the k-epsilon-kp two-phase turbulence model, the EBU-Arrhenius volatile and CO combustion model, the six-flux radiation model, coal devolatilization model and char combustion model. The blocking effect on NOx formation is discussed. In addition, the chemical equilibrium analysis is used to predict NOx concentration at different temperature. Results of CID simulation and chemical equilibrium analysis show that, optimizing air dynamic parameters can delay the NOx formation and decrease NOx emission, but it is effective only in a restricted range. In order to decrease NOx emission near to zero, the re-burning or other chemical methods must be used.

Processing Simulation of Virtual Manufacture in Laser Surface Treatmenthang

Processing simulation is at the bottom of the coral technology of VM and is also difficult due to the complexity of mechanism and diversity of parameters. Previously much research has been mainly carried out on the geometrical simulation or physical simulation respectively. The aim of this paper is to study the processing simulation in laser surface treatment based on the mechanism, put forward the architecture of the whole processing simulation and give the models of the processing. As a result the data structure layers in the whole simulation is presented.

Comparisons of static, quasi-static and dynamic 3D porous media scale network models for two-phase immiscible flow in porous media

A dynamic 3D pore-scale network model is formulated for investigating the effect of interfacial tension and oil-water viscosity during chemical flooding. The model takes into account both viscous and capillary forces in analyzing the impact of chemical properties on flow behavior or displacement configuration, while the static model with conventional invasion percolation algorithm incorporates the capillary pressure only. From comparisons of simulation results from these models. it indicates that the static pore scale network model can be used successfully when the capillary number is low. With the capillary increases due to the enhancement of water viscosity or decrease of interfacial tension, only the quasi-static and dynamic model can give insight into the displacement mechanisms.