A dynamic model with on-line identification for rotary sugar drying processes

A dynamic modelling methodology, which combines on-line variable estimation and parameter identification with physical laws to form an adaptive model for rotary sugar drying processes, is developed in this paper. In contrast to the conventional rate-based models using empirical transfer coefficients, the heat and mass transfer rates are estimated by using on-line measurements in the new model. Furthermore, a set of improved sectional solid transport equations with localized parameters is developed in this work to reidentified on-line using measurement data, the model is able to closely track the dynamic behaviour of rotary drying processes within a broad range of operational conditions. This adaptive model is validated against experimental data obtained from a pilot-scale rotary sugar dryer. The proposed modelling methodology can be easily incorporated into nonlinear model based control schemes to form a unified modelling and control framework.place the global correlation for the computation of solid retention time. Since a number of key model variables and parameters are identified on-line using measurement data, the model is able to closely track the dynamic behaviour of rotary drying processes within a broad range of operational conditions. This adaptive model is validated against experimental data obtained from a pilot-scale rotary sugar dryer. The proposed modelling methodology can be easily incorporated into nonlinear model based control schemes to form a unified modelling and control framework.

Design, modelling and analysis of a six component force balance for hypervelocity wind tunnel testing

A combination of modelling and analysis techniques was used to design a six component force balance. The balance was designed specifically for the measurement of impulsive aerodynamic forces and moments characteristic of hypervelocity shock tunnel testing using the stress wave force measurement technique. Aerodynamic modelling was used to estimate the magnitude and distribution of forces and finite element modelling to determine the mechanical response of proposed balance designs. Simulation of balance performance was based on aerodynamic loads and mechanical responses using convolution techniques. Deconvolution was then used to assess balance performance and to guide further design modifications leading to the final balance design. (C) 2001 Elsevier Science Ltd. All rights reserved.

Step size control for efficient discrete element simulation

The step size determines the accuracy of a discrete element simulation. The position and velocity updating calculation uses a pre-calculated table and hence the control of step size can not use the integration formulas for step size control. A step size control scheme for use with the table driven velocity and position calculation uses the difference between the calculation result from one big step and that from two small steps. This variable time step size method chooses the suitable time step size for each particle at each step automatically according to the conditions. Simulation using fixed time step method is compared with that of using variable time step method. The difference in computation time for the same accuracy using a variable step size (compared to the fixed step) depends on the particular problem. For a simple test case the times are roughly similar. However, the variable step size gives the required accuracy on the first run. A fixed step size may require several runs to check the simulation accuracy or a conservative step size that results in longer run times. (C) 2001 Elsevier Science Ltd. All rights reserved.

Differential perturbation of neuronal and glial glutamate transport systems in retinal ischaemia

Glutamate is the major excitatory neurotransmitter in the retina and is removed from the extracellular space by an energy-dependent process involving neuronal and glial cell transporters. The radial glial Muller cells express the glutamate transporter, GLAST, and preferentially accumulate glutamate. However, during an ischaemic episode, extracellular glutamate concentrations may rise to excitotoxic levels. Is this catastrophic rise in extracellular glutamate due to a failure of GLAST? Using immunocytochemistry, we monitored the transport of the glutamate transporter substrate, D-aspartate, in the retina under normal and ischaemic conditions. Two models of compromised retinal perfusion were compared: (1) Anaesthetised rats had their carotid arteries occluded for 7 days to produce a chronic reduction in retinal blood flow. Retinal function was assessed by electroretinography. D-aspartate was injected into the eye for 45 min, Following euthanasia, the retina was processed for D-aspartate. GLAST and glutamate immunocytochemistry. Although reduced retinal perfusion suppresses the electroretinogram b-wave, neither retinal histology, GLAST expression, nor the ability of Muller cells to uptake D-aspartate is affected. As this insult does not appear to cause excitotoxic neuronal damage, these data suggest that GLAST function and glutamate clearance are maintained during periods of reduced retinal perfusion. (2) Occlusion of the central retinal artery for 60 min abolishes retinal perfusion, inducing histological damage and electroretinogram suppression. Although GLAST expression appears to be normal. its ability to transport D-aspartate into Muller cells is greatly reduced. Interestingly, D-aspartate is transported into neuronal cells, i.e. photoreceptors, bipolar and ganglion cells. This suggests that while GLAST is vitally important for the clearance of excess extracellular glutamate, its capability to sustain inward transport is particularly susceptible to an acute ischaemic attack. Manipulation of GLAST function could alleviate the degeneration and blindness that result from ischaemic retinal disease. (C) 2001 Elsevier Science Ltd, All rights reserved.

Rearrangement of diflunisal acyl glucuronide into its beta-glucuronidase-resistant isomers facilitates transport through the small intestine to the colon of the rat

Many non-steroidal anti-inflammatory drugs (NSAIDs) which form acyl glucuronide conjugates as major metabolites have shown an antiproliferative effect on colorectal tumors. This study assesses the extent to which rearrangement of an acyl glucuronide metabolite of a model NSAID into beta -glucuronidase-resistant isomers facilitates its passage through the small intestine to reach the colon. Rats were dosed orally with diflunisal (DF), its acyl glucuronide (DAG) and a mixture of rearrangement isomers (iso-DAG) at 10 mg DF equivalents/kg. The parent drug DF appeared in plasma after all doses, with maximum concentrations of 20.5 +/- 2.5, 28.8 +/- 8.3 and 11.0 +/- 1.6 mug DF/ml respectively, obtained at 3.8 +/- 0.3, 3.6 +/- 1.8 and 7.5 +/- 0.9 hr after the DF, DAG and iso-DAG doses respectively. At 48 hr, 16.2 +/- 3.3, 19.8 +/- 0.8 and 42.9 +/- 10.1% of the doses respectively were recovered in feces, with less than or equal to 1% remaining in the intestine. About half of each dose was recovered as DF and metabolites in 48 hr urine: for DF and DAG doses, the majority was in the first 24 hr urine. whereas for iso-DAG doses, recoveries in the first and second 24 hr periods were similar. The results show that hydrolysis of both DAG and iso-DAG, and absorption of liberated DF, occur during passage through the gut, but that these processes occur more slowly and to a lesser degree for iso-DAG. The intrinsic hydrolytic capacities of various intestinal segments (including contents) towards DAG and iso-DAG were obtained by incubating homogenates under saturating concentrations of DAG/iso-DAG at 37 degreesC. Upper small intestine, lower small intestine, caecum and colon released 2400, 3200, 9200 and 22800 mug DF/hr/g tissue plus contents respectively from DAG substrate, and 18, 10, 140 and 120 mug DF/hr/g tissue plus contents respectively from iso-DAG substrate. The much greater resistance of iso-DAG to hydrolysis appears attributable to its resistance to beta -glucuronidases. The data suggest that in rats dosed with DF, DAG excreted in bile would be substantially hydrolysed in the small intestine and liberated DF reabsorbed, but that portion which rearranges to iso-DAG would likely reach the colon. (C) 2001 Elsevier Science Inc. All rights reserved.

A systematic approach to error isolation in computerized wastewater simulation models

Activated sludge models are used extensively in the study of wastewater treatment processes. While various commercial implementations of these models are available, there are many people who need to code models themselves using the simulation packages available to them, Quality assurance of such models is difficult. While benchmarking problems have been developed and are available, the comparison of simulation data with that of commercial models leads only to the detection, not the isolation of errors. To identify the errors in the code is time-consuming. In this paper, we address the problem by developing a systematic and largely automated approach to the isolation of coding errors. There are three steps: firstly, possible errors are classified according to their place in the model structure and a feature matrix is established for each class of errors. Secondly, an observer is designed to generate residuals, such that each class of errors imposes a subspace, spanned by its feature matrix, on the residuals. Finally. localising the residuals in a subspace isolates coding errors. The algorithm proved capable of rapidly and reliably isolating a variety of single and simultaneous errors in a case study using the ASM 1 activated sludge model. In this paper a newly coded model was verified against a known implementation. The method is also applicable to simultaneous verification of any two independent implementations, hence is useful in commercial model development.

Time-dependent master equation simulation of complex elementary reactions in combustion: Application to the reaction of (CH2)-C-1 with C2H2 from 300-2000 K

Computational simulations of the title reaction are presented, covering a temperature range from 300 to 2000 K. At lower temperatures we find that initial formation of the cyclopropene complex by addition of methylene to acetylene is irreversible, as is the stabilisation process via collisional energy transfer. Product branching between propargyl and the stable isomers is predicted at 300 K as a function of pressure for the first time. At intermediate temperatures (1200 K), complex temporal evolution involving multiple steady states begins to emerge. At high temperatures (2000 K) the timescale for subsequent unimolecular decay of thermalized intermediates begins to impinge on the timescale for reaction of methylene, such that the rate of formation of propargyl product does not admit a simple analysis in terms of a single time-independent rate constant until the methylene supply becomes depleted. Likewise, at the elevated temperatures the thermalized intermediates cannot be regarded as irreversible product channels. Our solution algorithm involves spectral propagation of a symmetrised version of the discretized master equation matrix, and is implemented in a high precision environment which makes hitherto unachievable low-temperature modelling a reality.

Using clusters of computers for large QU-GENE simulation experiments

The QU-GENE Computing Cluster (QCC) is a hardware and software solution to the automation and speedup of large QU-GENE (QUantitative GENEtics) simulation experiments that are designed to examine the properties of genetic models, particularly those that involve factorial combinations of treatment levels. QCC automates the management of the distribution of components of the simulation experiments among the networked single-processor computers to achieve the speedup.

Solution techniques for transport problems involving steep concentration gradients: application to noncatalytic fluid solid reactions

Some efficient solution techniques for solving models of noncatalytic gas-solid and fluid-solid reactions are presented. These models include those with non-constant diffusivities for which the formulation reduces to that of a convection-diffusion problem. A singular perturbation problem results for such models in the presence of a large Thiele modulus, for which the classical numerical methods can present difficulties. For the convection-diffusion like case, the time-dependent partial differential equations are transformed by a semi-discrete Petrov-Galerkin finite element method into a system of ordinary differential equations of the initial-value type that can be readily solved. In the presence of a constant diffusivity, in slab geometry the convection-like terms are absent, and the combination of a fitted mesh finite difference method with a predictor-corrector method is used to solve the problem. Both the methods are found to converge, and general reaction rate forms can be treated. These methods are simple and highly efficient for arbitrary particle geometry and parameters, including a large Thiele modulus. (C) 2001 Elsevier Science Ltd. All rights reserved.

Influence of seaward boundary condition on contaminant transport in unconfined coastal aquifers

Contaminant transport in coastal aquifers is complicated partly due to the conditions at the seaward boundary including seawater intrusion and tidal variations of sea level. Their inclusion in modelling this system will be computationally expensive. Therefore, it will be instructive to investigate the consequence of simplifying the seaward boundary condition by neglecting the seawater density and tidal variations in numerical predictions of contaminant transport in this zone. This paper presents a comparison of numerical predictions for a simplified seaward boundary condition with experimental results for a corresponding realistic one including a saltwater interface and tidal variations. Different densities for contaminants are considered. The comparison suggests that the neglect of the seawater intrusion and tidal variations does not affect noticeably the overall migration rate of the plume before it reaches the saltwater interface. However, numerical prediction shows that a more dense contaminant travels further seaward and part of the solute mass exits under the sea if the seawater density is not included. This is not consistent with the experimental result, which shows that the contaminant travels upwards to the shoreline along the saltwater interface. Neglect of seawater density, therefore, will result in an underestimation of the exit rate of solute mass around the coastline and fictitious migration paths under the seabed. For a less dense contaminant, neglect of seawater density has little effect on numerical prediction of migration paths. (C) 2001 Elsevier Science B.V. All rights reserved.

A simulation model for analysing patient activity in dermatology

We developed a general model to assess patient activity within the primary and secondary health-care sectors following a dermatology outpatient consultation. Based on observed variables from the UK teledermatology trial, the model showed that up to 11 doctor-patient interactions occurred before a patient was ultimately discharged from care. In a cohort of 1000 patients, the average number of health-care visits was 2.4 (range 1-11). Simulation analysis suggested that the most important parameter affecting the total number of doctor-patient Interactions is patient discharge from care following the initial consultation. This implies that resources should be concentrated in this area. The introduction of teledermatology (either realtime or store and forward) changes the values of the model parameters. The model provides a quantitative tool for planning the future provision of dermatology health-care.

Cloning and gastrointestinal expression of rat hephaestin: relationship to other iron transport proteins

The membrane-bound ceruloplasmin homolog hephaestin plays a critical role in intestinal iron absorption. The aims of this study were to clone the rat hephaestin gene and to examine its expression in the gastrointestinal tract in relation to other genes encoding iron transport proteins. The rat hephaestin gene was isolated from intestinal mRNA and was found to encode a protein 96% identical to mouse hephaestin. Analysis by ribonuclease protection assay and Western blotting showed that hephaestin was expressed at high levels throughout the small intestine and colon. Immunofluorescence localized the hephaestin protein to the mature villus enterocytes with little or no expression in the crypts. Variations in iron status had a small but nonsignificant effect on hephaestin expression in the duodenum. The high sequence conservation between rat and mouse hephaestin is consistent with this protein playing a central role in intestinal iron absorption, although its precise function remains to be determined.

Transport of endosomal early antigen I in the rat sciatic nerve and location in cultured neurons

Early endosomal antigen I (EEAI) is known to be a marker of early endosomes and in cultured hippocampal neurons it preferentially localizes to the dendritic but not the axonal compartment. We show in cultured dorsal root ganglia and superior cervical ganglia neurons that EEAI localizes to the cell bodies and the neurites of both sensory and sympathetic neurons. We then show in vivo using a ligated rat sciatic nerve that EEAI significantly accumulates on the proximal side and not on the distal side of the ligation. This suggests that EEAI is transported in the anterograde direction in axons either as part of the homeostatic process or to the nerve ligation site in response to nerve injury. NeuroReport 12:281-284 (C) 2001 Lippincott Williams & Wilkins.