Local displacement, time dependence

Besides the remarkable directionality of the motion, the images also demonstrate a periodic variation of the cluster from an elongated to a circular shape (Fig. 39). The diagrams in Fig. 39 depict the time dependence of the displacement and the cluster size. Until the cluster was finally trapped, the speed remained fairly constant as can be seen from the constant slope in Fig. 39 a. The oscillatory variation of the cluster shape is shown in Fig. 39b. Although a coarse model for the motion has been presented in Fig. 39, the actual cause of the motion remains unknown. The ratchet model proposed by J. Frost requires a non-equiUb-rium variation in the energetic potential to bias the Brownian motion of a molecule or particle under anisotropic boundary conditions [177]. Such local perturbations of the molecular structure are believed to be caused by the mechanical contact with the scaiming tip. A detailed and systematic study of this question is still in progress. 

A non-perturbative theory of the multiphonon relaxation of a localized vibrational mode, caused by a high-order anharmonic interaction with the nearest atoms of the crystal lattice, is proposed. It relates the rate of the process to the time-dependent non-stationary displacement correlation function of atoms. A non-linear integral equation for this function is derived and solved numerically for 3- and 4-phonon processes. We have found that the rate exhibits a critical behavior it sharply increases near a specific (critical) value(s) of the interaction. 

The relaxation of a local mode is characterized by the time-dependent anomalous correlations the rate of the relaxation is expressed through the non-stationary displacement correlation function. The non-linear integral equations for this function has been derived and solved numerically. In the physical meaning, the equation is the self-consistency condition of the time-dependent phonon subsystem. We found that the relaxation rate exhibits a critical behavior it is sharply increased near a specific (critical) value(s) of the interaction the corresponding dependence is characterized by the critical index k — 1, where k is the number of the created phonons. In the close vicinity of the critical point(s) the rate attains a very high value comparable to the frequency of phonons. 

Therefore, before describing the modification of the equilibrium FDT, we need to study in details the behavior of D(t). Note, however, that the integrated velocity correlation function [, Cvv(/) df takes on the meaning of a time-dependent diffusion coefficient only when the mean-square displacement increases without bounds (when the particle is localized, this quantity characterizes the relaxation of the mean square displacement Ax2 t) toward its finite limit Ax2(oo)). 

For numerical investigations of stress localizations in laminates, the discretizational effort can be reduced significantly if only the boundary needs to be discretized, as it is for e -ample the case in the classical boundary element method (BEM). But in this method a fundamental solution is needed which is in many cases difficult to achieve or even unknown. The Boundary Finite Element Method (BFEM) to be presented here does not require such a fundamental solution, because the element formulation is based on the finite element method (FEM), Thus the BFEM can be characterized to be a finite element based boundary discretization method. This method was originally developed from Wolf and Song [10] under the name Consistent Finite Element Cell Method for time-dependent problems in soil-mechanics. The basic assumption of this method is that a stiffness matrix describing the force-displacement relation at discrete degrees of freedom at the boundary of the continuum is scalable with respect to one point in three-dimensional space, the so-called similarity center, if similar contours within the continuum are considered. In contrast to this, the current work deals with the case of equivalent cross-sectional properties, i.e., that cross-sections parallel to the boundary can be described by the same stiffness matrix, which is the appropriate formulation for the case of the free-edge effect and the matrix crack problem. The boundary stiffness matrix results from a Matrix-Riccati equation. The field quantities inside of the continuum can be calculated from an ordinary differential equation. 

One by one, each of the N atoms is displaced randomly, and the closest local minima is determined. If the new structure has a lower total energy than the original one, this new one is kept, and the old one discarded. This is repeated approximately 500 1000 times (depending on cluster size). 

Lattice modes, or external modes, as well as internal modes produce harmonic displacements of the atoms. The internal modes result in local molecular deformation, whilst external modes are supposed to entrain the atoms when a molecule is rigidly displaced from its equilibrium position in the lattice. The time-dependent atomic position vector r can be expressed in terms of the internal displacement vector u nx, taken with respect to the molecular centre of mass, and the external displacement vector, Hext the displacement vectors have units of length. A, Eq. (A2.52). 

Note that similar methods can be applied also to the solution of systems arising from the time-dependent heat conduction - convection problems. There are also other space decomposition methods. Let us mention the displacement decomposition technique for solving the elasticity problems and composite grid technique, for solving problems, which need a local resolution. More details can be found in Blaheta (2002) and Blaheta et al. (2002b). 

The time dependence of local displacement c and displacement velocity dc/dt at the craze or crack initiation is shown schematically in Figure 1. The velocity drops to one half of its maximum value at ti/2 If one puts t = 0 at the maximum and assumes a symmetric time dependence of displacement velocity, the half width of velocity distribution at such an elementary act is 2ti/2. The Fourier transform (FT) of velocity yields the spectral distribution of the emitted acoustic burst. The intensity is a maximum at zero frequency and drops to half this value at (01/2 = 2mv /2> In the first approximation the product of the half width of velocity and frequency distribution of acoustic emission, 2ti/2 X 2vi/2, equals 0.8825. To have a substantial amount of energy available in the frequency range of 1 MHz ( —1 1/2) the displacement velocity curve vs. time must have a half width of 0.5 /xsec, i.e., the major part of the local displacement must occur within 0.5 /i.sec. The square of the maximum value of velocity times... 

The province of conventional dielectric measurements is here taken to be the determination of the relations of the polarization E and current density J. to the electric field in the macroscopic Maxwell equations. Proper theory should account for these relations in condensed phases as a function of state variables time dependence of applied fields and molecular parameters by appropriate statistical averaging over molecular displacements determined by the equations of motion in terms of molecular forces and fields. Simplifying assumptions and approximations are of course necessary. One kind often made and debated is use of an effective or mean local field at a molecule rather than the sum of microscopic... 

In addition to internal charge transport in transient cases the external current can also be compensated by dielectric effects (displacement current by local polarization). This transient charging current is characterized by electrical capacitances. They are considered in more detail in Section 7.3.3. Capacitive effects (charge storage) are generally responsible for time dependences. Apart from dielectric effects, storage phenomena can also occur if the stoichiometry changes by virtue of the current flow. Such chemical capacitances (see Section 6.7.4) will be treated in more detail in Section 7.3.4 (cf. also Section 7.4). 

For this algorithm, one can prove that detailed balance is guaranteed and the exact average of any configuration-dependent property over the accessible space is obtained. Two key issues determine the detailed balance. The first is the fact that the trial probability to pick the displacement vector Dfc to go from the fcth to the Zth e-sphere equals the trial probability to pick the displacement vector D fc for the reverse step. The second issue is that the trial probability for a local MC step that moves the walker from a point inside an e-sphere to a point outside that sphere is the same as for the reverse move i.e., (1 - / ) times what it would be in a walk restricted to local moves. 

Relaxation dispersion data for water on Cab-O-Sil, which is a monodis-perse silica fine particulate, are shown in Fig. 2 (45). The data are analyzed in terms of the model summarized schematically in Fig. 3. The y process characterizes the high frequency local motions of the liquid in the surface phase and defines the high field relaxation dispersion. There is little field dependence because the local motions are rapid. The p process defines the power-law region of the relaxation dispersion in this model and characterizes the molecular reorientations mediated by translational displacements on the length scale of the order of the monomer size, or the particle size. The a process represents averaging of molecular orientations by translational displacements on the order of the particle cluster size, which is limited to the long time or low frequency end by exchange with bulk or free water. This model has been discussed in a number of contexts and extended studies have been conducted (34,41,43). 

---