CXTANNEAL: an improved program for estimating solute transport parameters

CXTANNEAL is a program for analysing contaminant transport in soils. The code, written in Fortran 77, is a modified version of CXTFIT, a commonly used package for estimating solute transport parameters in soils. The improvement of the present code is that it includes simulated annealing as the optimization technique for curve fitting. Tests with hypothetical data show that CXTANNEAL performs better than the original code in searching for optimal parameter estimates. To reduce the computational time, a parallel version of CXTANNEAL (CXTANNEAL_P) was also developed. (C) 1999 Elsevier Science Ltd. All rights reserved.

Local and global error models to improve uncertainty quantification

In groundwater applications, Monte Carlo methods are employed to model the uncertainty on geological parameters. However, their brute-force application becomes computationally prohibitive for highly detailed geological descriptions, complex physical processes, and a large number of realizations. The Distance Kernel Method (DKM) overcomes this issue by clustering the realizations in a multidimensional space based on the flow responses obtained by means of an approximate (computationally cheaper) model; then, the uncertainty is estimated from the exact responses that are computed only for one representative realization per cluster (the medoid). Usually, DKM is employed to decrease the size of the sample of realizations that are considered to estimate the uncertainty. We propose to use the information from the approximate responses for uncertainty quantification. The subset of exact solutions provided by DKM is then employed to construct an error model and correct the potential bias of the approximate model. Two error models are devised that both employ the difference between approximate and exact medoid solutions, but differ in the way medoid errors are interpolated to correct the whole set of realizations. The Local Error Model rests upon the clustering defined by DKM and can be seen as a natural way to account for intra-cluster variability; the Global Error Model employs a linear interpolation of all medoid errors regardless of the cluster to which the single realization belongs. These error models are evaluated for an idealized pollution problem in which the uncertainty of the breakthrough curve needs to be estimated. For this numerical test case, we demonstrate that the error models improve the uncertainty quantification provided by the DKM algorithm and are effective in correcting the bias of the estimate computed solely from the MsFV results. The framework presented here is not specific to the methods considered and can be applied to other combinations of approximate models and techniques to select a subset of realizations

Sorbent characterization for circulating fluidized bed scrubbers

Harmful sulfur dioxide (SO2) emissions from power plants have increasingly been restricted since the 1970’s. Circulating fluidized bed (CFB) scrubber is a dry flue gas desulfurization method of absorbing SO2 out of the flue gas with sorbent. In current commercial plants, the used sorbent is commercial or on-site hydrated calcium hydroxide. The CFB scrubber process is characterized by a close but adequate approach to the flue gas saturation temperature that is achieved by spraying water to the absorber followed by a particulate control device. Very high SO2 removal is achieved along with a dry byproduct that is continuously recirculated back to the absorber for enhanced sorbent utilization. The aim of this work is to develop a method that would characterize the reactivity of sorbents used in CFB scrubbers and to conclude how different process parameters and sorbent properties affect the sulfur absorption. The developed characterization method is based on a fixed bed of sorbent and inert silica sand, through which an SO2 containing gas mixture is led. The reaction occurs in the bed and the SO2 concentration in the outlet as a function of time, a breakthrough curve, is obtained from the analyzer. Reactivity of the sorbents are evaluated by the absorbed sulfur amount. Results suggest that out of process parameters, lower SO2 concentration, lower temperature and higher moisture content enhance the desulfurization. Between different sorbents, specific surface area seems to be the most significant parameter. Large surface area linearly leads to more efficient desulfurization. Overall, the solid conversion levels in the tests were very low creating uncertainty to the validity of the results. New desing is being planned to overcome the problems of the device.

Avaliação da capacidade de adsorção do CO2 em zeólita 12X com gases sintéticos e originados da pirólise do lodo de esgoto

This work depicts a study of the adsorption of carbon dioxide on zeolite 13X. The activities were divided into four stages: study batch adsorption capacity of the adsorbent with synthetic CO2 (4%), fixed bed dynamic evaluation with the commercial mixture of gases (4% CO2, 1.11% CO, 1 2% H2, 0.233% CH4, 0.1% C3, 0.0233% C4 argon as inert closing balance), fixed bed dynamic modeling and evaluation of the breakthrough curve of CO2 originated from the pyrolysis of sewage sludge. The sewage sludge and the adsorbent were characterized by analysis TG / DTA, SEM, XRF and BET. Adsorption studies were carried out under the following operating conditions: temperature 40 °C (for the pyrolysis of the sludge T = 600 °C), pressures of 0.55 to 5.05 bar (batch process), flow rate of the gaseous mixture between 50 - 72 ml/min and the adsorbent masses of 10, 15 and 20 g (fixed bed process). The time for the adsorption batch was 7 h and on the fixed bed was around 180 min. The results of this study showed that in batch adsorption process step with zeolite 13X is efficient and the mass of adsorbed CO2 increases with the increases pressure, decreases with temperature increases and rises due the increase of activation temperature adsorbent. In the batch process were evaluated the breakthrough curves, which were compared with adsorption isotherms represented by the models of Langmuir, Freündlich and Toth. All models well adjusted to the experimental points, but the Langmuir model was chosen in view of its use in the dynamic model does not have implications for adsorption (indeterminacy and larger number of parameters such as occurred with others) in solving the equation. In the fixed bed dynamic study with the synthetic gas mixture, 20 g of mass adsorbent showed the maximum adsorption percentage 46.7% at 40 °C temperature and 50 mL/min of flow rate. The model was satisfactorily fitted to the three breakthrough curves and the parameters were: axial dispersion coefficient (0.0165 dm2/min), effective diffusivity inside the particle (dm2/min 0.0884) and external transfer coefficient mass (0.45 dm/min). The breakthrough curve for CO2 in the process of pyrolysis of the sludge showed a fast saturation with traces of aerosols presents in the gas phase into the fixed bed under the reaction process

Solute transport in crystalline rocks at Äspö — I: Geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted A¨ spo¨ Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours–days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption Kds are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.

Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish A¨ spo¨ Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na+, 85Sr2 +, 47Ca2 +and more strongly sorbing 86Rb+, 133Ba2 +, 137Cs+. Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted Kd values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning Kds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.

Determinación de parámetros hidráulicos en columnas experimentales de suelos del sudeste de la provincia de Buenos Aires

La cuenca del Río Quequén Grande constituye una región representativa de los sistemas agrícolas pampeanos. El acuífero en esta región es la principal fuente de abastecimiento de agua para todos los usos, lo que pone de manifiesto la importancia de desarrollar los elementos necesarios para el estudio de la movilidad de los contaminantes a través de los suelos hacia la zona saturada. El objetivo del presente trabajo fue obtener parámetros hidráulicos bajo condiciones de saturación en columnas intactas de distintos tipos de suelos de la región, contemplando suelos de uso agrícola y natural. Se estudiaron dos series de suelo: Serie Azul y Serie Semillero Buck, con muestras tomadas en dos sectores próximos a las localidades de Lobería y La Dulce, respectivamente. Las columnas fueron eluidas con una solución acuosa de iones cloruro y se midió el incremento progresivo de conductividad. A partir de las curvas de arribo obtenidas, se ajustó la ecuación de transporte, obteniendo los valores de dispersividad y porosidad. Dichos parámetros no presentaron diferencias significativas según el uso para los suelos de Lobería, de textura franco arcillosa, indicando que la labranza no produciría modificaciones evidentes en el perfil del suelo. En La Dulce, donde los suelos poseen una textura más arenosa, los parámetros variaron según el uso del suelo, presentando una mayor dispersividad y menor porosidad del suelo bajo explotación agrícola extensiva.

Migração de solutos em basalto fraturado: quantificação experimental em laboratório e validação matemática

A avaliação do risco a contaminação e a escolha de técnicas de remediação de poluentes em aquíferos fraturados depende da quantificação dos fenômenos envolvidos no transporte de solutos. A geometria da fratura, usualmente caracterizada pela abertura, é o principal parâmetro que indiretamente controla o transporte nos aquíferos fraturados. A simplificação mais comum desse problema é assumir que as fraturas são um par de placas planas e paralelas, isto é, com uma abertura constante. No entanto, por causa do limitado número de trabalhos experimentais, não está esclarecida a adequabilidade do uso de uma abertura constante para simular o transporte conservativo em fraturas do Aquífero Serra Geral (ASG), Brasil. O objetivo deste trabalho é avaliar a influência da abertura de uma fratura natural do Aquífero Serra Geral sob o transporte conservativo de solutos. Uma amostra natural de basalto fraturado foi usada em um experimento hidráulico e de transporte de um traçador conservativo (escala de laboratório). O campo de abertura foi medido usando a técnica avançada, de alta resolução e tridimensional, chamada microtomografia computadorizada de raios-X. A concentração de traçador medida foi utilizada para validar uma solução analítica unidimensional da Equação de Advecção-dispersão (ADE). O desemprenho do ajuste da ADE às curvas de passagem experimentais foi avaliado para quatro diferentes tipos de aberturas constantes. Os resultados mostraram que o escoamento de água e o transporte de contaminantes pode ocorrer através de fraturas micrométricas, ocasionando, eventualmente, a contaminação do ASG. A abertura de balanço de massa é a única que pode ser chamada propriamente de \"abertura equivalente\". O uso de aberturas constantes na ADE não permitiu representar completamente o formato das curvas de passagem porque o campo de velocidade não é uniforme e intrinsicamente bidimensional. Portanto, na simulação do transporte deve-se incorporar a heterogeneidade da abertura da fratura.

Results of Closed-Flow Column Experiments (zip-archive 2.3 MB)

The identification of transport parameters by inverse modeling often suffers from equifinality or parameter correlation when models are fitted to observations of the solute breakthrough in column outflow experiments. This parameters uncertainty can be approached by the application of multiple experimental designs such as column experiments in open-flow mode and the recently proposed closed-flow mode. Latter are characterized by the recirculation of the column effluent into the solution supply vessel that feeds the inflow. Depending on the experimental conditions, the solute concentration in the solution supply vessel and the effluent follows a damped sinusoidal oscillation. As a result, the closed-flow experiment provides additional observables in the breakthrough curve. The evaluation of these emergent features allows intrinsic control over boundary conditions and impacts the uncertainty of parameters in inverse modeling. We present a comprehensive sensitivity analysis to illustrate the potential application of closed-flow experiments. We show that the sensitivity with respect to the apparent dispersion can be controlled by the experimenter leading to a decrease in parameter uncertainty as compared to classical experiments by an order of magnitude for optimal settings. With these finding we are also able to reduce the equifinality found for situations, where rate-limited interactions impede a proper determination of the apparent dispersion and rate coefficients. Furthermore, we show the expected breakthrough curve for equilibrium and kinetic sorption, the latter showing strong similarities to the behavior found for completely mixed batch reactor experiments. This renders the closed-flow mode a useful complementary approach to classical column experiments.

Computational adjustment technique for digital mammographic images based on the digitizer characteristic curve

We evaluated the performance of a novel procedure for segmenting mammograms and detecting clustered microcalcifications in two types of image sets obtained from digitization of mammograms using either a laser scanner, or a conventional ""optical"" scanner. Specific regions forming the digital mammograms were identified and selected, in which clustered microcalcifications appeared or not. A remarkable increase in image intensity was noticed in the images from the optical scanner compared with the original mammograms. A procedure based on a polynomial correction was developed to compensate the changes in the characteristic curves from the scanners, relative to the curves from the films. The processing scheme was applied to both sets, before and after the polynomial correction. The results indicated clearly the influence of the mammogram digitization on the performance of processing schemes intended to detect microcalcifications. The image processing techniques applied to mammograms digitized by both scanners, without the polynomial intensity correction, resulted in a better sensibility in detecting microcalcifications in the images from the laser scanner. However, when the polynomial correction was applied to the images from the optical scanner, no differences in performance were observed for both types of images. (C) 2008 SPIE and IS&T [DOI: 10.1117/1.3013544]

Rotation curve bifurcations as indicators of close recent galaxy encounters

Context. Rotation curves of interacting galaxies often show that velocities are either rising or falling in the direction of the companion galaxy. Aims. We seek to reproduce and analyse these features in the rotation curves of simulated equal-mass galaxies suffering a one-to-one encounter as possible indicators of close encounters. Methods. Using simulations of major mergers in 3D, we study the time evolution of these asymmetries in a pair of galaxies during the first passage. Results. Our main results are: (a) the rotation curve asymmetries appear right at pericentre of the first passage, (b) the significant disturbed rotation velocities occur within a small time interval, of similar to 0.5 Gyr h(-1), and, therefore, the presence of bifurcation in the velocity curve could be used as an indicator of the pericentre occurrence. These results are in qualitative agreement with previous findings for minor mergers and flybys.

Analysis of the E-J Curve of HTS Tapes Under DC and AC Magnetic Fields at 77 K

The evaluation of the electrical characteristics of technical HTS tapes are of the key importance in determining the design and operational features of superconducting power apparatuses as well as to understand the external factors which affect the superconducting performance. In this work we report the systematic measurements of the electric field versus current density, E-J relation of short samples for three commercial HTS tapes (BSCCO-2223 tapes, with and without steel reinforcement, and YBCO-coated conductor) at 77 K. In order to get sensitive and noiseless voltage signals the measurements were carried out with DC transport current and subjecting the broad surface tape to DC (0-300 mT) and AC (0-62 mT, 60 Hz) magnetic fields. The voltage is measured by a sensitive nanovoltmeter and the applied magnetic field is monitored by a Hall sensor placed on the tape broad surface. The comparison between the results obtained from the three tapes was done by fitting a power-law equation for currents in the vicinity of the critical current. For the current regime below the critical one a linear correlation of the electric field against the current density is observed. The BSCCO samples presented the same behavior, i.e., a decreasing of n-index with the increasing DC and AC magnetic field strength. Under AC field the decreasing slope of n-index is steeper as compared to DC field. The n-index curve for the YBCO tape showed similar behavior for AC field, however under DC field in the 0-390 mT range exhibited a slight decreasing of the n-index.

Right adolescent idiopathic thoracic curve (Lenke 1 A and B): does cost of instrumentation and implant density improve radiographic and cosmetic parameters?

In adolescent idiopathic scoliosis (AIS) there has been a shift towards increasing the number of implants and pedicle screws, which has not been proven to improve cosmetic correction. To evaluate if increasing cost of instrumentation correlates with cosmetic correction using clinical photographs. 58 Lenke 1A and B cases from a multicenter AIS database with at least 3 months follow-up of clinical photographs were used for analysis. Cosmetic parameters on PA and forward bending photographs included angular measurements of trunk shift, shoulder balance, rib hump, and ratio measurements of waist line asymmetry. Pre-op and follow-up X-rays were measured for coronal and sagittal deformity parameters. Cost density was calculated by dividing the total cost of instrumentation by the number of vertebrae being fused. Linear regression and spearman`s correlation were used to correlate cost density to X-ray and photo outcomes. Three independent observers verified radiographic and cosmetic parameters for inter/interobserver variability analysis. Average pre-op Cobb angle and instrumented correction were 54A degrees (SD 12.5) and 59% (SD 25) respectively. The average number of vertebrae fused was 10 (SD 1.9). The total cost of spinal instrumentation ranged from $6,769 to $21,274 (Mean $12,662, SD $3,858). There was a weak positive and statistically significant correlation between Cobb angle correction and cost density (r = 0.33, p = 0.01), and no correlation between Cobb angle correction of the uninstrumented lumbar spine and cost density (r = 0.15, p = 0.26). There was no significant correlation between all sagittal X-ray measurements or any of the photo parameters and cost density. There was good to excellent inter/intraobserver variability of all photographic parameters based on the intraclass correlation coefficient (ICC 0.74-0.98). Our method used to measure cosmesis had good to excellent inter/intraobserver variability, and may be an effective tool to objectively assess cosmesis from photographs. Since increasing cost density only improves mildly the Cobb angle correction of the main thoracic curve and not the correction of the uninstrumented spine or any of the cosmetic parameters, one should consider the cost of increasing implant density in Lenke 1A and B curves. In the area of rationalization of health care expenses, this study demonstrates that increasing the number of implants does not improve any relevant cosmetic or radiographic outcomes.

Motor development curve from 0 to 12 months in infants born preterm

Aim: To trace a reference curve for motor development from birth up to 12 months of corrected chronological age in infants born preterm and low birth weight. Methods: This is a cross-sectional study with a sample of 308 preterm infants (53% boys) weighing < 2500 g at birth. The Alberta Infant Motor Scale (AIMS) was used for motor development assessment. Results: Comparing the motor performance of preterm infants with infants from a standardized sample on the AIMS, it was found that, except for the age group of the newborn, preterm infants showed lower motor development scores in comparison with the AIMS normative sample in all age groups between 1 and 12 months. The curve of motor development showed a continuous increase in the number of motor skills of preterm infants during their first 12 months of age. However, the average of motor acquisitions of preterm infants showed a nonlinear pattern with a standard indicator of stabilization between 8 and 10 months of age. Conclusion: Preterm infants, 1-12 months of age, showed motor development AIMS scores lower than the standards established in the normative sample. The findings may contribute as norm-reference for assessing the motor development of preterm infants in follow-up programmes in developing countries.