Transparent boundary conditions for wave propagation on unbounded domains

The numerical solution of the time dependent wave equation in an unbounded domain generally leads to a truncation of this domain, which requires the introduction of an artificial boundary with associated boundary conditions. Such nonreflecting conditions ensure the equivalence between the solution of the original problem in the unbounded region and the solution inside the artificial boundary. We consider the acoustic wave equation and derive exact transparent boundary conditions that are local in time and can be directly used in explicit methods. These conditions annihilate wave harmonics up to a given order on a spherical artificial boundary, and we show how to combine the derived boundary condition with a finite difference method. The analysis is complemented by a numerical example in two spatial dimensions that illustrates the usefulness and accuracy of transparent boundary conditions.

Improvement of compact terminal antenna performance by incorporating open-end slots in ground plane

This letter describes a new idea of increasing operational bandwidth of a compact planar inverted F antenna (PIFA) by introducing open-end slots in the ground plane under the radiating patch. The slots are not in the way of active modules of a wireless transceiver and thus the proposed antenna size reduction method is attractive from the point of view of practical implementation.

A distributed equivalent magnetic current based FDTD method for the calculation of E-fields induced by gradient coils

This paper evaluates a new, low-frequency finite-difference time-domain method applied to the problem of induced E-fields/eddy currents in the human body resulting from the pulsed magnetic field gradients in MRI. In this algorithm, a distributed equivalent magnetic current is proposed as the electromagnetic source and is obtained by quasistatic calculation of the empty coil's vector potential or measurements therein. This technique circumvents the discretization of complicated gradient coil geometries into a mesh of Yee cells, and thereby enables any type of gradient coil modelling or other complex low frequency sources. The proposed method has been verified against an example with an analytical solution. Results are presented showing the spatial distribution of gradient-induced electric fields in a multi-layered spherical phantom model and a complete body model. (C) 2004 Elsevier Inc. All rights reserved.

Comparative study of Li, Na, and K adsorptions on graphite by using ab initio method

A comprehensive study has been conducted to compare the adsorptions of alkali metals (including Li, Na, and K) on the basal plane of graphite by using molecular orbital theory calculations. All three metal atoms prefer to be adsorbed on the middle hollow site above a hexagonal aromatic ring. A novel phenomenon was observed, that is, Na, instead of Li or K, is the weakest among the three types of metal atoms in adsorption. The reason is that the SOMO (single occupied molecular orbital) of the Na atom is exactly at the middle point between the HOMO and the LUMO of the graphite layer in energy level. As a result, the SOMO of Na cannot form a stable interaction with either the HOMO or the LUMO of the graphite. On the other hand, the SOMO of Li and K can form a relatively stable interaction with either the HOMO or the LUMO of graphite. Why Li has a relatively stronger adsorption than K on graphite has also been interpreted on the basis of their molecular-orbital energy levels.

Distortional buckling of I-beams by finite element method

Distortional buckling, unlike the usual lateral-torsional buckling in which the cross-section remains rigid in its own plane, involves distortion of web in the cross-section. This type of buckling typically occurs in beams with slender web and stocky flanges. Most of the published studies assume the web to deform with a cubic shape function. As this assumption may limit the accuracy of the results, a fifth order polynomial is chosen here for the web displacements. The general line-type finite element model used here has two nodes and a maximum of twelve degrees of freedom per node. The model not only can predict the correct coupled mode but also is capable of handling the local buckling of the web.

Why H Atom Prefers the On-Top Site and Alkali Metals Favor the Middle Hollow Site on the Basal Plane of Graphite

In this work, the different adsorption properties of H and alkali metal atoms on the basal plane of graphite are studied and compared using a density functional method on the same model chemistry level. The results show that H prefers the on-top site while alkali metals favor the middle hollow site of graphite basal plane due to the unique electronic structures of H, alkali metals, and graphite. H has a higher electronegativity than carbon, preferring to form a covalent bond with C atoms, whereas alkaline metals have lower electronegativity, tending to adsorb on the highest electrostatic potential sites. During adsorption, there are more charges transferred from alkali metal to graphite than from H to graphite.

Modeling diffraction and imaging of laser beams by the mode-expansion method

We present a new method of modeling imaging of laser beams in the presence of diffraction. Our method is based on the concept of first orthogonally expanding the resultant diffraction field (that would have otherwise been obtained by the laborious application of the Huygens diffraction principle) and then representing it by an effective multimodal laser beam with different beam parameters. We show not only that the process of obtaining the new beam parameters is straightforward but also that it permits a different interpretation of the diffraction-caused focal shift in laser beams. All of the criteria that we have used to determine the minimum number of higher-order modes needed to accurately represent the diffraction field show that the mode-expansion method is numerically efficient. Finally, the characteristics of the mode-expansion method are such that it allows modeling of a vast array of diffraction problems, regardless of the characteristics of the incident laser beam, the diffracting element, or the observation plane. (C) 2005 Optical Society of America.

A wave equation technique for designing compact gradient coils

A primary purpose of this research is to design a gradient coil that is planar in construction and can be inserted within existing infrastructure. The proposed wave equation method for the design of gradient coils is novel within the field. it is comprehensively shown how this method can be used to design the planar x-, y-, and z-gradient wire windings to produce the required magnetic fields within a certain domain. The solution for the cylindrical gradient coil set is also elucidated. The wave equation technique is compared with the well-known target held method to gauge the quality of resultant design. In the case of the planar gradient coil design, it is shown that using the new method, a set of compact gradient coils with large field of view can be produced. The final design is considerably smaller in dimension when compared with the design obtained using the target field method, and therefore the manufacturing costs and materials required are somewhat reduced.

Numerical method for determination of the NMR frequency of the single-qubit operation in a silicon-based solid-state quantum computer

A numerical method is introduced to determine the nuclear magnetic resonance frequency of a donor (P-31) doped inside a silicon substrate under the influence of an applied electric field. This phosphorus donor has been suggested for operation as a qubit for the realization of a solid-state scalable quantum computer. The operation of the qubit is achieved by a combination of the rotation of the phosphorus nuclear spin through a globally applied magnetic field and the selection of the phosphorus nucleus through a locally applied electric field. To realize the selection function, it is required to know the relationship between the applied electric field and the change of the nuclear magnetic resonance frequency of phosphorus. In this study, based on the wave functions obtained by the effective-mass theory, we introduce an empirical correction factor to the wave functions at the donor nucleus. Using the corrected wave functions, we formulate a first-order perturbation theory for the perturbed system under the influence of an electric field. In order to calculate the potential distributions inside the silicon and the silicon dioxide layers due to the applied electric field, we use the multilayered Green's functions and solve an integral equation by the moment method. This enables us to consider more realistic, arbitrary shape, and three-dimensional qubit structures. With the calculation of the potential distributions, we have investigated the effects of the thicknesses of silicon and silicon dioxide layers, the relative position of the donor, and the applied electric field on the nuclear magnetic resonance frequency of the donor.

Linkage analyses of event-related potential slow wave phenotypes recorded in a working memory task

Working memory is an essential component of wide-ranging cognitive functions. It is a complex genetic trait probably influenced by numerous genes that individually have only a small influence. These genes may have an amplified influence on phenotypes closer to the gene action. In this study, event-related potential (ERP) phenotypes recorded during a working-memory task were collected from 656 adolescents from 299 families for whom genotypes were available. Univariate linkage analyses using the MERLIN variance-components method were conducted on slow wave phenotypes recorded at multiple sites while participants were required to remember the location of a target. Suggestive linkage (LOD > 2.2) was found on chromosomes 4, 5, 6, 10, 17, and 20. After correcting for multiple testing, suggestive linkage remained on chromosome 10. Empirical thresholds were computed for the most promising phenotypes. Those on chromosome 10 remained suggestive. A number of genes reported to regulate neural differentiation and function (i.e. NRP1, ANK3, and CHAT) were found under these linkage peaks and may influence the levels of neural activity occurring in individuals participating in a spatial working-memory task.

New insights into the interaction of hydrogen atoms with boron-substituted carbon

Boron substitution in carbon materials has been comprehensively investigated using the density functional theory method. It was found that there is a correlation between the stability of the graphene sheet, the distribution of T electrons, the electrostatic potential, and the capability for hydrogen-atom adsorption. Boron substitution destabilizes the graphene structure, increases the density of the electron wave around the substitutional boron atoms, and lowers the electrostatic potential, thus improving the hydrogen adsorption energy on carbon. However, this improvement is only ca. 10-20% instead of a factor of 4 or 5. Our calculations also show that two substitutional boron atoms provide consistent and reliable results, but one substitutional boron results in contradictory conclusions. This is a warning to other computational chemists who work on boron substitution that the conclusion from one substitutional boron might not be reliable.

Direct numerical simulation of turbulent flow and mass transfer around a rotating circular cylinder

Turbulent flow around a rotating circular cylinder has numerous applications including wall shear stress and mass-transfer measurement related to the corrosion studies. It is also of interest in the context of flow over convex surfaces where standard turbulence models perform poorly. The main purpose of this paper is to elucidate the basic turbulence mechanism around a rotating cylinder at low Reynolds numbers to provide a better understanding of flow fundamentals. Direct numerical simulation (DNS) has been performed in a reference frame rotating at constant angular velocity with the cylinder. The governing equations are discretized by using a finite-volume method. As for fully developed channel, pipe, and boundary layer flows, a laminar sublayer, buffer layer, and logarithmic outer region were observed. The level of mean velocity is lower in the buffer and outer regions but the logarithmic region still has a slope equal to the inverse of the von Karman constant. Instantaneous flow visualization revealed that the turbulence length scale typically decreases as the Reynolds number increases. Wavelet analysis provided some insight into the dependence of structural characteristics on wave number. The budget of the turbulent kinetic energy was computed and found to be similar to that in plane channel flow as well as in pipe and zero pressure gradient boundary layer flows. Coriolis effects show as an equivalent production for the azimuthal and radial velocity fluctuations leading to their ratio being lowered relative to similar nonrotating boundary layer flows.

Development of a non-intrusive system to monitor radial pulse wave velocity

Understanding arterial distensibility has shown to be important in the pathogenesis of cardiovascular abnormalities like hypertension. It is also known that arterial pulse wave velocity (PWV) is a measure of the elasticity or stiffness of peripheral arterial blood vessels. However, it generally requires complex instrumentations to have an accurate measurement and not suited for continual monitoring. In this paper, it describes a simple and non-intrusive method to detect the cardiovascular pulse from a human wrist above the radial artery and a fingertip. The main components of this proposed method are a piezoelectric transducer and a photo-plethysmography circuitry. 5 healthy adults (4 male) with age ranging from 25 to 38 years were recruited. The timing consistency of the detected pulsations is first evaluated and compared to that obtained from a commercial electrocardiogram. Furthermore, the derived PWV is then assessed by the predicted values attained from regression equations of two previous similar studies. The results show good correlations (p < 0.05) and similarities for the former and latter respectively. The simplicity and non-invasive nature of the proposed method can be attractive for even younger or badly disturbed patients. Moreover, it can be used for prolonged monitoring for the comfort of the patients.

Transition in homogeneously heated inclined plane parallel shear flows

The transition of internally heated inclined plane parallel shear flows is examined numerically for the case of finite values of the Prandtl number Pr. We show that as the strength of the homogeneously distributed heat source is increased the basic flow loses stability to two-dimensional perturbations of the transverse roll type in a Hopf bifurcation for the vertical orientation of the fluid layer, whereas perturbations of the longitudinal roll type are most dangerous for a wide range of the value of the angle of inclination. In the case of the horizontal inclination transverse roll and longitudinal roll perturbations share the responsibility for the prime instability. Following the linear stability analysis for the general inclination of the fluid layer our attention is focused on a numerical study of the finite amplitude secondary travelling-wave solutions (TW) that develop from the perturbations of the transverse roll type for the vertical inclination of the fluid layer. The stability of the secondary TW against three-dimensional perturbations is also examined and our study shows that for Pr=0.71 the secondary instability sets in as a quasi-periodic mode, while for Pr=7 it is phase-locked to the secondary TW. The present study complements and extends the recent study by Nagata and Generalis (2002) in the case of vertical inclination for Pr=0.

The vibrational spectra of rectangular plates

The wire drive pulse-echo system has been extensively used to excite and measure modes of vibration of thin rectangular plates. The frequency spectra of different modes have been investigated as a function of the material elastic moduli and the plate geometry. Most of the work was carried out on isotropic materials. For square plates a wide selection of materials were used. These were made isotropic in their in-plane dimensions where the displacements are taking place. The range of rnaterials enabled the dependence on Poisson's ratio to be investigated. A method of determining the value of Poisson's ratio resulted from this investigation. Certain modes are controlled principally by the shear modulus. Of these the fundamental has two nodal lines across the plate surface. One of them, which has nodes at the corners, (the Lame mode) is uniquely a pure shear mode where the diagonal is a full wave length. One controlled by the Young's modulus has been found. The precise harmonic relationship of the Lame mode series in square and rectangular plates was established. Use of the Rayleigh-Lamb equation has extended the theoretical support. The low order modes were followed over a wide range of sides ratios. Two fundamental types of modes have been recognised; These are the longitudinal modes where the frequency is controlled by the length of the plate only and the 2~f product has an asymptotic value approaching the rod velocity. The other type is the in-plane flexural modes (in effect a flexurally vibrating bar where the -2/w is the geometrical parameter). Where possible the experimental work was related to theory. Other modes controlled by the width dimension of the plate were followed. Anisotropic materials having rolled sheet elastic symmetry were investigated in terms of the appropriate theory. The work has been extended to examine materials from welds in steel plates.