Boundary continuation methods

Pseudomorphism received methodical study from about 1905. A micro-section taken across the interface between a substrate and an electrodeposit shows the grain boundaries of the former continue across the interface into the deposit (Fig. 12.5). As grain boundaries are internal faces of metal crystals, when they continue into the deposit the latter is displaying the form of the substrate. Hothersall s 1935 paper contains numerous excellent illustrations with substrates and deposits chosen from six different metals, crystallising in different lattice systems and with different equilibrium spacing. Grain boundary continuation and hence pseudomorphism is evident despite the differences. 

In this section we apply the adaptive boundary value solution procedure and the pseudo-arclength continuation method to a set of strained premixed hydrogen-air flames. Our goal is to predict accurately and efficiently the extinction behavior of these flames as a function of the strain rate and the equivalence ratio. Detailed transport and complex chemical kinetics are included in all of the calculations. The reaction mechanism for the hydrogen-air system is listed in Table... 

We have shown that multiple travelling front waves can occur in a reaction-diffusion-convection system. These waves can be studied in an unbounded system by using a wave transformation and solving a special boundary value problem with the use of continuation methods. These results provide various parameter dependences of the velocity of the wave. Moreover, in a bounded system the waves move back and forth through the system and form remarkable zig-zag patterns. 

By Fourier analysis one can easily show that for b=0 the linear mapping J is nonsingular and the implicit function theorem can be applied. As long as J is nonsingular we can apply a continuation method for b>0 and use the implicit function theorem to calculate the free boundary. 

Several procedures for workspace analysis of manipulators have been proposed iterative determination of reachable points by means of matrix formulation of direct problem, [14-16], or through probabilistic techniques, [17, 18], or continuation methods, [19] dynamical evaluation of extreme configuration, [20, 21] determination of boundary surfaces for Jacobian domain, [22-24] algebraic formulation for specific manipulators, [25, 26]. However, in order to facilitate the numerical solution for the design problem of Eqs. (7), (8) and (9) it can be useful to express the involved workspace characteristics by means of a suitable analytical formulation. 

A discontinuous element plays an important role in the boundary element method because the problem variables are not forced to be continuous across the elements. The major benefit of discontinuous elements is their ability to model discontinuous stress results. They are therefore very usefiil in fracture mechanics analysis for modelling the stress behaviour at a crack front. 

Methods for numerical analyses such as tlris can be obtained from commercial software, and the advent of the computer has considerably eased the work required to obtain numerical values for heat distribution and profiles in a short time, or even continuously if a monitor supplies the boundary values of heat content or temperature during an operation. 

These apparent restrictions in size and length of simulation time of the fully quantum-mechanical methods or molecular-dynamics methods with continuous degrees of freedom in real space are the basic reason why the direct simulation of lattice models of the Ising type or of solid-on-solid type is still the most popular technique to simulate crystal growth processes. Consequently, a substantial part of this article will deal with scientific problems on those time and length scales which are simultaneously accessible by the experimental STM methods on one hand and by Monte Carlo lattice simulations on the other hand. Even these methods, however, are too microscopic to incorporate the boundary conditions from the laboratory set-up into the models in a reahstic way. Therefore one uses phenomenological models of the phase-field or sharp-interface type, and finally even finite-element methods, to treat the diffusion transport and hydrodynamic convections which control a reahstic crystal growth process from the melt on an industrial scale. 

In Planck s investigation of equilibrium in dilute solutions, the law of Henry follows as a deduction, whereas in van t Hoff s theory, based on the laws of osmotic pressure ( 128), it must be introduced as a law of experience. The difference lies in the fact that in Planck s method the solution is converted continuously into a gas mixture of known potential, whilst in van t Hoff s method it stands in equilibrium with a gas of known potential, and the boundary condition (Henry s law) must be known as well. Planck (Thermodynamik, loc. cit.) also deduces the laws of osmotic pressure from the theory. 

The characteristic feature of solid—solid reactions which controls, to some extent, the methods which can be applied to the investigation of their kinetics, is that the continuation of product formation requires the transportation of one or both reactants to a zone of interaction, perhaps through a coherent barrier layer of the product phase or as a monomolec-ular layer across surfaces. Since diffusion at phase boundaries may occur at temperatures appreciably below those required for bulk diffusion, the initial step in product formation may be rapidly completed on the attainment of reaction temperature. In such systems, there is no initial delay during nucleation and the initial processes, perhaps involving monomolec-ular films, are not readily identified. The subsequent growth of the product phase, the main reaction, is thereafter controlled by the diffusion of one or more species through the barrier layer. Microscopic observation is of little value where the phases present cannot be unambiguously identified and X-ray diffraction techniques are more fruitful. More recently, the considerable potential of electron microprobe analyses has been developed and exploited. 

The boundary conditions are unchanged. The method of lines solution continues to use a second-order approximation for dajdr and merely adds a Vr term to the coefficients for the points at r Ar. 

The method of test functions is quite applicable in verifying convergence and determining the order of accuracy and is stipulated by a proper choice of the function I7(x). Such a function is free to be chosen in any convenient way so as to provide the validity of the continuity conditions at every discontinuity point of coefficients. By inserting it in equation (1) of Section 1 we are led to the right-hand side / = kU ) — qU and the boundary values jj, — U(0) and = U 1). The solution of such a problem relies on scheme (4) of Section 1 and then the difference solution will be compared with a known function U x) on various grids. 

---