Subcycling first- and second-order generalizations of the trapezoidal rule

Subcycling algorithms which employ multiple timesteps have been previously proposed for explicit direct integration of first- and second-order systems of equations arising in finite element analysis, as well as for integration using explicit/implicit partitions of a model. The author has recently extended this work to implicit/implicit multi-timestep partitions of both first- and second-order systems. In this paper, improved algorithms for multi-timestep implicit integration are introduced, that overcome some weaknesses of those proposed previously. In particular, in the second-order case, improved stability is obtained. Some of the energy conservation properties of the Newmark family of algorithms are shown to be preserved in the new multi-timestep extensions of the Newmark method. In the first-order case, the generalized trapezoidal rule is extended to multiple timesteps, in a simple way that permits an implicit/implicit partition. Explicit special cases of the present algorithms exist. These are compared to algorithms proposed previously. (C) 1998 John Wiley & Sons, Ltd.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

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.

Explicit/implicit partitioning and a new explicit form of the generalized alpha method

Most finite element packages use the Newmark algorithm for time integration of structural dynamics. Various algorithms have been proposed to better optimize the high frequency dissipation of this algorithm. Hulbert and Chung proposed both implicit and explicit forms of the generalized alpha method. The algorithms optimize high frequency dissipation effectively, and despite recent work on algorithms that possess momentum conserving/energy dissipative properties in a non-linear context, the generalized alpha method remains an efficient way to solve many problems, especially with adaptive timestep control. However, the implicit and explicit algorithms use incompatible parameter sets and cannot be used together in a spatial partition, whereas this can be done for the Newmark algorithm, as Hughes and Liu demonstrated, and for the HHT-alpha algorithm developed from it. The present paper shows that the explicit generalized alpha method can be rewritten so that it becomes compatible with the implicit form. All four algorithmic parameters can be matched between the explicit and implicit forms. An element interface between implicit and explicit partitions can then be used, analogous to that devised by Hughes and Liu to extend the Newmark method. The stability of the explicit/implicit algorithm is examined in a linear context and found to exceed that of the explicit partition. The element partition is significantly less dissipative of intermediate frequencies than one using the HHT-alpha method. The explicit algorithm can also be rewritten so that the discrete equation of motion evaluates forces from displacements and velocities found at the predicted mid-point of a cycle. Copyright (C) 2003 John Wiley Sons, Ltd.

Analysis of transient eddy currents in MRI using a cylindrical FDTD method

Most magnetic resonance imaging (MRI) spatial encoding techniques employ low-frequency pulsed magnetic field gradients that undesirably induce multiexponentially decaying eddy currents in nearby conducting structures of the MRI system. The eddy currents degrade the switching performance of the gradient system, distort the MRI image, and introduce thermal loads in the cryostat vessel and superconducting MRI components. Heating of superconducting magnets due to induced eddy currents is particularly problematic as it offsets the superconducting operating point, which can cause a system quench. A numerical characterization of transient eddy current effects is vital for their compensation/control and further advancement of the MRI technology as a whole. However, transient eddy current calculations are particularly computationally intensive. In large-scale problems, such as gradient switching in MRI, conventional finite-element method (FEM)-based routines impose very large computational loads during generation/solving of the system equations. Therefore, other computational alternatives need to be explored. This paper outlines a three-dimensional finite-difference time-domain (FDTD) method in cylindrical coordinates for the modeling of low-frequency transient eddy currents in MRI, as an extension to the recently proposed time-harmonic scheme. The weakly coupled Maxwell's equations are adapted to the low-frequency regime by downscaling the speed of light constant, which permits the use of larger FDTD time steps while maintaining the validity of the Courant-Friedrich-Levy stability condition. The principal hypothesis of this work is that the modified FDTD routine can be employed to analyze pulsed-gradient-induced, transient eddy currents in superconducting MRI system models. The hypothesis is supported through a verification of the numerical scheme on a canonical problem and by analyzing undesired temporal eddy current effects such as the B-0-shift caused by actively shielded symmetric/asymmetric transverse x-gradient head and unshielded z-gradient whole-body coils operating in proximity to a superconducting MRI magnet.

Quantitative analysis of vascular to cardiac selectivity of L- and T-type voltage-operated calcium channel antagonists in human tissues

1. Classical L-type voltage-operated calcium channel (VOCC) antagonists dilate blood vessels, depress myocardial contractility and slow cardiac conduction. 2. We compared four L-type VOCC antagonists and a novel tetralol derivative, mibefradil, reportedly 10-fold more selective for T- (transient) over L-type VOCC in two in vitro assays of human tissue, namely isolated small arteries from the aortic vasa vasorum in a myograph and right atrial trabeculae muscle under isometric force conditions. 3. In arteries contracted with K+ (62 mmol/L), the relaxation pIC(50) values for the VOCC antagonists felodipine, nifedipine, amlodipine, verapamil and mibefradil were 8.30, 7.78, 6.64, 6.26 and 6.22, respectively. In atrial trabeculae, the pIC(50) values to inhibit the inotropic response to a submaximal concentration of isoprenaline (6 nmol/L) for felodipine, nifedipine, verapamil, amlodipine and mibefradil were 7.21, 6.95, 6.91, 5.94 and 4.61, respectively. 4. Taking the anti-log (pIC(50) vessel - pIC(50) atrium) the vascular relaxation to cardiac depression potency ratios for mibefradil, felodipine, nifedipine, amlodipine and verapamil were 41, 12, 7, 5 and 0.22, respectively. 5. We conclude that, in human tissue assays, perhaps T- over L-type VOCC selectivity confers the most favourable vascular selectivity on mibefradil. Alternatively, splice variants of L-type VOCC in the vasculature (CaV1.2b) may be more sensitive to mibefradil than the splice variants in the heart (CaV1.2a).

Transient expression of a novel serine protease in the ectoderm of the ascidian Herdmania momus during development

We have studied gene expression during ascidian embryonic development using the technique of differential display and isolated partial cDNA sequences of 12 genes. Developmental regulation of these genes has been confirmed by northern hybridization analysis. Further cDNA cloning and sequence analysis of an mRNA that is present during gastrulation, neurulation and tailbud formation reveals that it encodes a novel serine protease containing a single kringle motif and catalytic domain. The spatial expression of this gene, designated Hmserp1, is restricted to precursor cells of the epidermis. The structure and expression of Hmsery1 is discussed in relation to possible functions during development.

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.

Outcomes of transient evoked otoacoustic emission testing in 6-year-old school children: A comparison with pure tone screening and tympanometry

Objectives: (1) To establish test performance measures for Transient Evoked Otoacoustic Emission testing of 6-year-old children in a school setting; (2) To investigate whether Transient Evoked Otoacoustic Emission testing provides a more accurate and effective alternative to a pure tone screening plus tympanometry protocol. Methods: Pure tone screening, tympanometry and transient evoked otoacoustic emission data were collected from 940 subjects (1880 ears), with a mean age of 6.2 years. Subjects were tested in non-sound-treated rooms within 22 schools. Receiver operating characteristics curves along with specificity, sensitivity, accuracy and efficiency values were determined for a variety of transient evoked otoacoustic emission/pure tone screening/tympanometry comparisons. Results: The Transient Evoked Otoacoustic Emission failure rate for the group was 20.3%. The failure rate for pure tone screening was found to be 8.9%, whilst 18.6% of subjects failed a protocol consisting of combined pure tone screening and tympanometry results. In essence, findings from the comparison of overall Transient Evoked Otoacoustic Emission pass/fail with overall pure tone screening pass/fail suggested that use of a modified Rhode Island Hearing Assessment Project criterion would result in a very high probability that a child with a pass result has normal hearing (true negative). However, the hit rate was only moderate. Selection of a signal-to-noise ratio (SNR) criterion set at greater than or equal to 1 dB appeared to provide the best test performance measures for the range of SNR values investigated. Test performance measures generally declined when tympanometry results were included, with the exception of lower false alarm rates and higher positive predictive values. The exclusion of low frequency data from the Transient Evoked Otoacoustic Emission SNR versus pure tone screening analysis resulted in improved performance measures. Conclusions: The present study poses several implications for the clinical implementation of Transient Evoked Otoacoustic Emission screening for entry level school children. Transient Evoked Otoacoustic Emission pass/fail criteria will require revision. The findings of the current investigation offer support to the possible replacement of pure tone screening with Transient Evoked Otoacoustic Emission testing for 6-year-old children. However, they do not suggest the replacement of the pure tone screening plus tympanometry battery. (C) 2001 Elsevier Science Ireland Ltd. All rights reserved.

Handedness effects in transient evoked otoacoustic emissions in school children

Handedness, as a potentially influencing, nonpathologic factor, has not been investigated in relation to transient evoked otoacoustic emissions (TEOAEs). The present study aimed to examine the effects of handedness on the TEOAE spectrum in entry-level schoolchildren, with attention also to possible ear asymmetry. A total of 228 subjects (114 males, 114 females, mean age = 6.3 years) were tested using the ILO292 Otodynamics Analyzer (Quickscreen mode) in quiet rooms in 22 schools. For statistical analysis, subjects were matched for factors such as handedness, gender, age, and history of recent ear infection. The results from subjects with passing TEOAE, pure-tone screening, and tympanometry revealed no significant handedness effect overall, although a significant ear asymmetry effect on the measurement parameters of AB difference, noise level, response level, whole-wave reproducibility, band reproducibility, and signal-to-noise ratios was found.

Transient evoked otoacoustic emissions in adults: a comparison between two test protocols

This study compares the performance of the Quickscreen and Default protocols of the ILO-96 Otodynamics Analyzer in recording transient evoked otoacoustic emissions (TEOAEs) from adults using clinical decision analysis. Data were collected from 25 males (mean age = 29.0 years, SD = 6.8) and 35 females (mean age = 28.1 years, SD = 9.6). The results showed that the mean signal-to-noise ratios obtained from the Quickscreen were significantly greater than those from the Default protocol at 1,2, and 4 kHz. The comparison of the performance of the two protocols, based on the results using receiver operating characteristics curves, revealed a higher performance of the Quickscreen than the Default protocol at 1 and 4 kHz but not at 2 kHz. In view of the enhanced performance of the Quickscreen over the Default protocol in general, the routine use of the Default protocol for testing adults in audiology clinics should be reconsidered.

Coupled thermal, structural and vibrational analysis of a hypersonic engine for flight test

This paper describes a relatively simple and quick method for implementing aerodynamic heating models into a finite element code for non-linear transient thermal-structural and thermal-structural-vibrational analyses of a Mach 10 generic HyShot scramjet engine. The thermal-structural-vibrational response of the engine was studied for the descent trajectory from 60 to 26 km. Aerodynamic heating fluxes, as a function of spatial position and time for varying trajectory points, were implemented in the transient heat analysis. Additionally, the combined effect of varying dynamic pressure and thermal loads with altitude was considered. This aero-thermal-structural analysis capability was used to assess the temperature distribution, engine geometry distortion and yielding of the structural material due to aerodynamic heating during the descent trajectory, and for optimising the wall thickness, nose radius of leading edge, etc. of the engine intake. A structural vibration analysis was also performed following the aero-thermal-structural analysis to determine the changes in natural frequencies of the structural vibration modes that occur at the various temperatures associated with the descent trajectory. This analysis provides a unique and relatively simple design strategy for predicting and mitigating the thermal-structural-vibrational response of hypersonic engines. (C) 2006 Elsevier SAS. All rights reserved.

Identification and analysis of novel genes expressed in the mouse embryonic facial primordia

Craniofacial anomalies are a common feature of human congenital dysmorphology syndromes, suggesting that genes expressed in the developing face are likely to play a wider role in embryonic development. To facilitate the identification of genes involved in embryogenesis, we previously constructed an enriched cDNA library by subtracting adult mouse liver cDNA from that of embryonic day (E)10.5 mouse pharyngeal arch cDNA. From this library, 273 unique clones were sequenced and known proteins binned into functional categories in order to assess enrichment of the library (1). We have now selected 31 novel and poorly characterised genes from this library and present bioinformatic analysis to predict proteins encoded by these genes, and to detect evolutionary conservation. Of these genes 61% (19/31) showed restricted expression in the developing embryo, and a subset of these was chosen for further in silico characterisation as well as experimental determination of subcellular localisation based on transient transfection of predicted full-length coding sequences into mammalian cell lines. Where a human orthologue of these genes was detected, chromosomal localisation was determined relative to known loci for human congenital disease.