Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Linking physiological and architectural models of cotton

Despite the strong influence of plant architecture on crop yield, most crop models either ignore it or deal with it in a very rudimentary way. This paper demonstrates the feasibility of linking a model that simulates the morphogenesis and resultant architecture of individual cotton plants with a crop model that simulates the effects of environmental factors on critical physiological processes and resulting yield in cotton. First the varietal parameters of the models were made concordant. Then routines were developed to allocate the flower buds produced each day by the crop model amongst the potential positions generated by the architectural model. This allocation is done according to a set of heuristic rules. The final weight of individual bolls and the shedding of buds and fruit caused by water, N, and C stresses are processed in a similar manner. Observations of the positions of harvestable fruits, both within and between plants, made under a variety of agronomic conditions that had resulted in a broad range of plant architectures were compared to those predicted by the model with the same environmental inputs. As illustrated by comparisons of plant maps, the linked models performed reasonably well, though performance of the fruiting point allocation and shedding algorithms could probably be improved by further analysis of the spatial relationships of retained fruit. (C) 2002 Elsevier Science Ltd. All rights reserved.

Knowledge-based system on optimum design of liquid retaining structures with genetic algorithms

This paper delineates the development of a prototype hybrid knowledge-based system for the optimum design of liquid retaining structures by coupling the blackboard architecture, an expert system shell VISUAL RULE STUDIO and genetic algorithm (GA). Through custom-built interactive graphical user interfaces under a user-friendly environment, the user is directed throughout the design process, which includes preliminary design, load specification, model generation, finite element analysis, code compliance checking, and member sizing optimization. For structural optimization, GA is applied to the minimum cost design of structural systems with discrete reinforced concrete sections. The design of a typical example of the liquid retaining structure is illustrated. The results demonstrate extraordinarily converging speed as near-optimal solutions are acquired after merely exploration of a small portion of the search space. This system can act as a consultant to assist novice designers in the design of liquid retaining structures.

Electropalatographic assessment of tongue-to-palate contacts exhibited in dysarthria following traumatic brain injury: Spatial characteristics

Consonant imprecision has been reported to be a common feature of the dysarthric speech disturbances exhibited by individuals who have sustained a traumatic brain injury (TBI). Inaccurate tongue placements against the hard palate during consonant articulation may be one factor underlying the imprecision. To investigate this hypothesis, electropalatography (EPG) was used to assess the spatial characteristics of the tongue-to-palate contacts exhibited by three males (aged 23-29 years) with dysarthria following severe TBI. Five nonneurologically impaired adults served as control subjects. Twelve single-syllable words of CV or CVC construction (where initial C = /t, d, S, z, k, g/, V=/i, a/) were read aloud three times by each subject while wearing an EPG palate. Spatial characteristics were analyzed in terms of the location, pattern, and amount of tongue-to-palate contact at the frame of maximum contact during production of each consonant. The results revealed that for the majority of consonants, the patterns and locations of contacts exhibited by the TBI subjects were consistent with the contacts generated by the group of control subjects. One notable exception was one subject's production of the alveolar fricatives in which complete closure across the palate was demonstrated, rather than the characteristic groove configuration. Major discrepancies were also noted in relation to the amount of tongue-to-palate contact exhibited, with two TBI subjects consistently demonstrating increased contacts compared to the control subjects. The implications of these findings for the development of treatment programs for dysarthric speech disorders subsequent to TBI are highlighted.

Spatial variability in nutritional status of arabic coffee based on dris index

The combined use of precision agriculture and the Diagnosis and Recommendation Integrated System (DRIS) allows the spatial monitoring of coffee nutrient balance to provide more balanced and cost-effective fertilizer recommendations. The objective of this work was to evaluate the spatial variability in the nutritional status of two coffee varieties using the Mean Nutritional Balance Index (NBIm) and its relationship with their respective yields. The experiment was conducted in eastern Minas Gerais in two areas, one planted with variety Catucaí and another with variety Catuaí. The NBIm of the two varieties and their yields were analyzed through geostatistics and, based on the models and parameters of the variograms, were interpolated to obtain their spatial distribution in the studied areas. Variety Catucai, with grater spatial variability, was more nutritional unbalanced than variety Catuai, and consequently produced lower yields. Excess of Fe and Mn makes these elements limiting yield factors.

Spatial variability of particle size fractions of an Oxisol cultivated with conilon coffee

Information on the spatial distribution of particle size fractions is essential for use planning and management of soils. The aim of this work to was to study the spatial variability of particle size fractions of a Typic Hapludox cultivated with conilon coffee. The soil samples were collected at depths of 0-0.20 and 0.20-0.40 m in the coffee canopy projection, totaling 109 georeferentiated points. At the depth of 0.2-0.4 m the clay fraction showed average value significantly higher, while the sand fraction showed was higher in the depth of 0-0.20 m. The silt showed no significant difference between the two depths. The particle size fractions showed medium and high spatial variability. The levels of total sand and clay have positive and negative correlation, respectively, with the altitude of the sampling points, indicating the influence of landscape configuration.

Effect of spatial arrangement on the production components and yield of sunflower

ABSTRACT The sunflower plant is an oilseed crop that has aroused a great interest in the Brazilian and international scenery especillay because of the possibility of using its oil for biodiesel production. The objective of this study was to evaluate productivity and yield components of Embrapa 122 sunflower according to the spatial arrangement. Treatments were arranged in 4 x 4 factorial arrangement, which are the four spacings between rows (0.30; 0.50; 0.70 and 0.90 m) and four sowing densities (30,000; 45,000; 60,000 and 75,000 plants ha-1). The experiment was carried out in a complete randomized block design with four replications. The experiments were conducted in the experimental area of the Plant Science Department in Fortaleza, State of Ceará-Brazil and on the Curu Vale Experimental Farm in Pentecoste, State of Ceará-Brazil. Productivity and the following production components were analyzed in the end of the crop cycle: harvested capitula, capitulum diameter, capitulum mass, achene mass per capitulum, mass of 100 achenes, achenes per capitulum, harvest index and oil content in the achenes. The experiments were analyzed jointly in relation to the cropping area and the data submitted to analysis of variance and quantitative factors tested by polynomial regression. The variables, spacing, density and cropping area did not interact with these variables and the spatial arrangement of the plants affected only the components. The cropping area influences the productive behavior of sunflower Embrapa 122. The spatial arrangement of the plants of sunflower of variety Embrapa 122 influences yield components but does not affect productivity.

Surgical correction of scoliosis: Numerical analysis and optimization of the procedure

A previously developed model is used to numerically simulate real clinical cases of the surgical correction of scoliosis. This model consists of one-dimensional finite elements with spatial deformation in which (i) the column is represented by its axis; (ii) the vertebrae are assumed to be rigid; and (iii) the deformability of the column is concentrated in springs that connect the successive rigid elements. The metallic rods used for the surgical correction are modeled by beam elements with linear elastic behavior. To obtain the forces at the connections between the metallic rods and the vertebrae geometrically, non-linear finite element analyses are performed. The tightening sequence determines the magnitude of the forces applied to the patient column, and it is desirable to keep those forces as small as possible. In this study, a Genetic Algorithm optimization is applied to this model in order to determine the sequence that minimizes the corrective forces applied during the surgery. This amounts to find the optimal permutation of integers 1, ... , n, n being the number of vertebrae involved. As such, we are faced with a combinatorial optimization problem isomorph to the Traveling Salesman Problem. The fitness evaluation requires one computing intensive Finite Element Analysis per candidate solution and, thus, a parallel implementation of the Genetic Algorithm is developed.

New p-i-n Si : H imager configuration for spatial resolution improvement

Amorphous glass/ZnO-Al/p(a-Si:H)/i(a-Si:H)/n(a-Si1-xCx:H)/Al imagers with different n-layer resistivities were produced by plasma enhanced chemical vapour deposition technique (PE-CVD). An image is projected onto the sensing element and leads to spatially confined depletion regions that can be readout by scanning the photodiode with a low-power modulated laser beam. The essence of the scheme is the analog readout, and the absence of semiconductor arrays or electrode potential manipulations to transfer the information coming from the transducer. The influence of the intensity of the optical image projected onto the sensor surface is correlated with the sensor output characteristics (sensitivity, linearity blooming, resolution and signal-to-noise ratio) are analysed for different material compositions (0.5 < x < 1). The results show that the responsivity and the spatial resolution are limited by the conductivity of the doped layers. An enhancement of one order of magnitude in the image intensity signal and on the spatial resolution are achieved at 0.2 mW cm(-2) light flux by decreasing the n-layer conductivity by the same amount. A physical model supported by electrical simulation gives insight into the image-sensing technique used.

GA topology optimization using random keys for tree encoding of structures

Topology optimization consists in finding the spatial distribution of a given total volume of material for the resulting structure to have some optimal property, for instance, maximization of structural stiffness or maximization of the fundamental eigenfrequency. In this paper a Genetic Algorithm (GA) employing a representation method based on trees is developed to generate initial feasible individuals that remain feasible upon crossover and mutation and as such do not require any repairing operator to ensure feasibility. Several application examples are studied involving the topology optimization of structures where the objective functions is the maximization of the stiffness and the maximization of the first and the second eigenfrequencies of a plate, all cases having a prescribed material volume constraint.

Optimização dos protocolos de reconstrução em imagens TC de AVC isquémico agudo: projecto

Mestrado de Radiações aplicadas às Tecnologias da Saúde. Área de especialização: Imagem Digital com Radiação X.

Incidence rate and spatio-temporal clustering of type 1 diabetes in Santiago, Chile, from 1997 to 1998

OBJECTIVE: To estimate the incidence rate of type 1 diabetes in the urban area of Santiago, Chile, from March 21, 1997 to March 20, 1998, and to assess the spatio-temporal clustering of cases during that period. METHODS: All sixty-one incident cases were located temporally (day of diagnosis) and spatially (place of residence) in the area of study. Knox's method was used to assess spatio-temporal clustering of incident cases. RESULTS: The overall incidence rate of type 1 diabetes was 4.11 cases per 100,000 children aged less than 15 years per year (95% confidence interval: 3.06--5.14). The incidence rate seems to have increased since the last estimate of the incidence calculated for the years 1986--1992 in the metropolitan region of Santiago. Different combinations of space-time intervals have been evaluated to assess spatio-temporal clustering. The smallest p-value was found for the combination of critical distances of 750 meters and 60 days (uncorrected p-value = 0.048). CONCLUSIONS: Although these are preliminary results regarding space-time clustering in Santiago, exploratory analysis of the data method would suggest a possible aggregation of incident cases in space-time coordinates.