Robust algorithms for solving stochastic partial differential equations

A robust semi-implicit central partial difference algorithm for the numerical solution of coupled stochastic parabolic partial differential equations (PDEs) is described. This can be used for calculating correlation functions of systems of interacting stochastic fields. Such field equations can arise in the description of Hamiltonian and open systems in the physics of nonlinear processes, and may include multiplicative noise sources. The algorithm can be used for studying the properties of nonlinear quantum or classical field theories. The general approach is outlined and applied to a specific example, namely the quantum statistical fluctuations of ultra-short optical pulses in chi((2)) parametric waveguides. This example uses a non-diagonal coherent state representation, and correctly predicts the sub-shot noise level spectral fluctuations observed in homodyne detection measurements. It is expected that the methods used wilt be applicable for higher-order correlation functions and other physical problems as well. A stochastic differencing technique for reducing sampling errors is also introduced. This involves solving nonlinear stochastic parabolic PDEs in combination with a reference process, which uses the Wigner representation in the example presented here. A computer implementation on MIMD parallel architectures is discussed. (C) 1997 Academic Press.

The formal representation of process system modelling assumptions and their implications

In order to analyse the effect of modelling assumptions in a formal, rigorous way, a syntax of modelling assumptions has been defined. The syntax of modelling assumptions enables us to represent modelling assumptions as transformations acting on the set of model equations. The notion of syntactical correctness and semantical consistency of sets of modelling assumptions is defined and methods for checking them are described. It is shown on a simple example how different modelling assumptions act on the model equations and their effect on the differential index of the resulted model is also indicated.

Applying linear time-varying constraints to econometric models: With an application to demand systems

When linear equality constraints are invariant through time they can be incorporated into estimation by restricted least squares. If, however, the constraints are time-varying, this standard methodology cannot be applied. In this paper we show how to incorporate linear time-varying constraints into the estimation of econometric models. The method involves the augmentation of the observation equation of a state-space model prior to estimation by the Kalman filter. Numerical optimisation routines are used for the estimation. A simple example drawn from demand analysis is used to illustrate the method and its application.

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.

The Sooty show laid bear - Children, puppets and make-believe

This article seeks to examine the complex nature of pretence as portrayed in a popular UK children's television puppet show. Central to the various theoretical concerns with make-believe philosophical, psychological and lucid - has been the phenomenological axis of the 'real-imaginary'. Because of the serial types of transformations undertaken by characters in the programme and the reliance placed on natural models of conversation, the problematic nature of this core boundary is highlighted. Despite its overt status as a fictional representation of animal behaviour, the very animality of the puppets is rendered opaque as their identities as children are linguistically accomplished. As a consequence we argue that the show as a piece of representational art is structured by moral and behavioural dictates typical of conventional adult-child interaction.

Assessment of factors that confound MRI and neuropathological correlation of human postmortem brain tissue

In spite of considerable technical advance in MRI techniques, the optical resolution of these methods are still limited. Consequently, the delineation of cytoarchitectonic fields based on probabilistic maps and brain volume changes, as well as small-scale changes seen in MRI scans need to be verified by neuronanatomical/neuropathological diagnostic tools. To attend the current interdisciplinary needs of the scientific community, brain banks have to broaden their scope in order to provide high quality tissue suitable for neuroimaging- neuropathology/anatomy correlation studies. The Brain Bank of the Brazilian Aging Brain Research Group (BBBABSG) of the University of Sao Paulo Medical School (USPMS) collaborates with researchers interested in neuroimaging-neuropathological correlation studies providing brains submitted to postmortem MRI in-situ. In this paper we describe and discuss the parameters established by the BBBABSG to select and to handle brains for fine-scale neuroimaging-neuropathological correlation studies, and to exclude inappropriate/unsuitable autopsy brains. We tried to assess the impact of the postmortem time and storage of the corpse on the quality of the MRI scans and to establish fixation protocols that are the most appropriate to these correlation studies. After investigation of a total of 36 brains, postmortem interval and low body temperature proved to be the main factors determining the quality of routine MRI protocols. Perfusion fixation of the brains after autopsy by mannitol 20% followed by formalin 20% was the best method for preserving the original brain shape and volume, and for allowing further routine and immunohistochemical staining. Taken to together, these parameters offer a methodological progress in screening and processing of human postmortem tissue in order to guarantee high quality material for unbiased correlation studies and to avoid expenditures by post-imaging analyses and histological processing of brain tissue.

On the validity of the dispersion model of hepatic drug elimination when intravascular transit time densities are long-tailed

The dispersion model with mixed boundary conditions uses a single parameter, the dispersion number, to describe the hepatic elimination of xenobiotics and endogenous substances. An implicit a priori assumption of the model is that the transit time density of intravascular indicators is approximated by an inverse Gaussian distribution. This approximation is limited in that the model poorly describes the tail part of the hepatic outflow curves of vascular indicators. A sum of two inverse Gaussian functions is proposed as ail alternative, more flexible empirical model for transit time densities of vascular references. This model suggests that a more accurate description of the tail portion of vascular reference curves yields an elimination rate constant (or intrinsic clearance) which is 40% less than predicted by the dispersion model with mixed boundary conditions. The results emphasize the need to accurately describe outflow curves in using them as a basis for determining pharmacokinetic parameters using hepatic elimination models. (C) 1997 Society for Mathematical Biology.

A two-dimensional fuzzy random model of soil pore structure

A new conceptual model for soil pore-solid structure is formalized. Soil pore-solid structure is proposed to comprise spatially abutting elements each with a value which is its membership to the fuzzy set ''pore,'' termed porosity. These values have a range between zero (all solid) and unity (all pore). Images are used to represent structures in which the elements are pixels and the value of each is a porosity. Two-dimensional random fields are generated by allocating each pixel a porosity by independently sampling a statistical distribution. These random fields are reorganized into other pore-solid structural types by selecting parent points which have a specified local region of influence. Pixels of larger or smaller porosity are aggregated about the parent points and within the region of interest by controlled swapping of pixels in the image. This creates local regions of homogeneity within the random field. This is similar to the process known as simulated annealing. The resulting structures are characterized using one-and two-dimensional variograms and functions describing their connectivity. A variety of examples of structures created by the model is presented and compared. Extension to three dimensions presents no theoretical difficulties and is currently under development.

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.

Evaluation of oral-motor movements and speech in patients with tetanus of a public service in Brazil

The characterisation of oral-motor movements and speech of patients with tetanus were investigated to determine the existence of possible signs that are characteristic of this pathology. Thirteen patients clinically diagnosed with tetanus (10 with severe tetanus and three with very severe tetanus) and admitted to an intensive care unit underwent clinical evaluation of oral-motor movements and speech. Statistical analysis indicated significant between-group differences for speech motor functions, suggesting that individuals with very severe tetanus present rigidity as a characteristic interfering in articulatory precision (P = 0 035) and movement rate (P = 0 038). For lip closure, tongue movement, palatal elevation, gag reflex and voice quality, no between-group differences were identified for the specific abnormal characteristics. The observed abnormal results indicate that muscle strength and functional status of the oral-motor system presented by most of the participants of the study did not ensure the necessary integrity for satisfactory performance. The characterisation of the oral myofunctional aspects of patients with tetanus provides medical teams, patients and families with a wider and better description of the clinical situation, giving support to the diagnosis, prognostics and treatment.

What needs to be done? Occupational therapy responsibilities and challenges regarding human rights

Background: Addressing human rights issues brings forth ethical and political responsibilities for occupational therapists and requires new epistemological and educational approaches. The way occupational therapists have faced these challenges has depended upon historical, cultural and social contexts. Aim and method: By means of literature review and historical analysis, this paper reflects on how occupational therapists have dealt with human matters issues and on the contemporary changes within the profession. Results and discussion: The paper portrays how Latin American occupational therapists have engaged in social transformation by choosing not to transform ethical and political problems into technical matters. Taking into account experiences and views from South Africa, Brazil and Chile, the paper outlines the importance of developing political literacy and interdisciplinary professional/postprofessional education to prepare the new generation of occupational therapists to engage in social transformation. Addressing issues of invisibility and lack of access to human rights, the paper reflects on the need of developing conceptual tools and strategies for change, and discusses the transformations being produced in contemporary occupational therapy. Conclusion: Occupational therapists and scientists need to be attentive to human rights issues. They also need to answer the call for interconnectedness in the present-day complex societies, and engage in networking and a cross-bordering dialogue. Nevertheless, although necessary and welcome, international cooperation requires a permanent exercise of cultural sensitivity, political awareness and self-awareness.

Virtual plants - integrating architectural and physiological models

Recent advances in computer technology have made it possible to create virtual plants by simulating the details of structural development of individual plants. Software has been developed that processes plant models expressed in a special purpose mini-language based on the Lindenmayer system formalism. These models can be extended from their architectural basis to capture plant physiology by integrating them with crop models, which estimate biomass production as a consequence of environmental inputs. Through this process, virtual plants will gain the ability to react to broad environmental conditions, while crop models will gain a visualisation component. This integration requires the resolution of the fundamentally different time scales underlying the approaches. Architectural models are usually based on physiological time; each time step encompasses the same amount of development in the plant, without regard to the passage of real time. In contrast, physiological models are based in real time; the amount of development in a time step is dependent on environmental conditions during the period. This paper provides a background on the plant modelling language, then describes how widely-used concepts of thermal time can be implemented to resolve these time scale differences. The process is illustrated using a case study. (C) 1997 Elsevier Science Ltd.

Model for Differential Nursing Diagnosis of Alterations in Urinary Elimination Based on Fuzzy Logic

Nursing diagnoses associated with alterations of urinary elimination require different interventions, Nurses, who are not specialists, require support to diagnose and manage patients with disturbances of urine elimination. The aim of this study was to present a model based on fuzzy logic for differential diagnosis of alterations in urinary elimination, considering nursing diagnosis approved by the North American Nursing Diagnosis Association, 2001-2002. Fuzzy relations and the maximum-minimum composition approach were used to develop the system. The model performance was evaluated with 195 cases from the database of a previous study, resulting in 79.0% of total concordance and 19.5% of partial concordance, when compared with the panel of experts. Total discordance was observed in only three cases (1.5%). The agreement between model and experts was excellent (kappa = 0.98, P < .0001) or substantial (kappa = 0.69, P < .0001) when considering the overestimative accordance (accordance was considered when at least one diagnosis was equal) and the underestimative discordance (discordance was considered when at least one diagnosis was different), respectively. The model herein presented showed good performance and a simple theoretical structure, therefore demanding few computational resources.

A Bayesian approach to fuzzy hypotheses testing for the estimation of optimal age for vaccination against measles

Fuzzy Bayesian tests were performed to evaluate whether the mother`s seroprevalence and children`s seroconversion to measles vaccine could be considered as ""high"" or ""low"". The results of the tests were aggregated into a fuzzy rule-based model structure, which would allow an expert to influence the model results. The linguistic model was developed considering four input variables. As the model output, we obtain the recommended age-specific vaccine coverage. The inputs of the fuzzy rules are fuzzy sets and the outputs are constant functions, performing the simplest Takagi-Sugeno-Kang model. This fuzzy approach is compared to a classical one, where the classical Bayes test was performed. Although the fuzzy and classical performances were similar, the fuzzy approach was more detailed and revealed important differences. In addition to taking into account subjective information in the form of fuzzy hypotheses it can be intuitively grasped by the decision maker. Finally, we show that the Bayesian test of fuzzy hypotheses is an interesting approach from the theoretical point of view, in the sense that it combines two complementary areas of investigation, normally seen as competitive. (C) 2007 IMACS. Published by Elsevier B.V. All rights reserved.

Assessment of Occupational Genotoxic Risk among Brazilian Hairdressers

Objectives: To evaluate the genotoxic risk to hairdressers exposed daily to chemical substances such as hair dyes, waving and straightening preparations and manicurists` products by the Comet assay test (single-cell gel electrophoresis). Methods: The Comet assay was performed on blood samples from 69 female hairdressers (36.4 +/- 10.7 years old) currently employed in 21 different beauty institutes in Sao Paulo, Brazil, and on 55 female control blood donors (32.6 +/- 10.0 years old) from the Sao Paulo University Clinical Hospital blood bank. All the control subjects had occupations other than hairdresser. Comet assays were performed by evaluating 100 blood lymphocytes per individual and graded by visual score according to comet tail length. Results: The hairdressers showed a higher frequency of DNA damage revealed by Comet Score (159.8 +/- 71) when compared to the control group (125.4 +/- 64.1), and the difference was statistically significant by the Student`s t-test (P = 0.005). Multiple regression analysis showed that in addition to the hairdressers` profession, tobacco use contributed to the higher frequency of cells with comets (P < 0.05). Conclusions: The observed DNA damage could be associated with the hairdressers` occupational environment, where different chemicals are chronically manipulated and inhaled. Considering that this profession in many countries, including Brazil, is not officially regulated, more attention should focus on these professionals not only by legislative bodies but also by multidisciplinary teams able to develop and implement risk prevention and control strategies for chemical, physical and biological agents to which hairdressers are exposed.