A Pre-processing Moving Mesh Method for Discontinuous Galerkin Approximations of Advection-Diffusion-Reaction Problems

We propose a pre-processing mesh re-distribution algorithm based upon harmonic maps employed in conjunction with discontinuous Galerkin approximations of advection-diffusion-reaction problems. Extensive two-dimensional numerical experiments with different choices of monitor functions, including monitor functions derived from goal-oriented a posteriori error indicators are presented. The examples presented clearly demonstrate the capabilities and the benefits of combining our pre-processing mesh movement algorithm with both uniform, as well as, adaptive isotropic and anisotropic mesh refinement.

An hr-adaptive discontinuous Galerkin method for advection-diffusion problems

We propose an adaptive mesh refinement strategy based on exploiting a combination of a pre-processing mesh re-distribution algorithm employing a harmonic mapping technique, and standard (isotropic) mesh subdivision for discontinuous Galerkin approximations of advection-diffusion problems. Numerical experiments indicate that the resulting adaptive strategy can efficiently reduce the computed discretization error by clustering the nodes in the computational mesh where the analytical solution undergoes rapid variation.

Age distribution of the central stars of galactic disk planetary nebulae

Context. Determination of the ages of central stars of planetary nebulae (CSPN) is a complex problem, and there is presently no single method that can be generally applied. We have developed several methods of estimating the ages of CSPN, based on both the observed nebular properties and some properties of the stars themselves. Aims. Our aim is to estimate the ages and the age distribution of CSPN and to compare the derived results with mass and age determinations of CSPN and white dwarfs based on empirical determinations of these quantities. Methods. We considered a sample of planetary nebulae in the galactic disk, most of which (similar to 69%) are located in the solar neighbourhood, within 3 kpc from the Sun. We discuss several methods of deriving the age distribution of CSPN, namely; (i) the use of an age-metallicity relation that also depends on the galactocentric distance; (ii) the use of an age-metallicity relation obtained for the galactic disk; and (iii) the determination of ages from the central star masses obtained from the observed nitrogen abundances. Results. We estimated the age distribution of CSPN with average uncertainties of 1-2 Gyr, and compared our results with the expected distribution based both on the observed mass distribution of white dwarfs and on the age distribution derived from available mass distributions of CSPN. Based on our derived age distributions, we conclude that most CSPN in the galactic disk have ages under 6 Gyr, and that the age distribution is peaked around 2-4 Gyr.

Chemical similarities between Galactic bulge and local thick disk red giants: O, Na, Mg, Al, Si, Ca, and Ti

Context. The formation and evolution of the Galactic bulge and its relationship with the other Galactic populations is still poorly understood. Aims. To establish the chemical differences and similarities between the bulge and other stellar populations, we performed an elemental abundance analysis of alpha- (O, Mg, Si, Ca, and Ti) and Z-odd (Na and Al) elements of red giant stars in the bulge as well as of local thin disk, thick disk and halo giants. Methods. We use high-resolution optical spectra of 25 bulge giants in Baade's window and 55 comparison giants (4 halo, 29 thin disk and 22 thick disk giants) in the solar neighborhood. All stars have similar stellar parameters but cover a broad range in metallicity (-1.5 < [Fe/H] < +0.5). A standard 1D local thermodynamic equilibrium analysis using both Kurucz and MARCS models yielded the abundances of O, Na, Mg, Al, Si, Ca, Ti and Fe. Our homogeneous and differential analysis of the Galactic stellar populations ensured that systematic errors were minimized. Results. We confirm the well-established differences for [alpha/Fe] at a given metallicity between the local thin and thick disks. For all the elements investigated, we find no chemical distinction between the bulge and the local thick disk, in agreement with our previous study of C, N and O but in contrast to other groups relying on literature values for nearby disk dwarf stars. For -1.5 < [Fe/H] < -0.3 exactly the same trend is followed by both the bulge and thick disk stars, with a star-to-star scatter of only 0.03 dex. Furthermore, both populations share the location of the knee in the [alpha/Fe] vs. [Fe/H] diagram. It still remains to be confirmed that the local thick disk extends to super-solar metallicities as is the case for the bulge. These are the most stringent constraints to date on the chemical similarity of these stellar populations. Conclusions. Our findings suggest that the bulge and local thick disk stars experienced similar formation timescales, star formation rates and initial mass functions, confirming thus the main outcomes of our previous homogeneous analysis of [O/Fe] from infrared spectra for nearly the same sample. The identical a-enhancements of thick disk and bulge stars may reflect a rapid chemical evolution taking place before the bulge and thick disk structures we see today were formed, or it may reflect Galactic orbital migration of inner disk/bulge stars resulting in stars in the solar neighborhood with thick-disk kinematics.

State protection under collective damping and diffusion

In this paper we provide a recipe for state protection in a network of oscillators under collective damping and diffusion. Our strategy is to manipulate the network topology, i.e., the way the oscillators are coupled together, the strength of their couplings, and their natural frequencies, in order to create a relaxation-diffusion-free channel. This protected channel defines a decoherence-free subspace (DFS) for nonzero-temperature reservoirs. Our development also furnishes an alternative approach to build up DFSs that offers two advantages over the conventional method: it enables the derivation of all the network-protected states at once, and also reveals, through the network normal modes, the mechanism behind the emergence of these protected domains.

Solution of a heterogeneous modeling of a horizontal-flow anaerobic immobilized biomass (HAIB) reactor by the sequencing method

A modeling study was completed to develop a methodology that combines the sequencing and finite difference methods for the simulation of a heterogeneous model of a tubular reactor applied in the treatment of wastewater. The system included a liquid phase (convection diffusion transport) and a solid phase (diffusion reaction) that was obtained by completing a mass balance in the reactor and in the particle, respectively. The model was solved using a pilot-scale horizontal-flow anaerobic immobilized biomass (HAIB) reactor to treat domestic sewage, with the concentration results compared with the experimental data. A comparison of the behavior of the liquid phase concentration profile and the experimental results indicated that both the numerical methods offer a good description of the behavior of the concentration along the reactor. The advantage of the sequencing method over the finite difference method is that it is easier to apply and requires less computational time to model the dynamic simulation of outlet response of HAIB.

Development and validation of HPLC method for imiquimod determination in skin penetration studies

A simple, rapid and sensitive analytical procedure for the measurement of imiquimod in skin samples after in vitro penetration studies has been developed and validated. In vitro penetration studies were carried out in Franz diffusion cells with porcine skin. Tape stripping technique was used to separate the stratum corneum (SC) from the viable epidermis and dermis. Imiquimod was extracted from skin samples using a 7:3 (v/v) methanol:acetate buffer (100 mm, pH 4.0) solution and ultrasonication. Imiquimod was analyzed by H-PLC using C(8) column and UV detection at 242 ran. The mobile phase used was acetonitrile:acetate buffer (pH 4.0, 100 mM):diethylamine (30:69.85:0.15, v/v) with flow rate 1 mL/min. Imiquimod eluted at 4.1 min and the running time was limited to 6.0 min. The procedure was linear across the following concentration ranges: 100-2500 ng/mL for both SC and tape-stripped skin and 20-800 ng/mL for receptor solution. Intra-day and inter-day accuracy and precision values were lower than 20% at the limit of quantitation. The recovery values ranged from 80 to 100%. The method is adequate to assay imiquimod from skin samples, enabling the determination of the cutaneous penetration profile of uniquimod by in vitro studies. Copyright (C) 2008 John Wiley & Sons, Ltd.

Short-term spatial diffusion in sigma o+ o- optical molasses

We have measured the spatial diffusion of atoms in a three-dimensional sigma(+)-sigma(-) optical molasses over twenty milliseconds timescale, starting from the initial interaction of the atoms with the molasses. We find that the diffusion constants agree well with a linear model for these short time scales and also compare favourably to other studies of diffusion made over longer time scales. These measurements enable us to quantify the detection method known as freezing molasses. We discuss this method, for detecting and measuring the momentum distribution of cold atoms, which relies on the slow diffusion of atoms in optical molasses to produce a freeze-frame of the spatial distribution of the atoms. This method enables a longer interrogation interval, providing a greatly increased signal-to-noise ratio. (C) 1998 Elsevier Science B.V.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Loss of connectivity in Alzheimer's disease: an evaluation of white matter tract integrity with colour coded MR diffusion tensor imaging

A novel MRI method-diffusion tensor imaging-was used to compare the integrity of several white matter fibre tracts in patients with probable Alzheimer's disease. Relative to normal controls, patients with probable Alzheimer's disease showed a highly significant reduction in the integrity of the association white matter fibre tracts, such as the splenium of the corpus callosum, superior longitudinal fasciculus, and cingulum. By contrast, pyramidal tract integrity seemed unchanged. This novel finding is consistent with the clinical presentation of probable Alzheimer's disease, in which global cognitive decline is a more prominent feature than motor disturbance.

Self-assembled films from WO(3): Electrochromism and lithium ion diffusion

WO(3)/chitosan and WO(3)/chitosan/poly(ethylene oxide) (PEO) films were prepared by the layer-by-layer method. The presence of chitosan enabled PEO to be carried into the self-assembled structure, contributing to an increase in the Li(+) diffusion rate. On the basis of the galvanostatic intermittent titration technique (GITT) and the quadratic logistic equation (QLE), a spectroelectrochemical method was used for determination of the ""optical"" diffusion coefficient (D(op)), enabling analysis of the Li(+) diffusion rate and, consequently, the coloration front rate in these host matrices. The D(op) values within the WO(3)/chitosan/PEO film were significantly higher than those within the WO(3)/chitosan film, mainly for higher values of injected charge. The presence of PEO also ensured larger accessibility to the electroactive sites, in accordance with the method employed here. Hence, this spectroelectrochemical method allowed us to separate the contribution of the diffusion process from the number of accessible electroactive sites in the materials, thereby aiding a better understanding of the useful electrochemical and electrochromic properties of these films for use in electrochromic devices. (C) 2010 Elsevier B.V. All rights reserved.

The states, diffusion, and concentration distribution of water in radiation-formed PVA/PVP hydrogels

Hydrogels with various compositions of polyvinyl alcohol (PVA) and poly(1-vinyl-2-pyrrolidinone) (PVP) were prepared by irradiating mixtures of PVA and PVP in aqueous solutions with gamma-rays from Co-60 sources at room temperature. The states of water in the hydrogels were characterized using DSC and NMR T-2 relaxation measurements and the kinetics of water diffusion in the hydrogels were studied by sorption experiments and NMR imaging. The DSC endothermic peaks in the temperature range -10 to +10 degrees C implied that there are at least two kinds of freezable water present in the matrix. The difference between the total water content and the freezable water content was refer-red to as bound water, which is not freezable. The weight fraction of water at which only nonfreezable water is present in a hydrogel with F-VP = 0.19 has been estimated to be g(H2O)/g(Polymer) = 0.375. From water sorption experiments, it was demonstrated that the early stage of the diffusion of water into the hydrogels was Fickian. A curve-fit of the early-stage experimental data to the Fickian model allowed determination of the water diffusion coefficient, which was found to lie between 1.5 x 10(-11) m(2) s(-1) and 4.5 x 10(-11) m(2) s(-1), depending on the polymer composition, the cross-link density, and the temperature. It was also found that the energy barrier for diffusion of water molecules into PVA/PVP hydrogels was approximate to 24 kJ mol(-1). Additionally, the diffusion coefficients determined from NMR imaging of the volumetric swelling of the gels agreed well with the results obtained by the mass sorption method.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

Measuring macromolecular diffusion using heteronuclear multiple-quantum pulsed-field-gradient NMR

We have previously shown that H-1 pulsed-field-gradient (PFG) NMR spectroscopy provides a facile method for monitoring protein self-association and can be used, albeit with some caveats, to measure the apparent molecular mass of the diffusant [Dingley et al. (1995) J. Biomol. NMR, 6, 321-328]. In this paper we show that, for N-15-labelled proteins, selection of H-1-N-15 multiple-quantum (MQ) coherences in PFG diffusion experiments provides several advantages over monitoring H-1 single-quantum (SQ) magnetization. First, the use of a gradient-selected MQ filter provides a convenient means of suppressing resonances from both the solvent and unlabelled solutes. Second, H-1-N-15 zero-quantum coherence dephases more rapidly than H-1 SQ coherence under the influence of a PFG. This allows the diffusion coefficients of larger proteins to be measured more readily. Alternatively, the gradient length and/or the diffusion delay may be decreased, thereby reducing signal losses from relaxation. In order to extend the size of macromolecules to which these experiments can be applied, we have developed a new MQ PFG diffusion experiment in which the magnetization is stored as longitudinal two-spin order for most of the diffusion period, thus minimizing sensitivity losses due to transverse relaxation and J-coupling evolution.

Evaluation of the buccal vestibule-palatal diffusion of 4% articaine hydrochloride in impacted maxillary third molar extractions

Aims: The aim of this study was to evaluate the vestibular-palatal diffusion of 4% articaine with epinephrine 1: 100,000 and 1: 200,000, in impacted maxillary third molar extractions, without palatal injection. Materials and Method: Two hundred teeth were selected from patients age 15 to 46. Patients were divided into 4 groups: 1A, were anesthetized with 4% articaine 1: 100,000 and the surgery was initiated 5 minutes following anesthesia. 1B, used 4% articaine 1: 100,000 but the surgery was started 10 minutes after anesthesia. 2A, used 4% articaine 1: 200,000 the surgery was started 5 minutes after. 2B, used 4% articaine 1: 200,000 but 10 minutes was allowed for anesthetic diffusion before the initiation of in groups (50 extractions each) only buccal vestibule anesthesia was initially administered (i.e. no palatal injections were used). Results: The rate of sufficient vestibule-palatal diffusion, as determined by the lack of necessity of supplemental palatal anesthesia, was: 1A(84%), 1B(98%), 2A(78%), 2B(82%). Chi-square (X2) and residual analyses showed that a higher vestibule-palatal diffusion was obtained using 4% articaine 1: 100,000 with a period of 10 minutes (p<0.05). Conclusions: Most of the extractions could be performed only with vestibule anesthesia. However, vasoconstrictor concentration and the time interval between administration of the anesthetic and initiation of surgery did influence buccal vestibule-palatal diffusion of 4% articaine in the extraction models used.