In the first part of this series of two papers, we proposed a theoretical formalism that enables one to model and categorize heterogeneous materials (media) via two-point correlation functions S2 and introduced an eff...
详细信息
In the first part of this series of two papers, we proposed a theoretical formalism that enables one to model and categorize heterogeneous materials (media) via two-point correlation functions S2 and introduced an efficient heterogeneous-medium (re)construction algorithm called the “lattice-point” algorithm. Here we discuss the algorithmic details of the lattice-point procedure and an algorithm modification using surface optimization to further speed up the (re)construction process. The importance of the error tolerance, which indicates to what accuracy the media are (re)constructed, is also emphasized and discussed. We apply the algorithm to generate three-dimensional digitized realizations of a Fontainebleau sandstone and a boron-carbide/aluminum composite from the two-dimensional tomographic images of their slices through the materials. To ascertain whether the information contained in S2 is sufficient to capture the salient structural features, we compute the two-point cluster functions of the media, which are superior signatures of the microstructure because they incorporate topological connectedness information. We also study the reconstruction of a binary laser-speckle pattern in two dimensions, in which the algorithm fails to reproduce the pattern accurately. We conclude that in general reconstructions using S2 only work well for heterogeneous materials with single-scale structures. However, two-point information via S2 is not sufficient to accurately model multiscale random media. Moreover, we construct realizations of hypothetical materials with desired structural characteristics obtained by manipulating their two-point correlation functions.
Kinetic Monte Carlo(KMC)is a stochastic model used to simulate crystal ***,most KMC models rely on a pre-defined lattice that neglects dislocations,lattice mismatch and strain *** this paper,we investigate the use of ...
详细信息
Kinetic Monte Carlo(KMC)is a stochastic model used to simulate crystal ***,most KMC models rely on a pre-defined lattice that neglects dislocations,lattice mismatch and strain *** this paper,we investigate the use of a 3D off-lattice KMC *** test this method by investigating impurity diffusion in a strained FCC *** faster than a molecular dynamics simulation,the most general implementation of off-lattice KMC is much slower than a lattice-based *** improved procedure is achieved for weakly strained systems by precomputing approximate saddle point locations based on unstrained lattice *** this way,one gives up some of the flexibility of the general method to restore some of the computational speed of lattice-based *** addition to providing an alternative approach to nano-materials simulation,this type of simulation will be useful for testing and calibrating methods that seek to parameterize the variation in the transition rates for lattice-based KMC using continuum modeling.
On March 11, 1944, the famous Eremitani Church in Padua (Italy) was destroyed in an Allied bombing along with the inestimable frescoes by Andrea Mantegna et al. contained in the Ovetari Chapel. In the last 60 years, s...
详细信息
ISBN:
(纸本)9783540728221
On March 11, 1944, the famous Eremitani Church in Padua (Italy) was destroyed in an Allied bombing along with the inestimable frescoes by Andrea Mantegna et al. contained in the Ovetari Chapel. In the last 60 years, several attempts have been made to restore the fresco fragments by traditional methods, but without much success. We have developed an efficient pattern recognition algorithm to map the original position and orientation of the fragments, based on comparisons with an old gray level image of the fresco prior to the damage. This innovative technique allowed for the partial reconstruction of the frescoes. Unfortunately, the surface covered by the fragments is only 77 m 2, while the original area was of several hundreds. This means that we can reconstruct only a fraction (less than 8%) of this inestimable artwork. In particular the original color of the blanks is not known. This begs the question of whether it is possible to estimate mathematically the original colors of the frescoes by making use of the potential information given by the available fragments and the gray level of the pictures taken before the damage. Can one estimate how faithful such restoration is? In this paper we retrace the development of the recovery of the frescoes as an inspiring and challenging real-life problem for the development of new mathematical methods. We introduce two models for the recovery of vector valued functions from incomplete data, with applications to the fresco recolorization problem. The models are based on the minimization of a functional which is formed by the discrepancy with respect to the data and additional regularization constraints. The latter refer to joint sparsity measures with respect to frame expansions for the first functional and functional total variation for the second. We establish the relations between these two models. As a byproduct we develop the basis of a theory of fidelity in color recovery, which is a crucial issue in art restoration and
A photonic quasicrystal consists of two or more dielectric materials arranged in a quasiperiodic pattern with noncrystallographic symmetry that has a photonic band gap. We use a novel method to find the pattern with t...
详细信息
A photonic quasicrystal consists of two or more dielectric materials arranged in a quasiperiodic pattern with noncrystallographic symmetry that has a photonic band gap. We use a novel method to find the pattern with the widest TM-polarized gap for two-component materials. Patterns are obtained by computing a finite sum of density waves, assigning regions where the sum exceeds a threshold to a material with one dielectric constant, ϵ1, and all other regions to another, ϵ0. Compared to optimized crystals, optimized quasicrystals have larger gaps at low constrasts ϵ1/ϵ0 and have gaps that are much more isotropic for all contrasts. For high contrasts, optimized hexagonal crystals have the largest gaps.
We describe a computational framework linking uncertainty quantification (UQ) methods for continuum problems depending on random parameters with equation-free (EF) methods for performing continuum deterministic numeri...
详细信息
We consider the construction of point processes from tilings, with equal-volume tiles, of d-dimensional Euclidean space Rd. We show that one can generate, with simple algorithms ascribing one or more points to each ti...
详细信息
We consider the construction of point processes from tilings, with equal-volume tiles, of d-dimensional Euclidean space Rd. We show that one can generate, with simple algorithms ascribing one or more points to each tile, point processes which are “superhomogeneous” (or “hyperuniform”)—i.e., for which the structure factor S(k) vanishes when the wave vector k tends to zero. The exponent γ characterizing the leading small-k behavior, S(k→0)∝kγ, depends in a simple manner on the nature of the correlation properties of the specific tiling and on the conservation of the mass moments of the tiles. Assigning one point to the center of mass of each tile gives the exponent γ=4 for any tiling in which the shapes and orientations of the tiles are short-range correlated. Smaller exponents in the range 4−d<γ<4 (and thus always superhomogeneous for d≤4) may be obtained in the case that the latter quantities have long-range correlations. Assigning more than one point to each tile in an appropriate way, we show that one can obtain arbitrarily higher exponents in both cases. We illustrate our results with explicit constructions using known deterministic tilings, as well as some simple stochastic tilings for which we can calculate S(k) exactly. Our results provide an explicit analytical construction of point processes with γ>4. Applications to condensed matter physics, and also to cosmology, are briefly discussed.
This paper gives a systematic introduction to HMM,the heterogeneous multiscale methods,including the fundamental design principles behind the HMM philosophy and the main obstacles that have to be overcome when using H...
详细信息
This paper gives a systematic introduction to HMM,the heterogeneous multiscale methods,including the fundamental design principles behind the HMM philosophy and the main obstacles that have to be overcome when using HMM for a particular *** is illustrated by examples from several application areas,including complex fluids,micro-fluidics,solids,interface problems,stochastic problems,and statistically self-similar *** is given to the technical tools,such as the various constrained molecular dynamics,that have been developed,in order to apply HMM to these *** of mathematical results on the error analysis of HMM are *** review ends with a discussion on some of the problems that have to be solved in order to make HMM a more powerful tool.
This paper presents a systematic approach for finding efficient boundary conditions for molecular dynamics simulations of crystalline solids. These boundary conditions effectively eliminate phonon reflection at the bo...
This paper presents a systematic approach for finding efficient boundary conditions for molecular dynamics simulations of crystalline solids. These boundary conditions effectively eliminate phonon reflection at the boundary and at the same time allow the thermal energy from the bath to be introduced to the system. Our starting point is the Mori-Zwanzig formalism [R. Zwanzig, J. Chem. Phys. 32, 1173 (1960); in Systems Far from Equilibrium, edited by L. Garrido (Interscience, New York, 1980); H. Mori, Prog. Theor. Phys. 33, 423 (1965)] for eliminating the thermal bath, but we take the crucial next step that goes beyond this formalism in order to obtain memory kernels that decay faster. An equivalent variational formulation allows us to find the optimal approximate boundary conditions, after specifying the spatial-temporal domain of dependence for the positions of the boundary atoms. Application to a one-dimensional chain, a two-dimensional Lennard-Jones system, and a three-dimensional model of α-iron with embedded atom potential is presented to demonstrate the effectiveness of this approach.
We use an “equation-free,” coarse-grained computational approach to accelerate molecular dynamics-based computations of demixing (segregation) of dissimilar particles subject to an upward gas flow (gas-fluidized bed...
详细信息
We use an “equation-free,” coarse-grained computational approach to accelerate molecular dynamics-based computations of demixing (segregation) of dissimilar particles subject to an upward gas flow (gas-fluidized beds). We explore the coarse-grained dynamics of these phenomena in gently fluidized beds of solid mixtures of different densities, typically a slow process for which reasonable continuum models are currently unavailable.
暂无评论