NMR shim coil design utilising a rapid spherical harmonic calculation method

A rapid spherical harmonic calculation method is used for the design of Nuclear Magnetic Resonance shim coils. The aim is to design each shim such that it generates a field described purely by a single spherical harmonic. By applying simulated annealing techniques, coil arrangements are produced through the optimal positioning of current-carrying circular arc conductors of rectangular cross-section. This involves minimizing the undesirable harmonies in relation to a target harmonic. The design method is flexible enough to be applied for the production of coil arrangements that generate fields consisting significantly of either zonal or tesseral harmonics. Results are presented for several coil designs which generate tesseral harmonics of degree one.

Large amplitude folding in finely layered viscoelastic rock structures

We analyze folding phenomena in finely layered viscoelastic rock. Fine is meant in the sense that the thickness of each layer is considerably smaller than characteristic structural dimensions. For this purpose we derive constitutive relations and apply a computational simulation scheme (a finite-element based particle advection scheme; see MORESI et al., 2001) suitable for problems involving very large deformations of layered viscous and viscoelastic rocks. An algorithm for the time integration of the governing equations as well as details of the finite-element implementation is also given. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered plate. The domain is compressed along the layer axis by prescribing velocities along the sides. First, for the viscous limit we consider the response to a series of harmonic perturbations of the director orientation. The Fourier spectra of the initial folding velocity are compared for different viscosity ratios. Turning to the nonlinear regime we analyze viscoelastic folding histories up to 40% shortening. The effect of layering manifests itself in that appreciable buckling instabilities are obtained at much lower viscosity ratios (1:10) as is required for the buckling of isotropic plates (1:500). The wavelength induced by the initial harmonic perturbation of the director orientation seems to be persistent. In the section of the parameter space considered here elasticity seems to delay or inhibit the occurrence of a second, larger wavelength. Finally, in a linear instability analysis we undertake a brief excursion into the potential role of couple stresses on the folding process. The linear instability analysis also provides insight into the expected modes of deformation at the onset of instability, and the different regimes of behavior one might expect to observe.

An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier B.V. All rights reserved.

Real exchange rate and elasticity of labour supply in a balance-of-payments-constrained macrodynamics

A macrodynamic model is proposed in which the real exchange rate and the elasticity of labour supply interact defining different trajectories of growth and income distribution in a developing economy. Growth depends on imports of capital goods which are paid with exports (there are no capital flows) and hence is constrained by equilibrium in current account. The role of the elasticity of labour supply is to prevent the real exchange rate from appreciating as the economy grows, thereby sustaining international competitiveness. The model allows for endogenous technological change and considers the impact of migration from the subsistence to the modern sector on the cumulative (Kaldor-Verdoorn) process of learning.

Structural change, balance-of-payments constraint, and economic growth: evidence from the multisectoral Thirlwall`s law

This paper contributes to the literature on balance-of-payments constrained growth by investigating how structural change identified with changes in the sectoral composition of exports and imports affects the external constraint We test both the original and a multisectoral version of Thirlwall`s law for a sample of Latin American and Asian countries The original Thirlwall s law is found to hold for all sample countries except South Korea, whereas the multisectoral analogue holds for all of them As the sectoral composition of exports and imports is found to matter for growth we analyze the evolution of each country`s weighted trade income elasticities

A spectral clustering algorithm for manufacturing cell formation

A graph clustering algorithm constructs groups of closely related parts and machines separately. After they are matched for the least intercell moves, a refining process runs on the initial cell formation to decrease the number of intercell moves. A simple modification of this main approach can deal with some practical constraints, such as the popular constraint of bounding the maximum number of machines in a cell. Our approach makes a big improvement in the computational time. More importantly, improvement is seen in the number of intercell moves when the computational results were compared with best known solutions from the literature. (C) 2009 Elsevier Ltd. All rights reserved.

A new algorithm and refined bounds for extended gcd computation

Extended gcd computation is interesting itself. It also plays a fundamental role in other calculations. We present a new algorithm for solving the extended gcd problem. This algorithm has a particularly simple description and is practical. It also provides refined bounds on the size of the multipliers obtained.

The Qu-Prolog unification algorithm: Formalisation and correctness

Qu-Prolog is an extension of Prolog which performs meta-level computations over object languages, such as predicate calculi and lambda-calculi, which have object-level variables, and quantifier or binding symbols creating local scopes for those variables. As in Prolog, the instantiable (meta-level) variables of Qu-Prolog range over object-level terms, and in addition other Qu-Prolog syntax denotes the various components of the object-level syntax, including object-level variables. Further, the meta-level operation of substitution into object-level terms is directly represented by appropriate Qu-Prolog syntax. Again as in Prolog, the driving mechanism in Qu-Prolog computation is a form of unification, but this is substantially more complex than for Prolog because of Qu-Prolog's greater generality, and especially because substitution operations are evaluated during unification. In this paper, the Qu-Prolog unification algorithm is specified, formalised and proved correct. Further, the analysis of the algorithm is carried out in a frame-work which straightforwardly allows the 'completeness' of the algorithm to be proved: though fully explicit answers to unification problems are not always provided, no information is lost in the unification process.

A 12-month randomized controlled trial of balance training in elderly women with ostebporosis: Improvement of quality of life

Objective: Physical and psychological incapacity, including fear of falling is related to decreased satisfaction with life in osteoporosis (OP). The impact of a balance exercise program on improving the quality of life is not well established. We have, therefore, investigated the effect of 12-month Balance Training Program in quality of life, functional balance and falls in elderly OP women. Methods: Sixty consecutive women with senile OP were randomized into a Balance Training Group (BT) of 30 patients and no intervention control group (CG) of 30 patients. The BT program included techniques to improve balance over a period of 12 months (1 h exercise session/week and home-based exercises). The quality of life was evaluated before and at the end of the trial using the Osteoporosis Assessment Questionnaire (OPAQ), functional balance was evaluated by Berg Balance Scale (BBS). Falls in the preceding year were noted and compared to the period of study. Results: The comparison of OPAQ variations (INITIAL-FINAL) revealed a significant improvement in quality of life in all parameters for BT compared to CG: well-being (1.61 +/- 1.44 vs. -1.46 +/- 1.32, p < 0001), physical function (1.30 +/- 1.33 vs. -0.36 +/- 0.82, p < 0.001), psychological status (1.58 +/- 1.36 vs. -1.02 +/- 0.83, p < 0.001), symptoms (2.76 +/- 1.96 vs. -0.63 +/- 0.87, p < 0.001), social interaction (1.01 +/- 1.51 vs. 0.35 +/- 1.08, p < 0.001). Of note, this overall benefit was paralleled by an improvement of BBS (-5.5 +/- 5.67 vs. +0.5 +/- 4.88 p < 0.001) and a reduction of falls in 50% in BT group vs. 26.6% for the CG (RR: 1.88, p < 0.025). Conclusion: The long-term Balance Training Program of OP women provides a striking overall health quality of life improvement in parallel with improving functional balance and reduced falls. (C) 2010 Elsevier Ireland Ltd. All rights reserved.

Analysis and implementation of a new constant acceleration subcycling algorithm

An algorithm for explicit integration of structural dynamics problems with multiple time steps is proposed that averages accelerations to obtain subcycle states at a nodal interface between regions integrated with different time steps. With integer time step ratios, the resulting subcycle updates at the interface sum to give the same effect as a central difference update over a major cycle. The algorithm is shown to have good accuracy, and stability properties in linear elastic analysis similar to those of constant velocity subcycling algorithms. The implementation of a generalised form of the algorithm with non-integer time step ratios is presented. (C) 1997 by John Wiley & Sons, Ltd.

Lack of galectin-3 alters the balance of innate immune cytokines and confers resistance to Rhodococcus equi infection

Galectin-3 is a p-galactoside-binding lectin implicated in the fine-tuning of innate immunity. Rhodococcus equi, a facultative intracellular bacterium of macrophages, causes severe granulomatous bronchopneumonia in young horses and immunocompromised humans. The aim of this study is to investigate the role of galectin-3 in the innate resistance mechanism against R. equi infection. The bacterial challenge of galectin-3-deficient mice (gal3(-/-)) and their wild-type counterpart (gal3(+/+)) revealed that the LD50 for the gal3(-/-) mice was about seven times higher than that for the gal3(+/+) mice. When challenged with a sublethal dose, gal3(-/-) mice showed lower bacteria counts and higher production of IL-12 and IFN-gamma production, besides exhibiting a delayed although increased inflammatory reaction. Gal3(-/-) macrophages exhibited a decreased frequency of bacterial replication and survival, and higher transcript levels of IL-1 beta, IL-6, IL-10, TLR2 and MyD88. R. equi-infected gal3(+/+) macrophages showed decreased expression of TLR2, whereas R. equi-infected gal3(-/-) macrophages showed enhanced expression of this receptor. Furthermore, galectin-3 deficiency in macrophages may be responsible for the higher IL-1 beta serum levels detected in infected gal3(-/-) mice. Therefore galectin-3 may exert a regulatory role in innate immunity by diminishing IL-1 beta production and thus affecting resistance to R. equi infection.

The subcycled Newmark algorithm

The popular Newmark algorithm, used for implicit direct integration of structural dynamics, is extended by means of a nodal partition to permit use of different timesteps in different regions of a structural model. The algorithm developed has as a special case an explicit-explicit subcycling algorithm previously reported by Belytschko, Yen and Mullen. That algorithm has been shown, in the absence of damping or other energy dissipation, to exhibit instability over narrow timestep ranges that become narrower as the number of degrees of freedom increases, making them unlikely to be encountered in practice. The present algorithm avoids such instabilities in the case of a one to two timestep ratio (two subcycles), achieving unconditional stability in an exponential sense for a linear problem. However, with three or more subcycles, the trapezoidal rule exhibits stability that becomes conditional, falling towards that of the central difference method as the number of subcycles increases. Instabilities over narrow timestep ranges, that become narrower as the model size increases, also appear with three or more subcycles. However by moving the partition between timesteps one row of elements into the region suitable for integration with the larger timestep these the unstable timestep ranges become extremely narrow, even in simple systems with a few degrees of freedom. As well, accuracy is improved. Use of a version of the Newmark algorithm that dissipates high frequencies minimises or eliminates these narrow bands of instability. Viscous damping is also shown to remove these instabilities, at the expense of having more effect on the low frequency response.

Simulated-annealing-based genetic algorithm for modeling the optical constants of solids

We propose a simulated-annealing-based genetic algorithm for solving model parameter estimation problems. The algorithm incorporates advantages of both genetic algorithms and simulated annealing. Tests on computer-generated synthetic data that closely resemble optical constants of a metal were performed to compare the efficiency of plain genetic algorithms against the simulated-annealing-based genetic algorithms. These tests assess the ability of the algorithms to and the global minimum and the accuracy of values obtained for model parameters. Finally, the algorithm with the best performance is used to fit the model dielectric function to data for platinum and aluminum. (C) 1997 Optical Society of America.