Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

Predictors of local recurrence after treatment of ductal carcinoma in situ - A meta-analysis

lBACKGROUND. Management of patients with ductal carcinoma in situ (DCIS) is a dilemma, as mastectomy provides nearly a 100% cure rate but at the expense of physical and psychologic morbidity. It would be helpful if we could predict which patients with DCIS are at sufficiently high risk of local recurrence after conservative surgery (CS) alone to warrant postoperative radiotherapy (RT) and which patients are at sufficient risk of local recurrence after CS + RT to warrant mastectomy. The authors reviewed the published studies and identified the factors that may be predictive of local recurrence after management by mastectomy, CS alone, or CS + RT. METHODS. The authors examined patient, tumor, and treatment factors as potential predictors for local recurrence and estimated the risks of recurrence based on a review of published studies. They examined the effects of patient factors (age at diagnosis and family history), tumor factors (sub-type of DCIS, grade, tumor size, necrosis, and margins), and treatment (mastectomy, CS alone, and CS + RT). The 95% confidence intervals (CI) of the recurrence rates for each of the studies were calculated for subtype, grade, and necrosis, using the exact binomial; the summary recurrence rate and 95% CI for each treatment category were calculated by quantitative meta-analysis using the fixed and random effects models applied to proportions. RESULTS, Meta-analysis yielded a summary recurrence rate of 22.5% (95% CI = 16.9-28.2) for studies employing CS alone, 8.9% (95% CI = 6.8-11.0) for CS + RT, and 1.4% (95% CI = 0.7-2.1) for studies involving mastectomy alone. These summary figures indicate a clear and statistically significant separation, and therefore outcome, between the recurrence rates of each treatment category, despite the likelihood that the patients who underwent CS alone were likely to have had smaller, possibly low grade lesions with clear margins. The patients with risk factors of presence of necrosis, high grade cytologic features, or comedo subtype were found to derive the greatest improvement in local control with the addition of RT to CS. Local recurrence among patients treated by CS alone is approximately 20%, and one-half of the recurrences are invasive cancers. For most patients, RT reduces the risk of recurrence after CS alone by at least 50%. The differences in local recurrence between CS alone and CS + RT are most apparent for those patients with high grade tumors or DCIS with necrosis, or of the comedo subtype, or DCIS with close or positive surgical margins. CONCLUSIONS, The authors recommend that radiation be added to CS if patients with DCIS who also have the risk factors for local recurrence choose breast conservation over mastectomy. The patients who may be suitable for CS alone outside of a clinical trial may be those who have low grade lesions with little or no necrosis, and with clear surgical margins. Use of the summary statistics when discussing outcomes with patients may help the patient make treatment decisions. Cancer 1999;85:616-28. (C) 1999 American Cancer Society.

Percolative fragmentation of char particles during gasification

Percolative fragmentation was confirmed to occur during gasification of three microporous coal chars. Indirect evidence obtained by the variation of electrical resistivity (ER) with conversion was supported by direct observation of numerous fragments during gasification. The resistivity increases slowly at low conversions and then sharply after a certain conversion value, which is a typical percolation phenomenon suggesting the occurrence of internal fragmentation at high conversion. Two percolation models are applied to interpret the experimental data and determine the percolation threshold. A percolation threshold of 0.02-0.07 was found, corresponding to a critical conversion of 92-96% for fragmentation. The electrical resistivity variation at high conversions is found to be very sensitive to diffusional effects during gasification. Partially burnt samples with a narrow initial particle size range were also observed microscopically, and found to yield a large number of small fragments even when the particles showed no disintegration and chemical control prevailed. It is proposed that this is due to the separation of isolated clusters from the particle surface. The particle size distribution of the fragments was essentially independent of the reaction conditions and the char type, and supported the prediction by percolation theory that the number fraction distribution varies linearly with mass in a log-log plot. The results imply that perimeter fragmentation would occur in practical combustion systems in which the reactions are strongly diffusion affected.

Beach water table fluctuations due to spring-neap tides: moving boundary effects

Tidal water table fluctuations in a coastal aquifer are driven by tides on a moving boundary that varies with the beach slope. One-dimensional models based on the Boussinesq equation are often used to analyse tidal signals in coastal aquifers. The moving boundary condition hinders analytical solutions to even the linearised Boussinesq equation. This paper presents a new perturbation approach to the problem that maintains the simplicity of the linearised one-dimensional Boussinesq model. Our method involves transforming the Boussinesq equation to an ADE (advection-diffusion equation) with an oscillating velocity. The perturbation method is applied to the propagation of spring-neap tides (a bichromatic tidal system with the fundamental frequencies wt and wt) in the aquifer. The results demonstrate analytically, for the first time, that the moving boundary induces interactions between the two primary tidal oscillations, generating a slowly damped water table fluctuation of frequency omega(1) - omega(2), i.e., the spring-neap tidal water table fluctuation. The analytical predictions are found to be consistent with recently published field observations. (C) 2000 Elsevier Science Ltd. All rights reserved.

Modelling solute transport in structured soils: Performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO3-, Cl-, PO43-) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed tone-way ANOVA, p < 0.001) indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale. (C) 2001 Elsevier Science Ltd. All rights reserved.

A context noise model of episodic word recognition

Item noise models of recognition assert that interference at retrieval is generated by the words from the study list. Context noise models of recognition assert that interference at retrieval is generated by the contexts in which the test word has appeared. The authors introduce the bind cue decide model of episodic memory, a Bayesian context noise model, and demonstrate how it can account for data from the item noise and dual-processing approaches to recognition memory. From the item noise perspective, list strength and list length effects, the mirror effect for word frequency and concreteness, and the effects of the similarity of other words in a list are considered. From the dual-processing perspective, process dissociation data on the effects of length, temporal separation of lists, strength, and diagnosticity of context are examined. The authors conclude that the context noise approach to recognition is a viable alternative to existing approaches.

Stochastic structure and individual-tree growth models

The majority of past and current individual-tree growth modelling methodologies have failed to characterise and incorporate structured stochastic components. Rather, they have relied on deterministic predictions or have added an unstructured random component to predictions. In particular, spatial stochastic structure has been neglected, despite being present in most applications of individual-tree growth models. Spatial stochastic structure (also called spatial dependence or spatial autocorrelation) eventuates when spatial influences such as competition and micro-site effects are not fully captured in models. Temporal stochastic structure (also called temporal dependence or temporal autocorrelation) eventuates when a sequence of measurements is taken on an individual-tree over time, and variables explaining temporal variation in these measurements are not included in the model. Nested stochastic structure eventuates when measurements are combined across sampling units and differences among the sampling units are not fully captured in the model. This review examines spatial, temporal, and nested stochastic structure and instances where each has been characterised in the forest biometry and statistical literature. Methodologies for incorporating stochastic structure in growth model estimation and prediction are described. Benefits from incorporation of stochastic structure include valid statistical inference, improved estimation efficiency, and more realistic and theoretically sound predictions. It is proposed in this review that individual-tree modelling methodologies need to characterise and include structured stochasticity. Possibilities for future research are discussed. (C) 2001 Elsevier Science B.V. All rights reserved.

Mathematical models in percutaneous absorption

A number of mathematical models have been used to describe percutaneous absorption kinetics. In general, most of these models have used either diffusion-based or compartmental equations. The object of any mathematical model is to a) be able to represent the processes associated with absorption accurately, b) be able to describe/summarize experimental data with parametric equations or moments, and c) predict kinetics under varying conditions. However, in describing the processes involved, some developed models often suffer from being of too complex a form to be practically useful. In this chapter, we attempt to approach the issue of mathematical modeling in percutaneous absorption from four perspectives. These are to a) describe simple practical models, b) provide an overview of the more complex models, c) summarize some of the more important/useful models used to date, and d) examine sonic practical applications of the models. The range of processes involved in percutaneous absorption and considered in developing the mathematical models in this chapter is shown in Fig. 1. We initially address in vitro skin diffusion models and consider a) constant donor concentration and receptor conditions, b) the corresponding flux, donor, skin, and receptor amount-time profiles for solutions, and c) amount- and flux-time profiles when the donor phase is removed. More complex issues, such as finite-volume donor phase, finite-volume receptor phase, the presence of an efflux. rate constant at the membrane-receptor interphase, and two-layer diffusion, are then considered. We then look at specific models and issues concerned with a) release from topical products, b) use of compartmental models as alternatives to diffusion models, c) concentration-dependent absorption, d) modeling of skin metabolism, e) role of solute-skin-vehicle interactions, f) effects of vehicle loss, a) shunt transport, and h) in vivo diffusion, compartmental, physiological, and deconvolution models. We conclude by examining topics such as a) deep tissue penetration, b) pharmacodynamics, c) iontophoresis, d) sonophoresis, and e) pitfalls in modeling.

ERP 'old/new' effects: memory strength and decisional factor(s)

Event-related potentials (ERPs) were recorded while subjects made old/new recognition judgments on new unstudied words and old words which had been presented at study either once ('weak') or three times ('strong'). The probability of an 'old' response was significantly higher for strong than weak words and significantly higher for weak than new words. Comparisons were made initially between ERPs to new, weak and strong words, and subsequently between ERPs associated with six strength-by-response conditions. The N400 component was found to be modulated by memory trace strength in a graded manner. Its amplitude was most negative in new word ERPs and most positive in strong word ERPs. This 'N400 strength effect' was largest at the left parietal electrode (in ear-referenced ERPs). The amplitude of the late positive complex (LPC) effect was sensitive to decision accuracy (and perhaps confidence). Its amplitude was larger in ERPs evoked by words attracting correct versus incorrect recognition decisions. The LPC effect had a left > right, centro-parietal scalp topography (in ear-referenced ERPs). Hence, whereas, the majority of previous ERP studies of episodic recognition have interpreted results from the perspective of dual-process models, we provide alternative interpretations of N400 and LPC old/new effects in terms of memory strength and decisional factor(s). (C) 2002 Elsevier Science Ltd. All rights reserved.

Effects of alcohol exposure during early life on neuron numbers in the rat hippocampus. I. Hilus neurons and granule cells

We previously showed that 16-day-old rats exposed to a relatively high dose of ethanol at 10-15 postnatal days of age have fewer neurons in the hilus region of the hippocampus compared with controls. Dentate gyrus granule cell numbers, however, showed no statistically significant changes attributable to the ethanol treatment. It is possible that some of the changes in brain morphology, brought about as a result of the exposure to ethanol during early life, may not be manifested until later in life. This question has been further addressed in an extension to our previous study. Wistar rats were exposed to a relatively high daily dose of ethanol on postnatal days 10-15 by placement in a chamber containing ethanol vapour, for 3 h/day. The blood ethanol concentration was found to be similar to430 mg/dl at the end of the period of exposure. Groups of ethanol-treated (ET), separation control (SC), and mother-reared control (MRC) rats were anaesthetised and killed either at 16 or 30 days of age by perfusion with phosphate-buffered 2.5% glutaraldehyde. The Cavalieri principle and the physical disector methods were used to estimate, respectively, the regional volumes and neuron cell numerical densities in the hilus and granule cell regions of the dentate gyrus. The total numbers of neurons in the hilus region and granule cell layer were computed from these estimates. It was found that 16-day-old animals had 398,000-441,000 granule cells, irrespective of group. The numbers of granule cells increased such that by 30 days of age, rats had 487,000-525,500 granule cells. However, there were no significant differences between ethanol-treated rats and their age-matched controls in granule cell numbers. In contrast, ethanol-treated rats had slightly but significantly fewer neurons in the hilus region than did control animals at 16 days of age, but not at 30 days of age. Therefore, it appears that a short period of ethanol exposure during early life can have effects on neuron numbers of some hippocampal neurons, but not others. The effects on hilar neuron numbers, observed as a result of such short periods of ethanol treatment, appeared to be transitory. (C) 2003 Wiley-Liss, Inc.

An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier B.V. All rights reserved.

Effects of age and alcohol exposure during early life on pyramidal cell numbers in the CA1-CA3 region of the rat hippocampus

We have previously shown that exposing rats to a relatively high dose of ethanol during early postnatal life can result in an alteration in spatial learning ability. The hippocampal formation is known to be involved in the control of this ability. The purpose of the present study was to determine whether exposure of rats to ethanol during early postnatal life had either immediate or delayed effects on the numbers of pyramidal cells in the CA1-CA3 subregion of the hippocampus. Wistar rats were exposed to a relatively high daily dose of ethanol at postnatal day 10-15 by placing them for 3 h/day in a chamber containing ethanol vapor. Groups of ethanol-treated (ET), separation control (SC), and mother-reared control (MRC) rats were anesthetized and killed at 16 and 30 days of age by perfusion with phosphate-buffered 2.5% glutaraldehyde. The Cavalieri principle was used to determine the volumes of the CA1 and CA2+CA3 regions. The physical disector method was used to estimate the numerical density of neurons in each of the subdivisions. The total number of pyramidal cells was calculated by multiplying the appropriate estimates of the numerical density by the volume. There were significant age-related reductions in the total numbers of pyramidal cells at 16-30 days of age irrespective of the groups examined. Ethanol treated rats were found to have slightly but significantly fewer pyramidal cell neurons than either the MRC or SC groups. These observations indicate that pyramidal cells in the hippocampus may be vulnerable to a relatively high dose of ethanol exposure during this short period of early postnatal life. (C) 2003 Wiley-Liss, Inc.

Regional Effects of Port Infrastructure: A Spatial CGE Application to Brazil

This article attempts to elucidate one of the mechanisms that link trade barriers, in the form of port costs, and subsequent growth and regional inequality. Prior attention has focused on inland or link costs, but port costs can be considered as a further barrier to enhancing trade liberalization and growth. In contrast to a highway link, congestion at a port may have severe impacts that are spread over space and time whereas highway link congestion may be resolved within several hours. Since a port is part of the transportation network, any congestion/disruption is likely to ripple throughout the hinterland. In this sense, it is important to model properly the role nodal components play in the context of spatial models and international trade. In this article, a spatial computable general equilibrium (CGE) model that is integrated to a transport network system is presented to simulate the impacts of increases in port efficiency in Brazil. The role of ports of entry and ports of exit are explicitly considered to grasp the holistic picture in an integrated interregional system. Measures of efficiency for different port locations are incorporated in the calibration of the model and used as the benchmark in our simulations. Three scenarios are evaluated: (1) an overall increase in port efficiency in Brazil to achieve international standards; (2) efficiency gains associated with decentralization in port management in Brazil; and (3) regionally differentiated increases in port efficiency to reach the boundary of the national efficiency frontier.

Modelling the Distribution of Ischaemic Stroke-Specific Survival Time Using an EM-based Mixture Approach with Random Effects Adjustment

A two-component survival mixture model is proposed to analyse a set of ischaemic stroke-specific mortality data. The survival experience of stroke patients after index stroke may be described by a subpopulation of patients in the acute condition and another subpopulation of patients in the chronic phase. To adjust for the inherent correlation of observations due to random hospital effects, a mixture model of two survival functions with random effects is formulated. Assuming a Weibull hazard in both components, an EM algorithm is developed for the estimation of fixed effect parameters and variance components. A simulation study is conducted to assess the performance of the two-component survival mixture model estimators. Simulation results confirm the applicability of the proposed model in a small sample setting. Copyright (C) 2004 John Wiley Sons, Ltd.

Maternal investment in egg production: environmental and population-specific effects on offspring performance in the snapping shrimp Alpheus nuttingi (Schmitt, 1924) (Decapoda, Alpheidae)

Among marine invertebrates, the overall biomass invested in egg production varies widely within populations, which can result from the interaction of endogenous and exogenous factors. Species that have constant reproduction throughout the year can be good models to study the influence of environmental factors on reproductive processes. We conducted a seasonal comparison of egg production in the intertidal snapping shrimp Alpheus nuttingi, which shows a continuous reproductive pattern, to examine the hypothesis that differences in egg production are driven by environmental conditions and population features. This population showed an uncommon strategy, characterized by females that produce eggs of varying sizes within their clutches, with reduced egg volume when the number of eggs is higher (Spring-Summer). In these seasons, higher temperatures and greater food availability may allow the production of more eggs compared to the Autumn-Winter seasons. Compared to other alpheid shrimps, this population produces small eggs, but in larger numbers. Despite the higher fecundity, the reproductive output is relatively low, this production being supported by the large size of females from the southern Atlantic region. Our findings showed that the egg production of A. nuttingi was greatly influenced by environmental factors. Therefore, this shrimp, and probably other decapods that possess continuous reproduction, adopt different reproductive strategies during the year. (C) Koninklijke Brill NV, Leiden, 2010.