Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. (c) 2004 Elsevier Ltd. All rights reserved.

Self-organization with traveling waves: A case for a convective torus

A traveling wave of BaSO4 in the chlorite-thiourea reaction has shown concentric precipitation patterns upon being triggered by the autocatalyst HOCl. The precipitation patterns show circular rings of alternate null and full precipitation regions. This self-organization appears to be the result of the formation of a convective torus. The formation of the convective torus can be described as a Benard-Marangoni instability with lateral heating.

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.

Homoclinic orbits and chaos in a multimode laser

We show experimentally that under certain conditions the chaotic intensity dynamics of an optically pumped NH3 bidirectional ring laser could be well described in terms of Shil'nikov homoclinic orbits and chaos. We found that the mechanism that resulted in this kind of dynamics of the laser is the competition between effects caused by the mode interaction between the forward and the backward modes of the laser and by the intrinsic single-mode dynamics of the interacting modes. (C) 1997 Optical Society of America.

Elevated striatal and decreased dorsolateral prefrontal cortical activity in response to emotional stimuli in euthymic bipolar disorder: no associations with psychotropic medication load

To examine abnormal patterns of frontal cortical-subcortical activity in response to emotional stimuli in euthymic individuals with bipolar disorder type I in order to identify trait-like, pathophysiologic mechanisms of the disorder. We examined potential confounding effects of total psychotropic medication load and illness variables upon neural abnormalities. We analyzed neural activity in 19 euthymic bipolar and 24 healthy individuals to mild and intense happy, fearful and neutral faces. Relative to healthy individuals, bipolar subjects had significantly increased left striatal activity in response to mild happy faces (p < 0.05, corrected), decreased right dorsolateral prefrontal cortical (DLPFC) activity in response to neutral, mild and intense happy faces, and decreased left DLPFC activity in response to neutral, mild and intense fearful faces (p < 0.05, corrected). Bipolar and healthy individuals did not differ in amygdala activity in response to either emotion. In bipolar individuals, there was no significant association between medication load and abnormal activity in these regions, but a negative relationship between age of illness onset and amygdala activity in response to mild fearful faces (p = 0.007). Relative to those without comorbidities, bipolar individuals with comorbidities showed a trend increase in left striatal activity in response to mild happy faces. Abnormally increased striatal activity in response to potentially rewarding stimuli and decreased DLPFC activity in response to other emotionally salient stimuli may underlie mood instabilities in euthymic bipolar individuals, and are more apparent in those with comorbid diagnoses. No relationship between medication load and abnormal neural activity in bipolar individuals suggests that our findings may reflect pathophysiologic mechanisms of the illness rather than medication confounds. Future studies should examine whether this pattern of abnormal neural activity could distinguish bipolar from unipolar depression.

Control of chaotic motion in a dual-spin spacecraft with nutational damping

Control of chaotic vibrations in a dual-spin spacecraft with an axial nutational damper is achieved using two techniques. The control methods are implemented on two realistic spacecraft parameter configurations that have been found to exhibit chaotic instability when a sinusoidally varying torque is applied to the spacecraft for a range of forcing amplitudes and frequencies. Such a torque, in practice, may arise under malfunction of the control system or from an unbalanced rotor. Chaotic instabilities arising from these torques could introduce uncertainties and irregularities into a spacecraft's attitude motion and, consequently, could have disastrous effects on its operation. The two control methods, recursive proportional feedback and continuous delayed feedback, are recently developed techniques for control of chaotic motion in dynamic systems. Each technique is outlined and the effectiveness on this model compared and contrasted. Numerical simulations are performed, and the results are studied by means of time history, phase space, Poincare map, Lyapunov characteristic exponents, and bifurcation diagrams.

Control of chaotic instability in a dual-spin spacecraft with dissipation using energy methods

Control of chaotic instability in a rotating multibody system in the form of a dual-spin spacecraft with an axial nutational damper is achieved using an algorithm derived using energy methods. The control method is implemented on two realistic spacecraft parameter configurations which have been found to exhibit chaotic instability when a sinusoidally varying torque is applied to the spacecraft for a range of forcing amplitudes and frequencies. Such a torque, in practice, may arise under malfunction of the control system or from an unbalanced rotor. Chaotic instabilities arising from these torques could introduce uncertainties and irregularities into a spacecraft's attitude and consequently impair pointing accuracy. The control method is formulated from nutational stability results derived using an energy sink approximation for a dual-spin spacecraft with an asymmetric platform and axisymmetric rotor. The effectiveness of the control method is shown numerically and the results are studied by means of time history, phase space, Poincare map, Lyapunov characteristic exponents and Bifurcation diagrams.

Review and future directions in the modelling and control of continuous drum granulation

Many granulation plants operate well below design capacity, suffering from high recycle rates and even periodic instabilities. This behaviour cannot be fully predicted using the present models. The main objective of the paper is to provide an overview of the current status of model development for granulation processes and suggest future directions for research and development. The end-use of the models is focused on the optimal design and control of granulation plants using the improved predictions of process dynamics. The development of novel models involving mechanistically based structural switching methods is proposed in the paper. A number of guidelines are proposed for the selection of control relevant model structures. (C) 2002 Published by Elsevier Science B.V.

The prediction of sharkskin instability observed during film blowing

The characteristics of sharkskin surface instability for linear low density polyethylene are studied as a function of film blowing processing conditions. By means of scanning electron microscopy and surface profilometry, is it found that for the standard industrial die geometry studied, sharkskin only occurs on the inside of the film bubble. Previous work suggests that this instability may be due to critical extensional stress levels at the exit of the die. Isothermal integral viscoelastic simulations of the annular extrusion process are reported, and confirm that the extensional stress at the die exit is large enough to cause local melt rupture. However the extensional stress level at the outer die wall predicts melt rupture of the outside bubble surface also, which contradicts the experimental findings. A significant temperature gradient is expected to exist across the die gap at the exit of the die, due to the external heating of the die and the low conductivity, of the polymer melt. It is shown that a gradient of 20 degreesC is required to cause sharkskin to only appear on the inner bubble surface.

Experimental investigation of contaminant transport in coastal groundwater

Contaminant transport in coastal aquifers is of increasing interest since, with the development of coastal areas, contaminants from surface sources may enter coastal aquifers and pollute the groundwater flow. Coastal groundwater flow is complicated because of the presence of a freshwater-saltwater diffusion zone and the tidal variation of sea level at the seaward end. This paper investigates experimentally the behaviour of contaminant plumes with different densities in an unconfined coastal aquifer. Experiments were performed in a flow tank filled with glass beads as the porous medium. Results show that the dense contaminant has a more diffusive front than the less dense one in the seaward direction towards the coastline. The plume becomes more diffusive when it travels closer to the saltwater interface. On the contrary, the less dense contaminant presents a relatively sharp outline. It tends to migrate in the upper portion of the aquifer and exits in a concentrated manner over a small discharge area at the coastline, not further seaward under the sea. Non-dimensional parameters show that instabilities occur in our experiments for a density difference of 1.2% or larger between the contaminant and the ambient water. The experimental results provide guidance for field monitoring and numerical modelling. (C) 2002 Elsevier Science Ltd. All rights reserved.

A unified model of flames as gasdynamic discontinuities

Viewed on a hydrodynamic scale, flames in experiments are often thin so that they may be described as gasdynamic discontinuities separating the dense cold fresh mixture from the light hot burned products. The original model of a flame as a gasdynamic discontinuity was due to Darrieus and to Landau. In addition to the fluid dynamical equations, the model consists of a flame speed relation describing the evolution of the discontinuity surface, and jump conditions across the surface which relate the fluid variables on the two sides of the surface. The Darrieus-Landau model predicts, in contrast to observations, that a uniformly propagating planar flame is absolutely unstable and that the strength of the instability grows with increasing perturbation wavenumber so that there is no high-wavenumber cutoff of the instability. The model was modified by Markstein to exhibit a high-wavenumber cutoff if a phenomenological constant in the model has an appropriate sign. Both models are postulated, rather than derived from first principles, and both ignore the flame structure, which depends on chemical kinetics and transport processes within the flame. At present, there are two models which have been derived, rather than postulated, and which are valid in two non-overlapping regions of parameter space. Sivashinsky derived a generalization of the Darrieus-Landau model which is valid for Lewis numbers (ratio of thermal diffusivity to mass diffusivity of the deficient reaction component) bounded away from unity. Matalon & Matkowsky derived a model valid for Lewis numbers close to unity. Each model has its own advantages and disadvantages. Under appropriate conditions the Matalon-Matkowsky model exhibits a high-wavenumber cutoff of the Darrieus-Landau instability. However, since the Lewis numbers considered lie too close to unity, the Matalon-Matkowsky model does not capture the pulsating instability. The Sivashinsky model does capture the pulsating instability, but does not exhibit its high-wavenumber cutoff. In this paper, we derive a model consisting of a new flame speed relation and new jump conditions, which is valid for arbitrary Lewis numbers. It captures the pulsating instability and exhibits the high-wavenumber cutoff of all instabilities. The flame speed relation includes the effect of short wavelengths, not previously considered, which leads to stabilizing transverse surface diffusion terms.

Stability of planar flames as gasdynamic discontinuities

The stability of a steadily propagating planar premixed flame has been the subject of numerous studies since Darrieus and Landau showed that in their model flames are unstable to perturbations of any wavelength. Moreover, the instability was shown to persist even for very small wavelengths, i.e. there was no high-wavenumber cutoff of the instability. In addition to the Darrieus-Landau instability, which results from thermal expansion, analysis of the diffusional thermal model indicates that premixed flames may exhibit cellular and pulsating instabilities as a consequence of preferential diffusion. However, no previous theory captured all the instabilities including a high-wavenumber cutoff for each. In Class, Matkowsky & Klimenko (2003) a unified theory is proposed which, in appropriate limits and under appropriate assumptions, recovers all the relevant previous theories. It also includes additional new terms, not present in previous theories. In the present paper we consider the stability of a uniformly propagating planar flame as a solution of the unified model. The results are then compared to those based on the models of Darrieus-Landau, Sivashinsky and Matalon-Matkowsky. In particular, it is shown that the unified model is the only model to capture the Darrieus-Landau, cellular and pulsating instabilities including a high-wavenumber cutoff for each.

Palaeomagnetic study of a subaerial volcanic ridge (Sao Jorge Island, Azores) for the past 1.3 Myr: evidence for the Cobb Mountain Subchron, volcano flank instability and tectonomagmatic implications

We present a palaeomagnetic study on 38 lava flows and 20 dykes encompassing the past 1.3 Myr on S. Jorge Island (Azores ArchipelagoNorth Atlantic Ocean). The sections sampled in the southeastern and central/western parts of the island record reversed and normal polarities, respectively. They indicate a mean palaeomagnetic pole (81.3 degrees N, 160.7 degrees E, K= 33 and A95= 3.4 degrees) with a latitude shallower than that expected from Geocentric Axial Dipole assumption, suggesting an effect of non-dipolar components of the Earth magnetic field. Virtual Geomagnetic Poles of eight flows and two dykes closely follow the contemporaneous records of the Cobb Mountain Subchron (ODP/DSDP programs) and constrain the age transition from reversed to normal polarity at ca. 1.207 +/- 0.017 Ma. Volcano flank instabilities, probably related to dyke emplacement along an NNWSSE direction, led to southwestward tilting of the lava pile towards the sea. Two spatially and temporally distinct dyke systems have been recognized on the island. The eastern is dominated by NNWSSE trending dykes emplaced before the end of the Matuyama Chron, whereas in the central/western parts the eruptive fissures oriented WNWESE controlled the westward growth of the S. Jorge Island during the Brunhes Chron. Both directions are consistent with the present-day regional stress conditions deduced from plate kinematics and tectonomorphology and suggest the emplacement of dykes along pre-existing fractures. The distinct timing and location of each dyke system likely results from a slight shift of the magmatic source.

Estudo geológico e geotécnico de taludes rochosos em S. Pedro da Afurada (Vila Nova de Gaia): contributos para a avaliação da estabilidade estrutural

Mestrado em Engenharia Geotécnica e Geoambiente

Geotecnia urbana da zona ribeirinha de Vila Nova de Gaia: criação de uma base de dados SIG geotécnica

O presente trabalho apresenta os resultados dos estudos geotécnicos e de uma base de dados da zona ribeirinha de Vila Nova de Gaia, com o objectivo de compreender melhor os aspectos geotécnicos em ambiente urbano numa área sensível com um registo histórico de instabilidade de taludes rochosos. Além disso, os escassos estudos científicos recentes de natureza geológica e geotécnica em Vila Nova de Gaia justificam o estudo exploratório da geotecnia urbana da zona ribeirinha de Vila Nova de Gaia. A importância de Vila Nova de Gaia como a terceira maior cidade portuguesa e como centro de intensa actividade económica e cultural despoleta uma constante necessidade de expansão. O aumento da densidade populacional acarreta a realização de projectos complexos de engenharia, utilizando o subsolo para a construção e, com frequência, em terrenos com características geotécnicas desfavoráveis. As cidades de Vila Nova de Gaia e do Porto foram sendo edificadas ao longo de encostas numa plataforma litoral caracterizada por uma vasta área aplanada, inclinando ligeiramente para Oeste. Esta plataforma foi cortada pelo Rio Douro num vale encaixado de vertentes abruptas, nas quais se localizam as zonas ribeirinhas das duas cidades. Este trabalho envolveu, inicialmente, uma caracterização topográfica, morfoestrutural, geotectónica e geomecânica da área de estudo e, numa fase posterior, o desenvolvimento duma base de dados geotécnica. Todos os dados geológicos e geotécnicos locais e os estudos geotécnicos levados a cabo in situ pelas diversas empresas e instituições foram representados cartograficamente numa base apoiada pelos Sistemas de Informação Geográfica (SIG). Esta metodologia inter‐disciplinar foi de grande valor para um melhor conhecimento dos riscos geológico‐geotécnicos ao longo das margens do Rio Douro. De facto, a cartografia geotécnica da zona ribeirinha de Vila Nova de Gaia deve constituir uma ferramenta importante para uma previsão mais rigorosa de futuras instabilidades de taludes e um bom instrumento para a gestão do espaço urbano.