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.

A note on the decomposition of the health concentration index

In recent work, the concentration index has been widely used as a measure of income-related health inequality. The purpose of this note is to illustrate two different methods for decomposing the overall health concentration index using data collected from a Short Form (SF-36) survey of the general Australian population conducted in 1995. For simplicity, we focus on the physical functioning scale of the SF-36. Firstly we examine decomposition 'by component' by separating the concentration index for the physical functioning scale into the ten items on which it is based. The results show that the items contribute differently to the overall inequality measure, i.e. two of the items contributed 13% and 5%, respectively, to the overall measure. Second, to illustrate the 'by subgroup' method we decompose the concentration index by employment status. This involves separating the population into two groups: individuals currently in employment; and individuals not currently employed. We find that the inequality between these groups is about five times greater than the inequality within each group. These methods provide insights into the nature of inequality that can be used to inform policy design to reduce income related health inequalities. Copyright (C) 2002 John Wiley Sons, Ltd.

Formal methods within a totally functional approach to programming

Taking functional programming to its extremities in search of simplicity still requires integration with other development (e.g. formal) methods. Induction is the key to deriving and verifying functional programs, but can be simplified through packaging proofs with functions, particularly folds, on data (structures). Totally Functional Programming avoids the complexities of interpretation by directly representing data (structures) as platonic combinators - the functions characteristic to the data. The link between the two simplifications is that platonic combinators are a kind of partially-applied fold, which means that platonic combinators inherit fold-theoretic properties, but with some apparent simplifications due to the platonic combinator representation. However, despite observable behaviour within functional programming that suggests that TFP is widely-applicable, significant work remains before TFP as such could be widely adopted.

Choice of macroinvertebrate metrics to evaluate stream conditions in Atlantic Forest, Brazil

The development of biomonitoring programs based on the macroinvertebrate community requires the understanding of species distribution patterns, as well as of the responses of the community to anthropogenic stressors. In this study, 49 metrics were tested as potential means of assessing the condition of 29 first- and second-order streams located in areas of differing types of land use in So Paulo State, Brazil. Of the sampled streams, 15 were in well-preserved regions in the Atlantic Forest, 5 were among sugarcane cultivations, 5 were in areas of pasture, and 4 were among eucalyptus plantations. The metrics were assessed against the following criteria: (1) predictable response to the impact of human activity; (2) highest taxonomic resolution, and (3) operational and theoretical simplicity. We found that 18 metrics were correlated with the environmental and spatial predictors used, and seven of these satisfied the selection criteria and are thus candidates for inclusion in a multimetric system to assess low-order streams in So Paulo State. These metrics are family richness; Ephemeroptera, Plecoptera and Trichoptera (EPT) richness; proportion of Megaloptera and Hirudinea; proportion of EPT; Shannon diversity index for genus; and adapted Biological Monitoring Work Party biotic index.

Transanal endorectal pull-through in children with Hirschsprung`s disease-technical refinements and comparison of results with the Duhamel procedure

Purpose: Transanal endorectal pull-through (TEPT) has drastically changed the treatment of Hirschsprung`s disease (HD). A short follow-up of children Submitted to TEPT reveals results that are similar to the classic transabdominal pull-through procedures. However, few reports compare the late results of TEPT with transabdominal pull-through procedures with respect to complication rates and the fecal continence. The aims of the present work are to describe some technical refinements that we introduced in the procedure and to compare the short and long-term outcome of TEPT with the outcomes of a group of patients with HD who previously underwent the Duhamel procedure. Methods: Thirty-five patients who underwent TEPT were prospectively studied and compared to a group of 29 patients who were treated with colostomy followed by a classical Duhamel pull-through. The main modifications introduced in the TEPT group were no preoperative colon preparation, operation conducted under general anesthesia in addition to regional sacral anesthesia, use of only one purse-string Suture in the rectal mucosa before transanal submucosal dissection, and no use of retractors and electrocautery during file submucosal dissection. Results: The most frequent early complications of TEPT group were perineal dermatitis (22.8%) and anastomotic strictures (8.6%). The comparison with patients who underwent Duhamel procedure revealed no difference in the incidence of preoperative enterocolitis, the patients of the TEPT group were younger at the time of diagnosis and of surgery, they had shorter operating times, and they began oral feeding more quickly after the operation. The incidence of wound infection was lower in the TEPT group. Moreover, the TEPT and Duhamel groups showed no difference in the incidences of mortality, postoperative partial continence, and total incontinence. Although the incidences of complete continence and postoperative enterocolitis were not different, a tendency to the increased incidence in the TEPT group was observed. Conclusions: This study further supports the technical advantages, the simplicity, and the decreased incidence of complications of a primary TEPT procedure when compared to a classical form of pull-through. Sonic technical refinements are described, and no preoperative colon preparation was necessary for the patients studied here. The results show that the long-term outcomes of the modified TEPT procedure are generally better than those obtained with classical approaches. Published by Elsevier Inc.

Discrete and continuous models for dry masonry columns

The dynamic response of dry masonry columns can be approximated with finite-difference equations. Continuum models follow by replacing the difference quotients of the discrete model by corresponding differential expressions. The mathematically simplest of these models is a one-dimensional Cosserat theory. Within the presented homogenization context, the Cosserat theory is obtained by making ad hoc assumptions regarding the relative importance of certain terms in the differential expansions. The quality of approximation of the various theories is tested by comparison of the dispersion relations for bending waves with the dispersion relation of the discrete theory. All theories coincide with differences of less than 1% for wave-length-block-height (L/h) ratios bigger than 2 pi. The theory based on systematic differential approximation remains accurate up to L/h = 3 and then diverges rapidly. The Cosserat model becomes increasingly inaccurate for L/h < 2 pi. However, in contrast to the systematic approximation, the wave speed remains finite. In conclusion, considering its relative simplicity, the Cosserat model appears to be the natural starting point for the development of continuum models for blocky structures.

A simple solution of the laminar dam break wave

The sudden release of a mass of fluid in a channel generates a highly unsteady flow motion, called dam break wave. While industrial fluids exhibit sometimes non-Newtonian behaviours, the viscous fluid flow assumption remains a useful approximation for simplified analyses. In this study, new solutions of laminar dam break wave are proposed for a semi-infinite reservoir based upon the method of characteristics. The solutions yield simple explicit expressions of the wave front location, wave front celerity and instantaneous free-surface profiles that compare favourably with experimental observations. Both horizontal and sloping channel configurations are treated. The simplicity of the equations may allow future extension to more complicated fluid flows.

Development and implementation of the population Fisher information matrix for the evaluation of population pharmacokinetic designs.

In population pharmacokinetic studies, the precision of parameter estimates is dependent on the population design. Methods based on the Fisher information matrix have been developed and extended to population studies to evaluate and optimize designs. In this paper we propose simple programming tools to evaluate population pharmacokinetic designs. This involved the development of an expression for the Fisher information matrix for nonlinear mixed-effects models, including estimation of the variance of the residual error. We implemented this expression as a generic function for two software applications: S-PLUS and MATLAB. The evaluation of population designs based on two pharmacokinetic examples from the literature is shown to illustrate the efficiency and the simplicity of this theoretic approach. Although no optimization method of the design is provided, these functions can be used to select and compare population designs among a large set of possible designs, avoiding a lot of simulations.

On the overabundance of light rare earth elements in terrestrial zircons and its implication for Earth's earliest magmatic differentiation

We present whole-rock and zircon rare earth element (REE) data from two early Archaean gneisses (3.81 Ga and 3.64 Ga) from the Itsaq gneiss complex, south-west Greenland. Both gneisses represent extremely rare examples of unaltered, fresh and relatively undeformed igneous rocks of such antiquity. Cathodoluminescence imaging of their zircons indicates a single crystallisation episode with no evidence for either later metamorphic and/or anatectic reworking or inheritance of earlier grains. Uniform, single-population U/Pb age data confirm the structural simplicity of these zircons. One sample, a 3.64 Ga granodioritic gneiss from the Gothabsfjord, yields a chondrite-normalised REE pattern with a positive slope from La to Lu as well as substantial positive Ce and slight negative Eu anomalies, features generally considered to be typical of igneous zircon. In contrast, the second sample, a 3.81 Ga tonalite from south of the Isua Greenstone Belt, has variable but generally much higher light REE abundances, with similar middle to heavy REE. Calculation of zircon/melt distribution coefficients (D-REE(zircon/melt)) from each sample yields markedly different values for the trivalent REE (i.e. Ce and Eu omitted) and simple application of one set of D-REE(zircon/melt) to model the melt composition for the other sample yields concentrations that are in error by up to two orders of magnitude for the light REE (La-Nd). The observed light REE overabundance in the 3.81 Ga tonalite is a commonly observed feature in terrestrial zircons for which a number of explanations ranging from lattice strain to disequilibrium crystallisation have been proposed and are further investigated herein. Regardless of the cause of light REE overabundance, our study shows that simple application of zircon/melt distribution coefficients is not an unambiguous method for ascertaining original melt composition. In this context, recent studies that use REE data to claim that > 4.3 Ga Hadean detrital zircons originally crystallised from an evolved magma, in turn suggesting the operation of geological processes in the early Earth analogous to those of the present day (e.g. subduction and melting of hydrated oceanic crust), must be regarded with caution. Indeed, comparison of terrestrial Hadean and > 3.9 Ga lunar highland zircons shows remarkable similarities in the light REE, even though subduction processes that have been used to explain the terrestrial zircons have never operated on the Moon. (C) 2002 Elsevier Science B.V. All rights reserved.

Vacancy solution theory for binary adsorption equilibria in heterogeneous carbon

A heterogeneous modified vacancy solution model of adsorption developed is evaluated. The new model considers the adsorption process through a mass-action law and is thermodynamically consistent, while maintaining the simplicity in calculation of multicomponent adsorption equilibria, as in the original vacancy solution theory. It incorporates the adsorbent heterogeneity through a pore-width-related potential energy, represented by Steele's 10-4-3 potential expression. The experimental data of various hydrocarbons, CO2 and SO2 on four different activated carbons - Ajax, Norit, Nuxit, and BPL - at multiple temperatures over a wide range of pressures were studied by the heterogeneous modified VST model to obtain the isotherm parameters and micropore-size distribution of carbons. The model successfully correlates the single-component adsorption equilibrium data for all compounds studied on various carbons. The fitting results for the vacancy occupancy parameter are consistent with the pressure change on different carbons, and the effect of pore heterogeneity is important in adsorption at elevated pressure. It predicts binary adsorption equilibria better than the IAST scheme, reflecting the significance of molecular size nonideality.

The Galilean turn in population ecology

The standard mathematical models in population ecology assume that a population's growth rate is a function of its environment. In this paper we investigate an alternative proposal according to which the rate of change of the growth rate is a function of the environment and of environmental change. We focus on the philosophical issues involved in such a fundamental shift in theoretical assumptions, as well as on the explanations the two theories offer for some of the key data such as cyclic populations. We also discuss the relationship between this move in population ecology and a similar move from first-order to second-order differential equations championed by Galileo and Newton in celestial mechanics.

A class B switch-mode assisted linear amplifier

A switch-mode assisted linear amplifier (SMALA) combining a linear (Class B) and a switch-mode (Class D) amplifier is presented. The usual single hysteretic controlled half-bridge current dumping stage is replaced by two parallel buck converter stages, in a parallel voltage controlled topology. These operate independently: one buck converter sources current to assist the upper Class B output device, and a complementary converter sinks current to assist the lower device. This topology lends itself to a novel control approach of a dead-band at low power levels where neither class D amplifier assists, allowing the class B amplifier to supply the load without interference, ensuring high fidelity. A 20 W implementation demonstrates 85% efficiency, with distortion below 0.08% measured across the full audio bandwidth at 15 W. The class D amplifier begins assisting at 2 W, and below this value, the distortion was below 0.03%. Complete circuitry is given, showing the simplicity of the additional class D amplifier and its corresponding control circuitry.

Modelagem da dispersão de odores em ambientes urbanos

O presente trabalho investigou o problema da modelagem da dispersão de compostos odorantes em presença de obstáculos (cúbicos e com forma complexa) sob condição de estabilidade atmosférica neutra. Foi empregada modelagem numérica baseada nas equações de transporte (CFD1) bem como em modelos algébricos baseados na pluma Gausseana (AERMOD2, CALPUFF3 e FPM4). Para a validação dos resultados dos modelos e a avaliação do seu desempenho foram empregados dados de experimentos em túnel de vento e em campo. A fim de incluir os efeitos da turbulência atmosférica na dispersão, dois diferentes modelos de sub-malha associados à Simulação das Grandes Escalas (LES5) foram investigados (Smagorinsky dinâmico e WALE6) e, para a inclusão dos efeitos de obstáculos na dispersão nos modelos Gausseanos, foi empregado o modelo PRIME7. O uso do PRIME também foi proposto para o FPM como uma inovação. De forma geral, os resultados indicam que o uso de CFD/LES é uma ferramenta útil para a investigação da dispersão e o impacto de compostos odorantes em presença de obstáculos e também para desenvolvimento dos modelos Gausseanos. Os resultados também indicam que o modelo FPM proposto, com a inclusão dos efeitos do obstáculo baseado no PRIME também é uma ferramenta muito útil em modelagem da dispersão de odores devido à sua simplicidade e fácil configuração quando comparado a modelos mais complexos como CFD e mesmo os modelos regulatórios AERMOD e CALPUFF. A grande vantagem do FPM é a possibilidade de estimar-se o fator de intermitência e a relação pico-média (P/M), parâmetros úteis para a avaliação do impacto de odores. Os resultados obtidos no presente trabalho indicam que a determinação dos parâmetros de dispersão para os segmentos de pluma, bem como os parâmetros de tempo longo nas proximidades da fonte e do obstáculo no modelo FPM pode ser melhorada e simulações CFD podem ser usadas como uma ferramenta de desenvolvimento para este propósito. Palavras chave: controle de odor, dispersão, fluidodinâmica computacional, modelagem matemática, modelagem gaussiana de pluma flutuante, simulação de grandes vórtices (LES).

Estudo da resistência não drenada de solo de baixa consistência por meio de ensaios de campo e laboratório

A variabilidade natural dos solos torna complexo o conhecimento de suas propriedades na elaboração de projetos geotécnicos, sendo a determinação da resistência ao cisalhamento não drenada um parâmetro importante nas análises de estabilidade de solos moles. Os ensaios de laboratório de cone e palheta, não convencionais, os ensaios de campo de palheta e piezocone e os ensaios de compressão simples e triaxial não adensado e não drenado foram utilizados para mensurar a resistência não drenada de uma camada de argila marinha mole localizada na planície costeira central brasileira. Os ensaios de laboratório foram realizados em amostras indeformadas coletadas com amostradores de pistão estacionário em vertical próxima à realização dos ensaios de campo. O sítio foi investigado preliminarmente por sondagens de simples reconhecimento, sendo apresentado o perfil estratigráfico por meio de modelagem computacional. Foram também realizados ensaios para caracterização física (análise granulométrica, teor de umidade, limites de liquidez e plasticidade, densidade real dos grãos) e mineralógica (difração de raios X), e ensaios de adensamento para obtenção do histórico de tensões e classificação de qualidade das amostras indeformadas. Os valores de resistência não drenada obtidos pelos ensaios de laboratório foram comparados ao perfil contínuo de resistência determinado empiricamente pelo ensaio de piezocone, com fator de cone Nkt calibrado pelo ensaio de palheta de campo, apresentando boa concordância, com a variabilidade natural do solo influenciando de forma preponderante a qualidade das amostras na variação entre os resultados. Os valores de resistência obtidos pelos ensaios de laboratório de cone e palheta foram comparados entre si, apresentando boa compatibilidade. Ambos, quando comparados ao ensaio de palheta de campo, não apresentaram boa concordância. Os resultados de resistência obtidos pelos ensaios de compressão simples e triaxial apresentaram boa compatibilidade com os resultados do ensaio de laboratório de cone, o que não ocorreu com os resultados do ensaio de laboratório de palheta. Na comparação entre a resistência normalizada pela tensão de sobreadensamento obtida pelos diversos métodos e algumas correlações empíricas da literatura internacional, foi observado para as amostras de solo com índice de plasticidade superior a 60% boa concordância com as correlações de Mesri (1975) e Jamiolkowski et al (1985). Os ensaios não convencionais apresentaram boa confiabilidade, que aliado a simplicidade e agilidade de execução, justificam a difusão destes na prática da investigação geotécnica brasileira como método alternativo para complementar e dar suporte às estimativas de resistência não drenada de solos moles.

Drying kinetics of jatropha seeds

Given the necessity of developing jatropha cultivation equipment, this work adjusted different mathematical models to experimental data obtained from the drying of jatropha seeds submitted to different drying conditions and selected the best model to describe the drying process. The experiment was carried out at the Federal Institute of Goiás - Rio Verde Campus. Seeds with initial moisture content of approximately 0.50 (kg water/kg dry matter) were dried in a forced air-ventilated oven, at temperatures of 45, 60, 75, 90 and 105°C to moisture content of 0.10 ± 0.005 (kg water/kg dry matter). The experimental data were adjusted to 11 mathematical models to represent the drying process of agricultural products. The models were compared using the coefficient of determination, chi-square test, relative mean error, estimated mean error and residual distribution. It was found that the increase in the air temperature caused a reduction in the drying time of seeds. The models Midilli and Two Terms were suitable to represent the drying process of Jatropha seeds and between them the use of the Midili model is recommended due to its greater simplicity.