Numerical evaluation of e-fields induced by body motion near high-field MRI scanner

In modern magnetic resonance imaging (MRI), both patients and radiologists are exposed to strong, nonuniform static magnetic fields inside or outside of the scanner, in which the body movement may be able to induce electric currents in tissues which could be possibly harmful. This paper presents theoretical investigations into the spatial distribution of induced E-fields in the human model when moving at various positions around the magnet. The numerical calculations are based on an efficient, quasistatic, finite-difference scheme and an anatomically realistic, full-body, male model. 3D field profiles from an actively-shielded 4 T magnet system are used and the body model projected through the field profile with normalized velocity. The simulation shows that it is possible to induce E-fields/currents near the level of physiological significance under some circumstances and provides insight into the spatial characteristics of the induced fields. The results are easy to extrapolate to very high field strengths for the safety evaluation at a variety of field strengths and motion velocities.

Finite difference calculations of permeability in large domains in a wide porosity range.

Determining effective hydraulic, thermal, mechanical and electrical properties of porous materials by means of classical physical experiments is often time-consuming and expensive. Thus, accurate numerical calculations of material properties are of increasing interest in geophysical, manufacturing, bio-mechanical and environmental applications, among other fields. Characteristic material properties (e.g. intrinsic permeability, thermal conductivity and elastic moduli) depend on morphological details on the porescale such as shape and size of pores and pore throats or cracks. To obtain reliable predictions of these properties it is necessary to perform numerical analyses of sufficiently large unit cells. Such representative volume elements require optimized numerical simulation techniques. Current state-of-the-art simulation tools to calculate effective permeabilities of porous materials are based on various methods, e.g. lattice Boltzmann, finite volumes or explicit jump Stokes methods. All approaches still have limitations in the maximum size of the simulation domain. In response to these deficits of the well-established methods we propose an efficient and reliable numerical method which allows to calculate intrinsic permeabilities directly from voxel-based data obtained from 3D imaging techniques like X-ray microtomography. We present a modelling framework based on a parallel finite differences solver, allowing the calculation of large domains with relative low computing requirements (i.e. desktop computers). The presented method is validated in a diverse selection of materials, obtaining accurate results for a large range of porosities, wider than the ranges previously reported. Ongoing work includes the estimation of other effective properties of porous media.

Numerical and analytical study of an atherosclerosis inflammatory disease model

We study a reaction–diffusion mathematical model for the evolution of atherosclerosis as an inflammation process by combining analytical tools with computer-intensive numerical calculations. The computational work involved the calculation of more than sixty thousand solutions of the full reaction–diffusion system and lead to the complete characterisation of the ωω-limit for every initial condition. Qualitative properties of the solution are rigorously proved, some of them hinted at by the numerical study

Analysis of the damping properties of high-transmission-ratio Cycloidal drives

Cycloidal drives are widely used in today’s industries for drives where large reduction ratios are required. Drive-train dynamics plays an important role in their design. This paper presents a new methodology for assessing damping characteristics of Cycloidal drives and compares the natural frequencies obtained from experiments and theoretical/numerical calculations using Fast-Fourier-Transforms.

A comparison of stability and bifurcation criteria for inflated spherical elastic shells

The nonlinear stability analysis introduced by Chen and Haughton [1] is employed to study the full nonlinear stability of the non-homogeneous spherically symmetric deformation of an elastic thick-walled sphere. The shell is composed of an arbitrary homogeneous, incompressible elastic material. The stability criterion ultimately requires the solution of a third-order nonlinear ordinary differential equation. Numerical calculations performed for a wide variety of well-known incompressible materials are then compared with existing bifurcation results and are found to be identical. Further analysis and comparison between stability and bifurcation are conducted for the case of thin shells and we prove by direct calculation that the two criteria are identical for all modes and all materials.

Silicon on silicon : self-organized nanotip arrays formed in reactive Ar+H2 plasmas

The formation of arrays of vertically aligned nanotips on a moderately heated (up to 500 degrees C) Si surface exposed to reactive low-temperature radio frequency (RF) Ar+H(2) plasmas is studied. It is demonstrated that the nanotip surface density, aspect ratio and height dispersion strongly depend on the substrate temperature, discharge power, and gas composition. It is shown that nanotips with aspect ratios from 2.0 to 4.0 can only be produced at a higher RF power density (41.7 mW cm(-3)) and a hydrogen content of about 60%, and that larger aspect ratios can be achieved at substrate temperatures of about 300 degrees C. The use of higher (up to 500 degrees C) temperatures leads to a decrease of the aspect ratio but promotes the formation of more uniform arrays with the height dispersion decreasing to 1.5. At lower (approximately 20 mW cm(-3)) RF power density, only semispherical nanodots can be produced. Based on these experimental results, a nanotip formation scenario is proposed suggesting that sputtering, etching, hydrogen termination, and atom/radical re-deposition are the main concurrent mechanisms for the nanostructure formation. Numerical calculations of the ion flux distribution and hydrogen termination profiles can be used to predict the nanotip shapes and are in a good agreement with the experimental results. This approach can be applied to describe the kinetics of low-temperature formation of other nanoscale materials by plasma treatment.

Effect of elemental composition and size on electron confinement in self-assembled SiC quantum dots : a combinatorial approach

A high level of control over quantum dot (QD) properties such as size and composition during fabrication is required to precisely tune the eventual electronic properties of the QD. Nanoscale synthesis efforts and theoretical studies of electronic properties are traditionally treated quite separately. In this paper, a combinatorial approach has been taken to relate the process synthesis parameters and the electron confinement properties of the QDs. First, hybrid numerical calculations with different influx parameters for Si1-x Cx QDs were carried out to simulate the changes in carbon content x and size. Second, the ionization energy theory was applied to understand the electronic properties of Si1-x Cx QDs. Third, stoichiometric (x=0.5) silicon carbide QDs were grown by means of inductively coupled plasma-assisted rf magnetron sputtering. Finally, the effect of QD size and elemental composition were then incorporated in the ionization energy theory to explain the evolution of the Si1-x Cx photoluminescence spectra. These results are important for the development of deterministic synthesis approaches of self-assembled nanoscale quantum confinement structures.

Sol → gel transition with rings in a finite polycondensing system in solution

We report a theoretical formulation for the mean cluster size distribution in a finite polycondensing system. Expressions for the mean number of n-mers with j bonds ( nj) are developed. Numerical calculations show that while the non-cyclic molecules make the dominant contribution to the small clusters, the large clusters are dominated by cyclic structures. The number of particles in ringless chains, n n,n-1, decays monotonically with n at all extents of reaction, but n n becomes bimodal near the gel point. We also find that the solvent plays an important role in the cluster size distribution.

Some studies on free vibration of composite laminates

Free vibration analysis is carried out to study the vibration characteristics of composite laminates using the modified shear deformation, layered, composite plate theory and employing the Rayleigh-Ritz energy approach. The analysis is presented in a unified form so as to incorporate all different combinations of laminate boundary conditions and with full coverage with regard to the various design parameters of a laminated plate. A parametric study is made using a beam characteristic function as the admissible function for the numerical calculations. The numerical results presented here are for an example case of fully clamped boundary conditions and are compared with previously published results. The effect of parameters, such as the aspect ratio of plates, ply-angle, number of layers and also the thickness ratios of plies in laminates on the frequencies of the laminate, is systematically studied. It is found that for anti-symmetric angle-ply or cross-ply laminates unique numerical values of the thickness ratios exist which improve the vibration characteristics of such laminates. Numerical values of the non-dimensional frequencies and nodal patterns, using the thickness ratio distribution of the plies, are then obtained for clamped laminates, fabricated out of various commonly used composite materials, and are presented in the form of the design curves.

Quantum spin nematics, dimerization, and deconfined criticality in Quasi-1D spin-one magnets

We study theoretically the destruction of spin nematic order due to quantum fluctuations in quasi-one-dimensional spin-1 magnets. If the nematic ordering is disordered by condensing disclinations, then quantum Berry phase effects induce dimerization in the resulting paramagnet. We develop a theory for a Landau-forbidden second order transition between the spin nematic and dimerized states found in recent numerical calculations. Numerical tests of the theory are suggested.

Theoretical and computational approaches on heterogeneous nucleation

Nucleation is the first step of a first order phase transition. A new phase is always sprung up in nucleation phenomena. The two main categories of nucleation are homogeneous nucleation, where the new phase is formed in a uniform substance, and heterogeneous nucleation, when nucleation occurs on a pre-existing surface. In this thesis the main attention is paid on heterogeneous nucleation. This thesis wields the nucleation phenomena from two theoretical perspectives: the classical nucleation theory and the statistical mechanical approach. The formulation of the classical nucleation theory relies on equilibrium thermodynamics and use of macroscopically determined quantities to describe the properties of small nuclei, sometimes consisting of just a few molecules. The statistical mechanical approach is based on interactions between single molecules, and does not bear the same assumptions as the classical theory. This work gathers up the present theoretical knowledge of heterogeneous nucleation and utilizes it in computational model studies. A new exact molecular approach on heterogeneous nucleation was introduced and tested by Monte Carlo simulations. The results obtained from the molecular simulations were interpreted by means of the concepts of the classical nucleation theory. Numerical calculations were carried out for a variety of substances nucleating on different substances. The classical theory of heterogeneous nucleation was employed in calculations of one-component nucleation of water on newsprint paper, Teflon and cellulose film, and binary nucleation of water-n-propanol and water-sulphuric acid mixtures on silver nanoparticles. The results were compared with experimental results. The molecular simulation studies involved homogeneous nucleation of argon and heterogeneous nucleation of argon on a planar platinum surface. It was found out that the use of a microscopical contact angle as a fitting parameter in calculations based on the classical theory of heterogeneous nucleation leads to a fair agreement between the theoretical predictions and experimental results. In the presented cases the microscopical angle was found to be always smaller than the contact angle obtained from macroscopical measurements. Furthermore, molecular Monte Carlo simulations revealed that the concept of the geometrical contact parameter in heterogeneous nucleation calculations can work surprisingly well even for very small clusters.

Vibration of simply-supported trapezoidal plates I. Symmetric trapezoids

A detailed investigation of the natural frequencies and mode shapes of simply supported symmetric trapezoidal plates is undertaken in this paper. For numerical calculations, the relationship that exists between the eigenvalue problem of a polygonal simply supported plate and the eigenvalue problem of polygonal membrane of the same shape is utilized with advantage. The deflection surface is expressed in terms of a Fourier sine series in transformed coordinates and the Galerkin method is used. Results are presented in the form of tables and graphs. Several features like the crossing of frequency curves and the metamorphosis of some of the nodal patterns are observed. By a suitable interpretation of the modes of those symmetric trapezoidal plates which have the median as the nodal line, the results for some of the modes of unsymmetrical trapezoidal plates are also deduced.

Vibration of simply-supported trapezoidal plates II. Unsymmetric trapezoids

This paper deals with the investigation of the vibration characteristics of simply-supported unsymmetric trapezoidal plates. For numerical calculations, the relationship between the eigenvalue problems of a polygonal simply-supported plate and polygonal membrane is again effectively utilized. The Galerkin method is applied, with the deflection surface expressed in terms of a Fourier sine series in transformed coordinates. Numerical values for the first seven to eight frequencies for different geometries of the unsymmetric trapezoid are presented in the form of tables. Also the nodal patterns for a few representative configurations are presented.

Penetration of hydromagnetic energy deep into the magnetosphere

Magnetometer data, acquired on spacecraft and simultaneously at high and low latitudes on the ground, are compared in order to study the propagation characteristics of hydromagnetic energy deep into the magnetosphere. Single events provide evidence that wave energy at L ∼ 3 can at times be only one order of magnitude lower than at L ∼ 13. In addition, statistical analyses of the H-component groundbased data obtained during local daytime hours of 17 July-3 August 1985 show that wave amplitudes at L ∼ 3 are generally 10-30 times lower than at L ∼ 13. The L-dependence of near-equator magnetic field fluctuations measured on ISEE-2 show a sharp drop in energy near the magnetopause and a more gradual fall-off of energy deeper inside the magnetosphere. Such high levels of wave power deep in the magnetosphere have not been quantitatively understood previously. Our initial attempt is to calculate the decay length of an evanescent wave generated at a thick magnetopause boundary. Numerical calculations show that fast magnetosonic modes (called magnetopause and inner mode) can be generated under very restrictive conditions for the field and plasma parameters. These fast compressional modes may have their energy reduced by only one order of magnitude over a penetration depth of about 8RE. More realistic numerical simulations need to be carried out to see whether better agreement with the data can be attained.

Theoretical Treatment of Interaction Parameters in Multicomponent Metallic Solutions

A new theoretical equation for interaction parameter in multicomponent metallic solutions is developed using the pseudopotential formalism coupled with the free energy of the hard sphere system. The approximate expression for the pseudopotential term is given in terms of the heat of solution at infinite dilution, to allow easy evaluation of the interaction parameter in various multicomponent systems. This theory has been applied to 23 non-ferrous alloys based on Pb, Sn, Bi and indium. Comparison with the results of previous theoretical calculations using only the hard sphere model suggests that the inclusion of the pseudopotential term yields a quantitatively more correct prediction of interaction parameters in multicomponent metallic solutions. Numerical calculations were also made for 320 Fe-base solutions relevant to steelmaking and the agreement between calculation and experimental data appears reasonable, with 90% reliability in predicting the correct sign.