Mixing-cell boundary conditions and apparent mass balance errors for advective-dispersive solute transport

In many field or laboratory situations, well-mixed reservoirs like, for instance, injection or detection wells and gas distribution or sampling chambers define boundaries of transport domains. Exchange of solutes or gases across such boundaries can occur through advective or diffusive processes. First we analyzed situations, where the inlet region consists of a well-mixed reservoir, in a systematic way by interpreting them in terms of injection type. Second, we discussed the mass balance errors that seem to appear in case of resident injections. Mixing cells (MC) can be coupled mathematically in different ways to a domain where advective-dispersive transport occurs: by assuming a continuous solute flux at the interface (flux injection, MC-FI), or by assuming a continuous resident concentration (resident injection). In the latter case, the flux leaving the mixing cell can be defined in two ways: either as the value when the interface is approached from the mixing-cell side (MC-RT -), or as the value when it is approached from the column side (MC-RT +). Solutions of these injection types with constant or-in one case-distance-dependent transport parameters were compared to each other as well as to a solution of a two-layer system, where the first layer was characterized by a large dispersion coefficient. These solutions differ mainly at small Peclet numbers. For most real situations, the model for resident injection MC-RI + is considered to be relevant. This type of injection was modeled with a constant or with an exponentially varying dispersion coefficient within the porous medium. A constant dispersion coefficient will be appropriate for gases because of the Eulerian nature of the usually dominating gaseous diffusion coefficient, whereas the asymptotically growing dispersion coefficient will be more appropriate for solutes due to the Lagrangian nature of mechanical dispersion, which evolves only with the fluid flow. Assuming a continuous resident concentration at the interface between a mixing cell and a column, as in case of the MC-RI + model, entails a flux discontinuity. This flux discontinuity arises inherently from the definition of a mixing cell: the mixing process is included in the balance equation, but does not appear in the description of the flux through the mixing cell. There, only convection appears because of the homogeneous concentration within the mixing cell. Thus, the solute flux through a mixing cell in close contact with a transport domain is generally underestimated. This leads to (apparent) mass balance errors, which are often reported for similar situations and erroneously used to judge the validity of such models. Finally, the mixing cell model MC-RI + defines a universal basis regarding the type of solute injection at a boundary. Depending on the mixing cell parameters, it represents, in its limits, flux as well as resident injections. (C) 1998 Elsevier Science B.V. All rights reserved.

Combined effect of heterogeneity, anisotropy and saturation on steady state flow and transport: Structure recognition and numerical simulation

1 Natural soil profiles may be interpreted as an arrangement of parts which are characterized by properties like hydraulic conductivity and water retention function. These parts form a complicated structure. Characterizing the soil structure is fundamental in subsurface hydrology because it has a crucial influence on flow and transport and defines the patterns of many ecological processes. We applied an image analysis method for recognition and classification of visual soil attributes in order to model flow and transport through a man-made soil profile. Modeled and measured saturation-dependent effective parameters were compared. We found that characterizing and describing conductivity patterns in soils with sharp conductivity contrasts is feasible. Differently, solving flow and transport on the basis of these conductivity maps is difficult and, in general, requires special care for representation of small-scale processes.

Do optic nerve head and visual field parameters in patients with obstructive sleep apnea syndrome differ from those in control individuals?

BACKGROUND It has been suggested that sleep apnea syndrome may play a role in normal-tension glaucoma contributing to optic nerve damage. The purpose of this study was to evaluate if optic nerve and visual field parameters in individuals with sleep apnea syndrome differ from those in controls. PATIENTS AND METHODS From the records of the sleep laboratory at the University Hospital in Bern, Switzerland, we recruited consecutive patients with severe sleep apnea syndrome proven by polysomnography, apnea-hypopnea index >20, as well as no sleep apnea controls with apnea-hypopnea index <10. Participants had to be unknown to the ophtalmology department and had to have no recent eye examination in the medical history. All participants underwent a comprehensive eye examination, scanning laser polarimetry (GDx VCC, Carl Zeiss Meditec, Dublin, California), scanning laser ophthalmoscopy (Heidelberg Retina Tomograph II, HRT II), and automated perimetry (Octopus 101 Programm G2, Haag-Streit Diagnostics, Koeniz, Switzerland). Mean values of the parameters of the two groups were compared by t-test. RESULTS The sleep apnea group consisted of 69 eyes of 35 patients; age 52.7 ± 9.7 years, apnea-hypopnea index 46.1 ± 24.8. As controls served 38 eyes of 19 patients; age 45.8 ± 11.2 years, apnea-hypopnea index 4.8 ± 1.9. A difference was found in mean intraocular pressure, although in a fully overlapping range, sleep apnea group: 15.2 ± 3.1, range 8-22 mmHg, controls: 13.6 ± 2.3, range 9-18 mmHg; p<0.01. None of the extended visual field, optic nerve head (HRT) and retinal nerve fiber layer (GDx VCC) parameters showed a significant difference between the groups. CONCLUSION Visual field, optic nerve head, and retinal nerve fiber layer parameters in patients with sleep apnea did not differ from those in the control group. Our results do not support a pathogenic relationship between sleep apnea syndrome and glaucoma.

Solution structure and interaction of cupiennin 1a, a spider venom peptide, with phospholipid bilayers

The solution structure of cupiennin 1a, a 35 residue, basic antibacterial peptide isolated from the venom of the spider Cupiennius salei, has been determined by nuclear magnetic resonance (NMR) spectroscopy. The peptide was found to adopt a helix−hinge−helix structure in a membrane mimicking solvent. The hinge may play a role in allowing the amphipathic N-terminal helix and polar C-terminal helix to orient independently upon membrane binding, in order to achieve maximal antibacterial efficacy. Solid-state 31P and 2H NMR was used to further study the effects of cupiennin 1a on the dynamic properties of lipid membranes, using zwitterionic chain deuterated dimyristoylphosphatidylcholine (d54-DMPC) and anionic dimyristoylphosphatidylglycerol (DMPG) multilamellar vesicles. In d54-DMPC alone, cupiennin 1a caused a decrease in the 31P chemical shift anisotropy, indicating some interaction with the lipid head groups, and a decrease in order over the entire acyl chain. In contrast, for the mixed (d54-DMPC/DMPG) lipid system cupiennin 1a appeared to induce lateral separation of the two lipids as evidenced by the 31P spectra, in which the peptide preferentially interacted with DMPG. Little effect was observed on the deuterated acyl chain order parameters in the d54-DMPC/DMPG model membranes. Furthermore, 31P NMR relaxation measurements confirmed a differential effect on the lipid motions depending upon the membrane composition. Therefore, subtle differences are likely in the mechanism by which cupiennin 1a causes membrane lysis in either prokaryotic or eukaryotic cells, and may explain the specific spectrum of activity.

Unique mixed phenotype and unexpected functional effect revealed by novel compound heterozygosity mutations involving SCN5A.

BACKGROUND Functional characterization of mutations involving the SCN5A-encoded cardiac sodium channel has established the pathogenic mechanisms for type 3 long QT syndrome and type 1 Brugada syndrome and has provided key insights into the physiological importance of essential structure-function domains. OBJECTIVE This study sought to present the clinical and biophysical phenotypes discerned from compound heterozygosity mutations in SCN5A on different alleles in a toddler diagnosed with QT prolongation and fever-induced ventricular arrhythmias. METHODS A 22-month-old boy presented emergently with fever and refractory ventricular tachycardia. Despite restoration of sinus rhythm, the infant sustained profound neurological injury and died. Using polymerase chain reaction, denaturing high-performance liquid chromatography, and direct DNA sequencing, comprehensive open-reading frame/splice mutational analysis of the 12 known long QT syndrome susceptibility genes was performed. RESULTS The infant had 2 SCN5A mutations: a maternally inherited N-terminal frame shift/deletion (R34fs/60) and a paternally inherited missense mutation, R1195H. The mutations were engineered by site-directed mutagenesis and heterologously expressed transiently in HEK293 cells. As expected, the frame-shifted and prematurely truncated peptide, SCN5A-R34fs/60, showed no current. SCN5A-R1195H had normal peak and late current but abnormal voltage-dependent gating parameters. Surprisingly, co-expression of SCN5A-R34fs/60 with SCN5A-R1195H elicited a significant increase in late sodium current, whereas co-expression of SCN5A-WT with SCN5A-R34fs/60 did not. CONCLUSIONS A severe clinical phenotype characterized by fever-induced monomorphic ventricular tachycardia and QT interval prolongation emerged in a toddler with compound heterozygosity involving SCN5A: R34fs/60, and R1195H. Unexpectedly, the 94-amino-acid fusion peptide derived from the R34fs/60 mutation accentuated the late sodium current of R1195H-containing Na(V)1.5 channels in vitro.

Microscopic theory of cooperative spin crossover: Interaction of molecular modes with phonons

In this article, we present a new microscopic theoretical approach to the description of spin crossover in molecular crystals. The spin crossover crystals under consideration are composed of molecular fragments formed by the spin-crossover metal ion and its nearest ligand surrounding and exhibiting well defined localized (molecular) vibrations. As distinguished from the previous models of this phenomenon, the developed approach takes into account the interaction of spin-crossover ions not only with the phonons but also a strong coupling of the electronic shells with molecular modes. This leads to an effective coupling of the local modes with phonons which is shown to be responsible for the cooperative spin transition accompanied by the structural reorganization. The transition is characterized by the two order parameters representing the mean values of the products of electronic diagonal matrices and the coordinates of the local modes for the high- and low-spin states of the spin crossover complex. Finally, we demonstrate that the approach provides a reasonable explanation of the observed spin transition in the [Fe(ptz)6](BF4)2 crystal. The theory well reproduces the observed abrupt low-spin → high-spin transition and the temperature dependence of the high-spin fraction in a wide temperature range as well as the pronounced hysteresis loop. At the same time within the limiting approximations adopted in the developed model, the evaluated high-spin fraction vs. T shows that the cooperative spin-lattice transition proves to be incomplete in the sense that the high-spin fraction does not reach its maximum value at high temperature.

Triple isotope (δD, δ17O, δ18O) study on precipitation, drip water and speleothem fluid inclusions for a Western Central European cave (NW Switzerland)

Deuterium (δD) and oxygen (δ18O) isotopes are powerful tracers of the hydrological cycle and have been extensively used for paleoclimate reconstructions as they can provide information on past precipitation, temperature and atmospheric circulation. More recently, the use of δ17O excess derived from precise measurement of δ17O and δ18O gives new and additional insights in tracing the hydrological cycle whereas uncertainties surround this proxy. However, 17O excess could provide additional information on the atmospheric conditions at the moisture source as well as about fractionations associated with transport and site processes. In this paper we trace water stable isotopes (δD,δ17O and δ18O) along their path from precipitation to cave drip water and finally to speleothem fluid inclusions for Milandre cave in northwestern Switzerland. A two year-long daily resolved precipitation isotope record close to the cave site is compared to collected cave drip water (3 months average resolution) and fluid inclusions of modern and Holocene stalagmites. Amount weighted mean δD,δ18O and δ17O are -71.0‰, -9.9‰, -5.2‰ for precipitation, -60.3‰, -8.7‰, -4.6‰ for cave drip water and -61.3‰, -8.3‰, -4.7‰ for recent fluid inclusions respectively. Second order parameters have also been derived in precipitation and drip water and present similar values with 18 per meg for 17O excess whereas d-excess is 1.5‰ more negative in drip water. Furthermore, the atmospheric signal is shifted towards enriched values in the drip water and fluid inclusions (Δ of ~ + 10‰ for δD). The isotopic composition of cave drip water exhibits a weak seasonal signal which is shifted by around 8 - 10 months (groundwater residence time) when compared to the precipitation. Moreover, we carried out the first δ17O measurement in speleothem fluid inclusions, as well as the first comparison of the δ17 O behaviour from the meteoric water to the fluid inclusions entrapment in speleothems. This study on precipitation, drip water and fluid inclusions will be used as a speleothem proxy calibration for Milandre cave in order to reconstruct paleotemperatures and moisture source variations for Western Central Europe.

Threshold-dependent effects on peripheral nerve in vivo excitability properties in the rat

Various factors, including maturity, have been shown to influence peripheral nerve excitability measures, but little is known about differences in these properties between axons with different stimulation thresholds. Multiple nerve excitability tests were performed on the caudal motor axons of immature and mature female rats, recording from tail muscles at three target compound muscle action potential (CMAP) levels: 10%, 40% ("standard" level), and 60% of the maximum CMAP amplitude. Compared to lower target levels, axons at high target levels have the following characteristics: lower strength-duration time constant, less threshold reduction during depolarizing currents and greater threshold increase to hyperpolarizing currents, most notably to long hyperpolarizing currents in mature rats. Threshold-dependent effects on peripheral nerve excitability properties depend on the maturation stage, especially inward rectification (Ih), which becomes inversely related to threshold level. Performing nerve excitability tests at different target levels is useful in understanding the variation in membrane properties between different axons within a nerve. Because of the threshold effects on nerve excitability and the possibility of increased variability between axons and altered electric recruitment order in disease conditions, excitability parameters measured only at the "standard" target level should be interpreted with caution, especially the responses to hyperpolarizing currents.

Orientation dependent molecular friction on organic layer compound crystals

High resolution friction force maps of the benzylammonium terminated crystalline surface of a layer compound are presented. The lateral force map acquired with an atomic force microscope, reveals a significant contrast between different molecular orientations yielding molecular rows which differ from their neighboring ones. The single crystals are formed by stacks of copper oxalate sheets sandwiched between stereoregular organic cations, resulting in highly organized surface structures. Single molecular defects are observed at small loads. The experimental results are compared with numerical calculations which indicate a transition from an unperturbed state at small loads to a distorted state at higher loads. (C) 2011 American Institute of Physics.

Time-dependent interactions of the two porphyrinic compounds chlorin e6 and mono-L-aspartyl-chlorin e6 with phospholipid vesicles probed by NMR spectroscopy

The distribution processes of chlorin e6 (CE) and monoaspartyl-chlorin e6 (MACE) between the outer and inner phospholipid monolayers of 1,2-dioleoyl-phosphatidylcholine (DOPC) vesicles were monitored by 1H NMR spectroscopy through analysis of chemical shifts and line widths of the DOPC vesicle resonances. Chlorin adsorption to the outer vesicle monolayer induced changes in the DOPC 1H NMR spectrum. Most pronounced was a split of the N-methyl choline resonance, allowing for separate analysis of inner and outer vesicle layers. Transbilayer distribution of the chlorin compounds was indicated by time-dependent characteristic spectral changes of the DOPC resonances. Kinetic parameters for the flip-flop processes, that is, half-lives and rate constants, were obtained from the experimental data points. In comparison to CE, MACE transbilayer movement was significantly reduced, with MACE remaining more or less attached to the outer membrane layer. The distribution coefficients for CE and MACE between the vesicular and aqueous phase were determined. Both CE and MACE exhibited a high affinity for the vesicular phase. For CE, a positive correlation was found between transfer rate and increasing molar ratio CE/DOPC. Enhanced membrane rigidity induced by increasing amounts of cholesterol into the model membrane was accompanied by a decrease of CE flip-flop rates across the membrane. The present study shows that the movement of porphyrins across membranes can efficiently be investigated by 1H NMR spectroscopy and that small changes in porphyrin structure can have large effects on membrane kinetics.

Relative Predicativity and dependent recursion in second-order set theory and higher-orders theories

This article reports that some robustness of the notions of predicativity and of autonomous progression is broken down if as the given infinite total entity we choose some mathematical entities other than the traditional ω. Namely, the equivalence between normal transfinite recursion scheme and new dependent transfinite recursion scheme, which does hold in the context of subsystems of second order number theory, does not hold in the context of subsystems of second order set theory where the universe V of sets is treated as the given totality (nor in the contexts of those of n+3-th order number or set theories, where the class of all n+2-th order objects is treated as the given totality).

On the initiation of boudinage and folding instabilities in layered elasto-visco-plastic solids - Insights from natural microstructures and numerical simulations

The present understanding of the initiation of boudinage and folding structures is based on viscosity contrasts and stress exponents, considering an intrinsically unstable state of the layer. The criterion of localization is believed to be prescribed by geometry-material interactions, which are often encountered in natural structures. An alternative localization phenomenon has been established for ductile materials, in which instability emerges for critical material parameters and loading rates from homogeneous conditions. In this thesis, conditions are sought under which this type of instability prevails and whether localization in geological materials necessarily requires a trigger by geometric imperfections. The relevance of critical deformation conditions, material parameters and the spatial configuration of instabilities are discussed in a geological context. In order to analyze boudinage geometries, a numerical eigenmode analysis is introduced. This method allows determining natural frequencies and wavelengths of a structure and inducing perturbations on these frequencies. In the subsequent coupled thermo-mechanical simulations, using a grain size evolution and end-member flow laws, localization emerges when material softening through grain size sensitive viscous creep sets in. Pinch-and-swell structures evolve along slip lines through a positive feedback between the matrix response and material bifurcations inside the layer, independent from the mesh-discretization length scale. Since boudinage and folding are considered to express the same general instability, both structures should arise independently of the sign of the loading conditions and for identical material parameters. To this end, the link between material to energy instabilities is approached by means of bifurcation analyses of the field equations and finite element simulations of the coupled system of equations. Boudinage and folding structures develop at the same critical energy threshold, where dissipative work by temperature-sensitive creep overcomes the diffusive capacity of the layer. This finding provides basis for a unified theory for strain localization in layered ductile materials. The numerical simulations are compared to natural pinch-and-swell microstructures, tracing the adaption of grain sizes, textures and creep mechanisms in calcite veins. The switch from dislocation to diffusion creep relates to strain-rate weakening, which is induced by dissipated heat from grain size reduction, and marks the onset of continuous necking. The time-dependent sequence uncovers multiple steady states at different time intervals. Microstructurally and mechanically stable conditions are finally expressed in the pinch-and-swell end members. The major outcome of this study is that boudinage and folding can be described as the same coupled energy-mechanical bifurcation, or as one critical energy attractor. This finding allows the derivation of critical deformation conditions and fundamental material parameters directly from localized structures in the field.