Robust model-order reduction of complex biological processes

This paper addresses robust model-order reduction of a high dimensional nonlinear partial differential equation (PDE) model of a complex biological process. Based on a nonlinear, distributed parameter model of the same process which was validated against experimental data of an existing, pilot-scale BNR activated sludge plant, we developed a state-space model with 154 state variables in this work. A general algorithm for robustly reducing the nonlinear PDE model is presented and based on an investigation of five state-of-the-art model-order reduction techniques, we are able to reduce the original model to a model with only 30 states without incurring pronounced modelling errors. The Singular perturbation approximation balanced truncating technique is found to give the lowest modelling errors in low frequency ranges and hence is deemed most suitable for controller design and other real-time applications. (C) 2002 Elsevier Science Ltd. All rights reserved.

Finite difference time domain (FDTD) method for modeling the effect of switched gradients on the human body in MRI

Numerical modeling of the eddy currents induced in the human body by the pulsed field gradients in MRI presents a difficult computational problem. It requires an efficient and accurate computational method for high spatial resolution analyses with a relatively low input frequency. In this article, a new technique is described which allows the finite difference time domain (FDTD) method to be efficiently applied over a very large frequency range, including low frequencies. This is not the case in conventional FDTD-based methods. A method of implementing streamline gradients in FDTD is presented, as well as comparative analyses which show that the correct source injection in the FDTD simulation plays a crucial rule in obtaining accurate solutions. In particular, making use of the derivative of the input source waveform is shown to provide distinct benefits in accuracy over direct source injection. In the method, no alterations to the properties of either the source or the transmission media are required. The method is essentially frequency independent and the source injection method has been verified against examples with analytical solutions. Results are presented showing the spatial distribution of gradient-induced electric fields and eddy currents in a complete body model.

Does a reading lexicon provide orthographic representations for spelling?

University students spelled low-frequency words to dictation and subsequently made lexical decisions to them. In Experiment I, lexical decisions were slower on words students had spelled incorrectly relative to words they had spelled correctly, and there A as a larger repetition benefit 101 incorrectly spelled words. In experiment 2, the latency advantage for items spelled correctly was replicated when words were presented for only 200 ms and also in a spelling recognition task, In Experiment 3. masked identity and form priming effects were similar for words that had been spelled correctly and incorrectly, Item spelling accuracy tracked word frequency effects in the way chat it combined with repetition and priming effects. we inter that an individuals learning with a word's orthography underlies word frequency and item spelling accuracy effects and that a single orthographic lexicon serves visual word recognition and spelling. (C) 2000 Elsevier Science (USA).

An investigation of magnetic antennas for ground penetrating radar

For ground penetrating radar (GPR), smaller antennas would provide considerable practical advantages. Some of which are: portability; ease of use; and higher spatial sampling. A theoretical comparison of the fundamental limits of a small electric field antenna and a small magnetic field antenna shows that the minimum Q constraints are identical. Furthermore, it is shown that only the small magnetic loop antenna can be constructed to approach, arbitrarily closely, the fundamental minimum Q limit. This is achieved with the addition of a high permeability material which reduces energy stored in the magnetic fields. This is of special interest to some GPR applications. For example, applications requiring synthetic aperture data collection would benefit from the increased spatial sampling offered by electrically smaller antennas. Low frequency applications may also benefit, in terms of reduced antenna dimensions, by the use of electrically small antennas. Under these circumstances, a magnetic type antenna should be considered in preference to the typical electric field antenna. Numerical modeling data supports this assertion.

Experimental manipulation of prior experience: Effects on item and associative recognition

Frequency of exposure to very low- and high-frequency words was manipulated in a three-phase (familiarisation, study, and test) design. During familiarisation, words were presented with their definition (once, four times, or not presented). One week (Experiment 1) or one day (Experiment 2) later, participants studied a list of homogeneous pairs (i.e., pair members were matched on background and familiarisation frequency). Item and associative recognition of high- and very low-frequency words presented in intact, rearranged, old-new, or new-new pairs were tested in Experiment 1. Associative recognition of very low-frequency words was tested in Experiment 2. Results showed that prior familiaris ation improved associative recognition of very low-frequency pairs, but had no effect on high-frequency pairs. The role of meaning in the formation of item-to-item and item-to-context associations and the implications for current models of memory are discussed.

Layered silicate nanocomposites based on various high-functionality epoxy resins. Part II: The influence of an organoclay on the rheological behavior of epoxy prepolymers

The effect of an organically surface modified layered silicate on the viscosity of various epoxy resins of different structures and different functionalities was investigated. Steady and dynamic shear viscosities of the epoxy resins containing 0-10 wt% of the organoclay were determined using parallel plate rheology. Viscosity results were compared with those achieved through addition of a commonly used micron-sized CaCO3 filler. It was found that changes in viscosities due to the different fillers were of the same order, since the layered silicate was only dispersed on a micron-sized scale in the monomer (prior to reaction), as indicated by X-ray diffraction measurements. Flow activation energies at a low frequency were determined and did not show any significant changes due to the addition of organoclay or CaCO3. Comparison between dynamic and steady shear experiments showed good agreement for low layered silicate concentrations below 7.5 wt%, i.e. the Cox-Merz rule can be applied. Deviations from the Cox-Merz rule appeared at and above 10 wt%, although such deviations were only slightly above experimental error. Most resin organoclay blends were well predicted by the Power Law model, only concentrations of 10 wt% and above requiring the Herschel-Buckley (yield stress) model to achieve better fits. Wide-angle X-ray measurements have shown that the epoxy resin swells the layered silicate with an increase in the interlayer distance of approximately 15 Angstrom, and that the rheology behavior is due to the lateral, micron-size of these swollen tactoids.

Reflection of structural waves at a solid/liquid interface

This paper investigates the reflection characteristics of structural or guided waves in rods at a solid/liquid interface. Structural waves, whose wavelengths are much larger than the diameter of the rod, are described in a first approximation by classical one-dimensional wave theory. The reflection characteristics of such waves at a solid/liquid (melting) interface has been reported by two different ultrasonic measurement techniques: first, measuring the fast regression rate of a melting interface during the burning of metal rod samples in an oxygen-enriched environment, and second, monitoring the propagation of the solid/liquid interface during the slow melting and solidification of a rod sample in a furnace. The second work clearly shows that the major reflection occurs from the solid/liquid interface and not the liquid/gas interface as predicted by plane longitudinal wave reflectivity theory. The present work confirms this observation by reporting on the results of some specially designed experiments to identify the main interface of reflection for structural waves in rods. Hence, it helps in explaining the fundamental discrepancy between the reflection characteristics at a solid/liquid interface between low frequency structural waves and high frequency bulk waves, and confirms that the detected echo within a burning metallic rod clearly represents a reflection from the solid/liquid interface. (C) 2003 Elsevier Science B.V. All rights reserved.

Isolation and characterization of microsatellite loci from Helicoverpa armigera Hübner (Lepidoptera: Noctuidae)

Five microsatellite loci are presented for Helicoverpa armigera. These microsatellite loci were obtained through the construction of enriched libraries, overcoming previously reported difficulties with obtaining microsatellites from H. armigera and other Lepidoptera due to the low frequency of microsatellites in their genomes. The description of a further five microsatellite loci for H. armigera makes microsatellite based population genetics studies feasible.

Life history parameters and biocontrol potential of the mealybug parasitoid Coccidoxenoides peregrinus (Timberlake) (Hymenoptera: Encyrtidae): Asexuality, fecundity and ovipositional patterns

Properties relevant to the ovipositional activity and lifetime productivity of Coccidoxenoides peregrinus (Timberlake) were assessed in the laboratory, to determine the potential of this species as a biocontrol agent against the citrus mealybug, Planococcus citri (Risso). In general, this species has not performed well in orchards, except for a few localities on different continents. The mode of reproduction of C peregrinus is almost entirely thelytokous, with males produced sporadically and at low frequency. The females have both pro-ovigenic and synovigenic traits, which raises questions of the utility of this distinction. The females have a high reproductive potential with 10-20 eggs per day available within the first two days (after a short (12 h) pre-oviposition period), and 80-150 eggs per day thereafter until death at about eight days. Mean lifetime fecundity was 239.2 +/- 34.3 eggs. C peregrinus oviposits across a range of P. citri instars, but productivity relies predominantly on second instar hosts. Second stage (N2) hosts received most eggs in choice (about 52%) and no-choice (about 50%) tests. Most eggs deposited into N2 hosts (82%) reached adult stage whereas only a few of those deposited into N1 and N3 (about 5% each) developed successfully. The haemolymph of parasitised reproductive mealybugs contained granular structures and no parasitoid eggs were found 24 h after exposure to ovipositing wasps. Also, no wasps emerged from parasitised adult hosts that were kept alive. Parasitoid eggs deposited into adult hosts were presumed encapsulated and destroyed, as control mealybugs (not exposed to female wasps) had no granular structures in their haemolymph. Wasps exposed to an abundance of hosts soon started ovipositing, but only for a relatively short time each day (about 2.5 h out of a 7 h exposure). They stopped ovipositing despite eggs judged to be mature in their ovaries. The reproductive output of C peregrinus is discussed in relation to the ecological factors that could influence this output, and the implications for biocontrol are discussed. (C) 2003 Elsevier Inc. All rights reserved.

Isolation and characterization of microsatellite loci from Chiasmia assimilis (Warren, 1899) (Lepidoptera: Geometridae)

Twelve microsatellite loci are presented for the biological control agent Chiasmia assimilis (Warren, 1899). These microsatellite loci were obtained through the construction of an enriched library, overcoming previous reported difficulties with obtaining microsatellites from other Lepidoptera due to the low frequency of microsatellites in their genomes.

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.

An FDTD model for calculation of gradient-induced eddy currents in MRI system

In most magnetic resonance imaging (MRI) systems, pulsed magnetic gradient fields induce eddy currents in the conducting structures of the superconducting magnet. The eddy currents induced in structures within the cryostat are particularly problematic as they are characterized by long time constants by virtue of the low resistivity of the conductors. This paper presents a three-dimensional (3-D) finite-difference time-domain (FDTD) scheme in cylindrical coordinates for eddy-current calculation in conductors. This model is intended to be part of a complete FDTD model of an MRI system including all RF and low-frequency field generating units and electrical models of the patient. The singularity apparent in the governing equations is removed by using a series expansion method and the conductor-air boundary condition is handled using a variant of the surface impedance concept. The numerical difficulty due to the asymmetry of Maxwell equations for low-frequency eddy-current problems is circumvented by taking advantage of the known penetration behavior of the eddy-current fields. A perfectly matched layer absorbing boundary condition in 3-D cylindrical coordinates is also incorporated. The numerical method has been verified against analytical solutions for simple cases. Finally, the algorithm is illustrated by modeling a pulsed field gradient coil system within an MRI magnet system. The results demonstrate that the proposed FDTD scheme can be used to calculate large-scale eddy-current problems in materials with high conductivity at low frequencies.

The relationship between the auditory brain-stem response and its reconstructed waveforms following discrete wavelet transformation

Objective: To examine the relationship between the auditory brain-stem response (ABR) and its reconstructed waveforms following discrete wavelet transformation (DWT), and to comment on the resulting implications for ABR DWT time-frequency analysis. Methods: ABR waveforms were recorded from 120 normal hearing subjects at 90, 70, 50, 30, 10 and 0 dBnHL, decomposed using a 6 level discrete wavelet transformation (DWT), and reconstructed at individual wavelet scales (frequency ranges) A6, D6, D5 and D4. These waveforms were then compared for general correlations, and for patterns of change due to stimulus level, and subject age, gender and test ear. Results: The reconstructed ABR DWT waveforms showed 3 primary components: a large-amplitude waveform in the low-frequency A6 scale (0-266.6 Hz) with its single peak corresponding in latency with ABR waves III and V; a mid-amplitude waveform in the mid-frequency D6 scale (266.6-533.3 Hz) with its first 5 waves corresponding in latency to ABR waves 1, 111, V, VI and VII; and a small-amplitude, multiple-peaked waveform in the high-frequency D5 scale (533.3-1066.6 Hz) with its first 7 waves corresponding in latency to ABR waves 1, 11, 111, IV, V, VI and VII. Comparisons between ABR waves 1, 111 and V and their corresponding reconstructed ABR DWT waves showed strong correlations and similar, reliable, and statistically robust changes due to stimulus level and subject age, gender and test ear groupings. Limiting these findings, however, was the unexplained absence of a small number (2%, or 117/6720) of reconstructed ABR DWT waves, despite their corresponding ABR waves being present. Conclusions: Reconstructed ABR DWT waveforms can be used as valid time-frequency representations of the normal ABR, but with some limitations. In particular, the unexplained absence of a small number of reconstructed ABR DWT waves in some subjects, probably resulting from 'shift invariance' inherent to the DWT process, needs to be addressed. Significance: This is the first report of the relationship between the ABR and its reconstructed ABR DWT waveforms in a large normative sample. (C) 2004 International Federation of Clinical Neurophysiology. Published by Elsevier Ireland Ltd. All rights reserved.

Risk of non-small cell lung cancer and the cytochrome P4501A1 Ile462Val polymorphism

Objective: The Ile462Val substitution in the cytochrome P450 1A1 gene (CYP1A1) results in increased enzymatic activity. Preliminary data suggesting a link between this polymorphism and lung cancer risk in Caucasians are inconsistent, reflecting small sample sizes and the relatively low frequency of the variant. Methods: The data set consisted of 1050 primary non-small cell lung cancer cases and 581 controls, a large homogenous population designed specifically to address previous inconsistencies. Patients were genotyped using a PCR-RFLP technique. Results: Carriers of the valine allele, CYP1A1*2C, (Ile/Val or Val/Val genotypes) were significantly over-represented in non-small cell lung cancer compared to controls (OR=1.9; 95% CI=1.2-2.9; p=0.005) when adjusted for confounders, particularly in women (OR=4.6; 95% CI=1.7-12.4; p=0.003). The valine variant was statistically significantly over-represented in cases of lung cancer younger than the median age (64 years) (OR=2.5; 95% CI=1.3-4.8; p=0.005) and cases with less than the median cumulative tobacco-smoke exposure (46 pack-years) (OR=2.4; 95% CI=1.3-4.7; p=0.007). Conclusions: These new data establish an association between the CYP1A1 Ile462Val polymorphism and the risk of developing non-small cell lung cancer, especially among women.

An object-oriented designed finite-difference time-domain simulator for electromagnetic analysis and design in MRI - application to high field analyses

This paper presents a finite-difference time-domain (FDTD) simulator for electromagnetic analysis and design applications in MRI. It is intended to be a complete FDTD model of an MRI system including all RF and low-frequency field generating units and electrical models of the patient. The pro-ram has been constructed in an object-oriented framework. The design procedure is detailed and the numerical solver has been verified against analytical solutions for simple cases and also applied to various field calculation problems. In particular, the simulator is demonstrated for inverse RF coil design, optimized source profile generation, and parallel imaging in high-frequency situations. The examples show new developments enabled by the simulator and demonstrate that the proposed FDTD framework can be used to analyze large-scale computational electromagnetic problems in modern MRI engineering. (C) 2004 Elsevier Inc. All rights reserved.