An adaptive multiscale method for density-driven instabilities

Accurate modeling of flow instabilities requires computational tools able to deal with several interacting scales, from the scale at which fingers are triggered up to the scale at which their effects need to be described. The Multiscale Finite Volume (MsFV) method offers a framework to couple fine-and coarse-scale features by solving a set of localized problems which are used both to define a coarse-scale problem and to reconstruct the fine-scale details of the flow. The MsFV method can be seen as an upscaling-downscaling technique, which is computationally more efficient than standard discretization schemes and more accurate than traditional upscaling techniques. We show that, although the method has proven accurate in modeling density-driven flow under stable conditions, the accuracy of the MsFV method deteriorates in case of unstable flow and an iterative scheme is required to control the localization error. To avoid large computational overhead due to the iterative scheme, we suggest several adaptive strategies both for flow and transport. In particular, the concentration gradient is used to identify a front region where instabilities are triggered and an accurate (iteratively improved) solution is required. Outside the front region the problem is upscaled and both flow and transport are solved only at the coarse scale. This adaptive strategy leads to very accurate solutions at roughly the same computational cost as the non-iterative MsFV method. In many circumstances, however, an accurate description of flow instabilities requires a refinement of the computational grid rather than a coarsening. For these problems, we propose a modified iterative MsFV, which can be used as downscaling method (DMsFV). Compared to other grid refinement techniques the DMsFV clearly separates the computational domain into refined and non-refined regions, which can be treated separately and matched later. This gives great flexibility to employ different physical descriptions in different regions, where different equations could be solved, offering an excellent framework to construct hybrid methods.

Data-driven clustering : new methods and applications

Abstract : This work is concerned with the development and application of novel unsupervised learning methods, having in mind two target applications: the analysis of forensic case data and the classification of remote sensing images. First, a method based on a symbolic optimization of the inter-sample distance measure is proposed to improve the flexibility of spectral clustering algorithms, and applied to the problem of forensic case data. This distance is optimized using a loss function related to the preservation of neighborhood structure between the input space and the space of principal components, and solutions are found using genetic programming. Results are compared to a variety of state-of--the-art clustering algorithms. Subsequently, a new large-scale clustering method based on a joint optimization of feature extraction and classification is proposed and applied to various databases, including two hyperspectral remote sensing images. The algorithm makes uses of a functional model (e.g., a neural network) for clustering which is trained by stochastic gradient descent. Results indicate that such a technique can easily scale to huge databases, can avoid the so-called out-of-sample problem, and can compete with or even outperform existing clustering algorithms on both artificial data and real remote sensing images. This is verified on small databases as well as very large problems. Résumé : Ce travail de recherche porte sur le développement et l'application de méthodes d'apprentissage dites non supervisées. Les applications visées par ces méthodes sont l'analyse de données forensiques et la classification d'images hyperspectrales en télédétection. Dans un premier temps, une méthodologie de classification non supervisée fondée sur l'optimisation symbolique d'une mesure de distance inter-échantillons est proposée. Cette mesure est obtenue en optimisant une fonction de coût reliée à la préservation de la structure de voisinage d'un point entre l'espace des variables initiales et l'espace des composantes principales. Cette méthode est appliquée à l'analyse de données forensiques et comparée à un éventail de méthodes déjà existantes. En second lieu, une méthode fondée sur une optimisation conjointe des tâches de sélection de variables et de classification est implémentée dans un réseau de neurones et appliquée à diverses bases de données, dont deux images hyperspectrales. Le réseau de neurones est entraîné à l'aide d'un algorithme de gradient stochastique, ce qui rend cette technique applicable à des images de très haute résolution. Les résultats de l'application de cette dernière montrent que l'utilisation d'une telle technique permet de classifier de très grandes bases de données sans difficulté et donne des résultats avantageusement comparables aux méthodes existantes.

Moving from a Mechanical Microkeratome to a Femtosecond Laser for LASIK to Correct Astigmatic Patients: Clinical Outcomes of a Retrospective, Consecutive, Comparative Study.

Background: To evaluate outcomes after optimized laser in situ keratomileusis (LASIK) for astigmatism correction with flap created by a mechanical microkeratome or a femtosecond laser. Patients and Methods: In this retrospective study, a total of 102 eyes of 71 consecutive patients were enrolled undergoing optimized LASIK treatments using the Allegretto laser system (WaveLight Laser Technologie AG, Erlangen, Germany). A mechanical microkeratome for flap creation was used (One Use, Moria®) in 46 eyes (31 patients, spherical equivalent [SE] -4.44 D ± 2.4) and a femtosecond laser was used (LDV, Ziemer®) in 56 eyes (40 patients, spherical equivalent [SE] -3.07 D ± 3.3). The two groups were matched for inclusion criteria and were operated under similar conditions by the same surgeon. Results: Overall, the preoperative spherical equivalent was -9.5 diopters (D) to +3.37 D; the preoperative manifest astigmatism was between -1.5 D and -3.5 D. At 6 months postoperatively, the mean postoperative uncorrected distance visual acuity (UDVA) was 0.93 ± 0.17 (range 0.4 to 1.2) in the Moria group and 1.0 ± 0.21 (range 0.6 to 1.6) in the Femto group, which was statistically significant (p = 0.003). Comparing the cylinder power there was a statistical difference between the two groups (p = 0.0015). Conclusions: This study shows that the method of flap creation has a significant impact on postoperative astigmatism with a significantly better postoperative UDVA in the Femto group. These findings suggest that the femtosecond laser provides a better platform for LASIK treatment of astigmatism than the commonly used microkeratome.

Mean first-passage time of continuous non-Markovian processes driven by colored noise

An equation for mean first-passage times of non-Markovian processes driven by colored noise is derived through an appropriate backward integro-differential equation. The equation is solved in a Bourret-like approximation. In a weak-noise bistable situation, non-Markovian effects are taken into account by an effective diffusion coefficient. In this situation, our results compare satisfactorily with other approaches and experimental data.

Dynamics of a paramagnetic colloidal particle driven on a magnetic-bubble lattice

We present a theoretical study of the recently observed dynamical regimes of paramagnetic colloidal particles externally driven above a regular lattice of magnetic bubbles [P. Tierno, T. H. Johansen, and T. M. Fischer, Phys. Rev. Lett. 99, 038303 (2007)]. An external precessing magnetic field alters the potential generated by the surface of the film in such a way to either drive the particle circularly around one bubble, ballistically through the array, or in triangular orbits on the interstitial regions between the bubbles. In the ballistic regime, we observe different trajectories performed by the particles phase locked with the external driving. Superdiffusive motion, which was experimentally found bridging the localized and delocalized dynamics, emerge only by introducing a certain degree of randomness into the bubbles size distribution.

Trace element partitioning in HP-LT metamorphic assemblages during subduction-related metamorphism, Ile de Groix, France: a detailed LA-ICPMS study

Devolatilization reactions and subsequent transfer of fluid from subducted oceanic crust into the overlying mantle wedge are important processes, which are responsible for the specific geochemical characteristics of subduction-related metamorphic rocks, as well as those of arc magmatism. To better understand the geochemical fingerprint induced by fluid mobilization during dehydration and rehydration processes related to subduction zone metamorphism, the trace element and rare earth element (REE) distribution patterns in HP-LT metamorphic assemblages in eclogite-, blueschist- and greenschist-facies rocks of the Ile de Groix were obtained by laser ablation inductively coupled plasma mass spectrometry (LA-ICPMS) analysis. This study focuses on 10 massive basic rocks representing former hydrothermally altered mid-ocean ridge basalts (MORB), four banded basic rocks of volcano-sedimentary origin and one micaschist. The main hosts for incompatible trace elements are epidote (REE, Th, U, Pb, Sr), garnet [Y, heavy REE (HREE)], phengite (Cs, Rb, Ba, B), titanite [Ti, Nb, Ta, REE; HREE > LREE (light REE)], rutile (Ti, Nb, Ta) and apatite (REE, Sr). The trace element contents of omphacite, amphibole, albite and chlorite are low. The incompatible trace element contents of minerals are controlled by the stable metamorphic mineral assemblage and directly related to the appearance, disappearance and reappearance of minerals, especially epidote, garnet, titanite, rutile and phengite, during subduction zone metamorphism. Epidote is a key mineral in the trace element exchange process because of its large stability field, ranging from lower greenschist- to blueschist- and eclogite-facies conditions. Different generations of epidote are generally observed and related to the coexisting phases at different stages of the metamorphic cycle (e.g. lawsonite, garnet, titanite). Epidote thus controls most of the REE budget during the changing P-T conditions along the prograde and retrograde path. Phengite also plays an important role in determining the large ion lithophile element (LILE) budget, as it is stable to high P-T conditions. The breakdown of phengite causes the release of LILE during retrogression. A comparison of trace element abundances in whole-rocks and minerals shows that the HP-LT metamorphic rocks largely retain the geochemical characteristics of their basic, volcano-sedimentary and pelitic protoliths, including a hydrothermal alteration overprint before the subduction process. A large part of the incompatible trace elements remained trapped in the rocks and was recycled within the various metamorphic assemblages stable under changing metamorphic conditions during the subduction process, indicating that devolatilization reactions in massive basic rocks do not necessarily imply significant simultaneous trace element and REE release.

Characterization of surface functional groups present on laboratory-generated and ambient aerosol particles by means of heterogeneous titration reactions

A Knudsen flow reactor has been used to quantify surface functional groups on aerosols collected in the field. This technique is based on a heterogeneous titration reaction between a probe gas and a specific functional group on the particle surface. In the first part of this work, the reactivity of different probe gases on laboratory-generated aerosols (limonene SOA, Pb(NO3)2, Cd(NO3)2) and diesel reference soot (SRM 2975) has been studied. Five probe gases have been selected for the quantitative determination of important functional groups: N(CH3)3 (for the titration of acidic sites), NH2OH (for carbonyl functions), CF3COOH and HCl (for basic sites of different strength), and O3 (for oxidizable groups). The second part describes a field campaign that has been undertaken in several bus depots in Switzerland, where ambient fine and ultrafine particles were collected on suitable filters and quantitatively investigated using the Knudsen flow reactor. Results point to important differences in the surface reactivity of ambient particles, depending on the sampling site and season. The particle surface appears to be multi-functional, with the simultaneous presence of antagonistic functional groups which do not undergo internal chemical reactions, such as acid-base neutralization. Results also indicate that the surface of ambient particles was characterized by a high density of carbonyl functions (reactivity towards NH2OH probe in the range 0.26-6 formal molecular monolayers) and a low density of acidic sites (reactivity towards N(CH3)3 probe in the range 0.01-0.20 formal molecular monolayer). Kinetic parameters point to fast redox reactions (uptake coefficient ?0>10-3 for O3 probe) and slow acid-base reactions (?0<10-4 for N(CH3)3 probe) on the particle surface. [Authors]

Mutual control of membrane fission and fusion proteins.

Membrane fusion and fission are antagonistic reactions controlled by different proteins. Dynamins promote membrane fission by GTP-driven changes of conformation and polymerization state, while SNAREs fuse membranes by forming complexes between t- and v-SNAREs from apposed vesicles. Here, we describe a role of the dynamin-like GTPase Vps1p in fusion of yeast vacuoles. Vps1p forms polymers that couple several t-SNAREs together. At the onset of fusion, the SNARE-activating ATPase Sec18p/NSF and the t-SNARE depolymerize Vps1p and release it from the membrane. This activity is independent of the SNARE coactivator Sec17p/alpha-SNAP and of the v-SNARE. Vps1p release liberates the t-SNAREs for initiating fusion and at the same time disrupts fission activity. We propose that reciprocal control between fusion and fission components exists, which may prevent futile cycles of fission and fusion.

Formes sévères d'hypersensibilité médicamenteuse retardée [Severe delayed drug hypersensitivity reactions].

Although most delayed drug hypersensitivity reactions are mild and show rapid improvement after drug discontinuation, there are severe systemic and/or cutaneous drug reactions which may be life-threatening. These entities are discussed here, namely DRESS syndrome (Drug Reaction with Eosinophilia and Systemic Symptoms), acute generalized exanthematous pustulosis (AGEP), Stevens-Johnson syndrome and toxic epidermal necrolysis (TEN). Early detection of warning signs and symptoms may help to take appropriate measures precociously.

Detection of millimetric deformation using a terrestrial laser scanner: experiment and application to a rockfall event

Terrestrial laser scanning (TLS) is one of the most promising surveying techniques for rockslope characterization and monitoring. Landslide and rockfall movements can be detected by means of comparison of sequential scans. One of the most pressing challenges of natural hazards is combined temporal and spatial prediction of rockfall. An outdoor experiment was performed to ascertain whether the TLS instrumental error is small enough to enable detection of precursory displacements of millimetric magnitude. This consists of a known displacement of three objects relative to a stable surface. Results show that millimetric changes cannot be detected by the analysis of the unprocessed datasets. Displacement measurement are improved considerably by applying Nearest Neighbour (NN) averaging, which reduces the error (1¿) up to a factor of 6. This technique was applied to displacements prior to the April 2007 rockfall event at Castellfollit de la Roca, Spain. The maximum precursory displacement measured was 45 mm, approximately 2.5 times the standard deviation of the model comparison, hampering the distinction between actual displacement and instrumental error using conventional methodologies. Encouragingly, the precursory displacement was clearly detected by applying the NN averaging method. These results show that millimetric displacements prior to failure can be detected using TLS.

Rockfall monitoring by Terrestrial Laser Scanning ¿ case study of the basaltic rock face at Castellfollit de la Roca (Catalonia, Spain)

This case study deals with a rock face monitoring in urban areas using a Terrestrial Laser Scanner. The pilot study area is an almost vertical, fifty meter high cliff, on top of which the village of Castellfollit de la Roca is located. Rockfall activity is currently causing a retreat of the rock face, which may endanger the houses located at its edge. TLS datasets consist of high density 3-D point clouds acquired from five stations, nine times in a time span of 22 months (from March 2006 to January 2008). The change detection, i.e. rockfalls, was performed through a sequential comparison of datasets. Two types of mass movement were detected in the monitoring period: (a) detachment of single basaltic columns, with magnitudes below 1.5 m3 and (b) detachment of groups of columns, with magnitudes of 1.5 to 150 m3. Furthermore, the historical record revealed (c) the occurrence of slab failures with magnitudes higher than 150 m3. Displacements of a likely slab failure were measured, suggesting an apparent stationary stage. Even failures are clearly episodic, our results, together with the study of the historical record, enabled us to estimate a mean detachment of material from 46 to 91.5 m3 year¿1. The application of TLS considerably improved our understanding of rockfall phenomena in the study area.

Stochastic delay equations with non-negativity constraints driven by fractional Brownian motion

In this note we prove an existence and uniqueness result for the solution of multidimensional stochastic delay differential equations with normal reflection. The equations are driven by a fractional Brownian motion with Hurst parameter H > 1/2. The stochastic integral with respect to the fractional Brownian motion is a pathwise Riemann¿Stieltjes integral.

Stochastic delay differential equations driven by fractional Brownian motion with Hurst parameter H 1/2

We consider the Cauchy problem for a stochastic delay differential equation driven by a fractional Brownian motion with Hurst parameter H>¿. We prove an existence and uniqueness result for this problem, when the coefficients are sufficiently regular. Furthermore, if the diffusion coefficient is bounded away from zero and the coefficients are smooth functions with bounded derivatives of all orders, we prove that the law of the solution admits a smooth density with respect to Lebesgue measure on R.