Rotating-frame nuclear magnetic relaxation of spins diffusing on a disordered lattice: a Monte Carlo model

The rotating-frame nuclear magnetic relaxation rate of spins diffusing on a disordered lattice has been calculated by Monte Carlo methods. The disorder includes not only variation in the distances between neighbouring spin sites but also variation in the hopping rate associated with each site. The presence of the disorder, particularly the hopping rate disorder, causes changes in the time-dependent spin correlation functions which translate into asymmetry in the characteristic peak in the temperature dependence of the dipolar relaxation rate. The results may be used to deduce the average hopping rate from the relaxation but the effect is not sufficiently marked to enable the distribution of the hopping rates to be evaluated. The distribution, which is a measure of the degree of disorder, is the more interesting feature and it has been possible to show from the calculation that measurements of the relaxation rate as a function of the strength of the radiofrequency spin-locking magnetic field can lead to an evaluation of its width. Some experimental data on an amorphous metal - hydrogen alloy are reported which demonstrate the feasibility of this novel approach to rotating-frame relaxation in disordered materials.

Simulation of transverse and longitudinal magnetic ripple structures induced by surface anisotropy

Micromagnetic ripple structures on the surfaces of thick specimens of ultra-soft magnetic material having strong surface anisotropy Ks favouring out-of-surface magnetization have been calculated. These ripples have wavelengths of the order of 0.1 μm and extend to a depth ∼ √A/Ms, where A is the exchange constant and Ms is the saturation magnetization. The wave-vectors of the ripple structures are either transverse or parallel to the bulk magnetization. Both structures have lower energy than the one-dimensional structure discussed by O'Handley and Woods, and they exhibit stronger normal magnetization. The transverse structure requires a surface anisotropy Ks ≥ 0.80K0, where is that required for the one-dimensional structure. The threshold for longitudinal ripples is 0.84K0. It is suggested that the transverse structure probably constitutes the ground state. The magnitudes of Ks and A should be obtainable from measurements of the ripple wavelength and amplitude, and Ms.

The implementation of an aeronautical CFD flow code onto distributed memory parallel systems

The parallelization of an industrially important in-house computational fluid dynamics (CFD) code for calculating the airflow over complex aircraft configurations using the Euler or Navier–Stokes equations is presented. The code discussed is the flow solver module of the SAUNA CFD suite. This suite uses a novel grid system that may include block-structured hexahedral or pyramidal grids, unstructured tetrahedral grids or a hybrid combination of both. To assist in the rapid convergence to a solution, a number of convergence acceleration techniques are employed including implicit residual smoothing and a multigrid full approximation storage scheme (FAS). Key features of the parallelization approach are the use of domain decomposition and encapsulated message passing to enable the execution in parallel using a single programme multiple data (SPMD) paradigm. In the case where a hybrid grid is used, a unified grid partitioning scheme is employed to define the decomposition of the mesh. The parallel code has been tested using both structured and hybrid grids on a number of different distributed memory parallel systems and is now routinely used to perform industrial scale aeronautical simulations. Copyright © 2000 John Wiley & Sons, Ltd.

Parallelisation of electromagnetic simulation codes

In this paper results obtained from the parallelisation of existing 3D electromagnetic Finite Element codes within the ESPRIT HPCN project PARTEL are presented. The parallelisation procedure, based on the Bulk Synchronous Parallel approach, is outlined and the encouraging results obtained in terms of speed-up on some industrially significant test cases are described and discussed.

Multigrid for American option pricing with stochastic volatility

The paper describes an implicit finite difference approach to the pricing of American options on assets with a stochastic volatility. A multigrid procedure is described for the fast iterative solution of the discrete linear complementarity problems that result. The accuracy and performance of this approach is improved considerably by a strike-price related analytic transformation of asset prices and adaptive time-stepping.

Phase-accuracy comparisons and improved far-field estimates for 3-D edge elements on tetrahedral meshes

Edge-element methods have proved very effective for 3-D electromagnetic computations and are widely used on unstructured meshes. However, the accuracy of standard edge elements can be criticised because of their low order. This paper analyses discrete dispersion relations together with numerical propagation accuracy to determine the effect of tetrahedral shape on the phase accuracy of standard 3-D edgeelement approximations in comparison to other methods. Scattering computations for the sphere obtained with edge elements are compared with results obtained with vertex elements, and a new formulation of the far-field integral approximations for use with edge elements is shown to give improved cross sections over conventional formulations.

pKa values of aqua ligands of platinum (II) anticancer complexes: [1H, 15N] and 195Pt NMR Studies of cis- and trans- [PtCl2 (NH3) (Cyclohexylamine)]

Syntheses and NMR studies are reported of two 15N-labelled Pt(II) complexes of anticancer interest: cis-PtCl2(15NH3)(c-C6H1115NH2), a metabolite of the orally-active Pt(IV) complex cis,trans,cis-[PtCl2(acetate)2(c-C6H11NH2)(NH3), and trans-[PtCl2(15NH3)(c-C6H1115NH2), a reduction product of the active Pt(IV) complex trans,trans,trans-[PtCl2(OH)2(c-C6H11NH2). For cis-[PtCl2(15NH3)(c-C6H1115NH2), hydrolysis was faster for the chloride ligand trans to cyclohexylamine, and the pKa values determined by [1H, 15N NMR spectroscopy for the two cis monoaqua isomers were the same (6.73). The trans monoaqua complex was a stronger acid with pKa of 5.4 (determined by 195Pt NMR). For the cis diaqua complex, pKa values of 5.68 and 7.68 were determined.

The thermal desulfination of allylic sulfonyl halides

Allylic sulfonyl halides can be generated by halogenolysis of the corresponding triorganotin sulfinates. Allylic sulfonyl bromides and iodides undergo a first order, thermal desulfination with allylic rearrangement to yield the corresponding allylic halides. The desulfination of a cyclic allylic sulfonyl bromide is stereospecific, proceeding with T-syn bromine migration.

Building temporal constraints into knowledge bases for process control: an examination [Reply to author's reply]

Reply to author's reply

Conservative semi-Lagrangian finite volume schemes

Semi-Lagrangian finite volume schemes for the numerical approximation of linear advection equations are presented. These schemes are constructed so that the conservation properties are preserved by the numerical approximation. This is achieved using an interpolation procedure based on area-weighting. Numerical results are presented illustrating some of the features of these schemes.

A semi-Lagrangian finite volume method for Newtonian contraction flows

A new finite volume method for solving the incompressible Navier--Stokes equations is presented. The main features of this method are the location of the velocity components and pressure on different staggered grids and a semi-Lagrangian method for the treatment of convection. An interpolation procedure based on area-weighting is used for the convection part of the computation. The method is applied to flow through a constricted channel, and results are obtained for Reynolds numbers, based on half the flow rate, up to 1000. The behavior of the vortex in the salient corner is investigated qualitatively and quantitatively, and excellent agreement is found with the numerical results of Dennis and Smith [Proc. Roy. Soc. London A, 372 (1980), pp. 393-414] and the asymptotic theory of Smith [J. Fluid Mech., 90 (1979), pp. 725-754].

Modelling technology to predict flip-chip assembly

This paper describes modelling technology and its use in providing data governing the assembly of flip-chip components. Details are given on the reflow and curing stages as well as the prediction of solder joint shapes. The reflow process involves the attachment of a die to a board via solder joints. After a reflow process, underfill material is placed between the die and the substrate where it is heated and cured. Upon cooling the thermal mismatch between the die, underfill, solder bumps, and substrate will result in a nonuniform deformation profile across the assembly and hence stress. Shape predictions then thermal solidification and stress prediction are undertaken on solder joints during the reflow process. Both thermal and stress calculations are undertaken to predict phenomena occurring during the curing of the underfill material. These stresses may result in delamination between the underfill and its surrounding materials leading to a subsequent reduction in component performance and lifetime. Comparisons between simulations and experiments for die curvature will be given for the reflow and curing process