High-order accurate equations describing vibrations of thin bars

A method for deriving one-dimensional wave propagation equations in thin inhomogeneous anisotropic bars based on the mathematical homogenization theory for periodic media is used to obtain equations governing the longitudinal and transverse vibrations of a homogeneous circular bar. The equations are derived up to O(ε⁸) terms and take into account variable body forces and surface loads. Here, ε is the ratio of the bar’s typical thickness to the typical wavelength.     

薄膜外延生长的计算机模拟

以Cu膜为例，用Monte-Carlo算法模拟了薄膜生长的随机过程，并提出了更加完善的模型.在合理选择原子间相互作用计算方法的基础上，考虑了原子的吸附、在生长表面的迁移及迁移所引起的近邻原子连带效应、从生长表面的脱附等过程.模拟计算了薄膜的早期成核情况以及表面粗糙度和相对密度.结果表明，随着衬底温度的升高或入射率的降低，沉积在衬底上的原子逐步由离散型分布向聚集状态过渡形成一些岛核，并且逐步由二维岛核向三维岛核过渡.在一定的原子入射率下，存在三个优化温度，成核率最高时的最大成核温度Tn、薄膜的表面粗糙度最低 关键词： Monte-Carlo算法 计算机模拟 薄膜生长     

A relaxed model and its homogenization for nematic liquid crystals in composite materials

We analyse a model for equilibrium configurations of composite systems of nematic liquid crystal with polymer inclusions, in the presence of an external magnetic field. We assume that the system has a periodic structure, and consider the relaxed problem on the unit length constraint of the nematic director field. The relaxation of the Oseen–Frank energy functional is carried out by including bulk as well as surface energy penalty terms, rendering the problem fully non‐linear. We employ two‐scale convergence methods to obtain effective configurations of the system, as the size of the polymeric inclusions tends to zero. We discuss the minimizers of the effective energies for, both, the constrained as well as the unconstrained models. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling of an underground waste disposal site by upscaling

In this paper, we study the global behaviour of an underground waste disposal in order to have an accurate upscaled model suitable for the computations involved in safety assessment processes. We start from a detailed model describing the transport of pollutant leaking from a high number of units. Using the method of homogenization, going to the limit, we obtain first a macroscopic model where the sources are now appearing globally. Then we compute a first‐order matched asymptotic expansion and we give the error estimates for this approximation. Copyright © 2004 John Wiley & Sons, Ltd.     

ΓD‐convergence for a class of quasilinear elliptic equations in thin structures

We study the asymptotic behaviour of the solution of a stationary quasilinear elliptic problem posed in a domain Ω^(ε) of asymptotically degenerating measure, i.e. meas Ω^(ε) → 0 as ε → 0, where ε is the parameter that characterizes the scale of the microstructure. We obtain the convergence of the solution and the homogenized model of the problem is constructed using the notion of convergence in domains of degenerating measure. Proofs are given using the method of local characteristics of the medium Ω^(ε) associated with our problem in a variational form. Copyright © 2005 John Wiley & Sons, Ltd.     

Estimating tail probabilities of heavy tailed distributions with asymptotically zero relative error

Efficient estimation of tail probabilities involving heavy tailed random variables is amongst the most challenging problems in Monte-Carlo simulation. In the last few years, applied probabilists have achieved considerable success in developing efficient algorithms for some such simple but fundamental tail probabilities. Usually, unbiased importance sampling estimators of such tail probabilities are developed and it is proved that these estimators are asymptotically efficient or even possess the desirable bounded relative error property. In this paper, as an illustration, we consider a simple tail probability involving geometric sums of heavy tailed random variables. This is useful in estimating the probability of large delays in M/G/1 queues. In this setting we develop an unbiased estimator whose relative error decreases to zero asymptotically. The key idea is to decompose the probability of interest into a known dominant component and an unknown small component. Simulation then focuses on estimating the latter ‘residual’ probability. Here we show that the existing conditioning methods or importance sampling methods are not effective in estimating the residual probability while an appropriate combination of the two estimates it with bounded relative error. As a further illustration of the proposed ideas, we apply them to develop an estimator for the probability of large delays in stochastic activity networks that has an asymptotically zero relative error.      

Bounds for non‐periodic mixtures of infinitely many materials

We prove bounds on the homogenized coefficients for general non‐periodic mixtures of an arbitrary number of isotropic materials, in the heat conduction framework. The component materials and their proportions are given through the Young measure associated to the sequence of coefficient functions. Upper and lower bounds inequalities are deduced in terms of algebraic relations between this Young measure and the eigenvalues of the H‐limit matrix. The proofs employ arguments of compensated compactness and fine properties of Young measures. When restricted to the periodic case, we recover known bounds. Copyright © 2005 John Wiley & Sons, Ltd.     

Slowing down in the three-dimensional three-state Potts glass with nearest neighbor exchange ± J: A Monte Carlo study

,Static and dynamic properties of the Potts model on the simple cubic lattice with nearest neighbor -interaction are obtained from Monte Carlo simulations in a temperature range where full thermal equilibrium still can be achieved (). For a lattice size L = 16, in this range finite size effects are still negligible, but the data for the spin glass susceptibility agree with previous extrapolations based on finite size scaling of very small lattices. While the static properties are compatible with a zero temperature transition, they certainly do not prove it. Unlike the Ising spin glass, the decay of the time-dependent order parameter is compatible with a simple Kohlrausch function, , while a power law prefactor cannot be distinguished. The Kohlrausch exponent y ( T ) decreases from at [0pt] to at [0pt] however. The relaxation time is compatible with the exponential divergence postulated by McMillan for spin glasses at their lower critical dimension, but the exponent that can be extracted still differs significantly from the theoretical value, . Thus the present results support the conclusion that the Potts spin glass in d = 3 dimensions differs qualitatively from the Ising spin glass. Received: 8 October 1997 / Accepted: 27 November 1997     

Dynamic Scaling of Ramified Clusters Formed on Liquid Surfaces

A comprehensive simulation model -deposition,diffusion, rotation, reaction and aggregation model is presented to simulate the formation processes of ramified clusters on liquid surfaces, where clusters can diffuse and rotate easily. The mobility (including diffusion and rotation) of clusters is related to its mass, which is given by Dm = Dos-γD and θm =′θos-γθ, respectively. The influence of the reaction probability on the kinetics and structure formation is included in the simulation model. We concentrate on revealing dynamic scaling during ramified cluster formation. For this purpose, the time evolution of the cluster density and the weight-average cluster size as well as the cluster-size distribution scaling function at different time are determined for various conditions. The dependence of the cluster density on the deposition flux and time-dependence of fractal dimension are also investigated. The obtained results are helpful in understanding the formation of clusters or thin film growth on liquid surfaces.     

Scattering spectra and colour of disperse systems of low-absorptive particles

In this article, the results of calculations of scattering spectra of low-absorptive spherical particles are presented. It was obtained using different approximate and numerical methods. The comparative evaluations of application of single scattering approximation, the average field approximation and of the methods of computer modeling on the basis of Monte-Carlo method on various scattering multiplicity were performed. The comparison demonstrates the good agreement of calculated results with experimental scattering spectra of eye lens. By using the obtained spectral data the colour coordinates are calculated. The regularities of changes in colour characteristics in dependence of disperse system parameters are analyzed.