Hybrid Multiscale Finite Volume method for two-phase flow in porous media

We present a novel hybrid (or multiphysics) algorithm, which couples pore-scale and Darcy descriptions of two-phase flow in porous media. The flow at the pore-scale is described by the Navier?Stokes equations, and the Volume of Fluid (VOF) method is used to model the evolution of the fluid?fluid interface. An extension of the Multiscale Finite Volume (MsFV) method is employed to construct the Darcy-scale problem. First, a set of local interpolators for pressure and velocity is constructed by solving the Navier?Stokes equations; then, a coarse mass-conservation problem is constructed by averaging the pore-scale velocity over the cells of a coarse grid, which act as control volumes; finally, a conservative pore-scale velocity field is reconstructed and used to advect the fluid?fluid interface. The method relies on the localization assumptions used to compute the interpolators (which are quite straightforward extensions of the standard MsFV) and on the postulate that the coarse-scale fluxes are proportional to the coarse-pressure differences. By numerical simulations of two-phase problems, we demonstrate that these assumptions provide hybrid solutions that are in good agreement with reference pore-scale solutions and are able to model the transition from stable to unstable flow regimes. Our hybrid method can naturally take advantage of several adaptive strategies and allows considering pore-scale fluxes only in some regions, while Darcy fluxes are used in the rest of the domain. Moreover, since the method relies on the assumption that the relationship between coarse-scale fluxes and pressure differences is local, it can be used as a numerical tool to investigate the limits of validity of Darcy's law and to understand the link between pore-scale quantities and their corresponding Darcy-scale variables.

Les enjeux spatiaux et sociaux d'une dynamique interethnique en transition dans les quartiers de centre-ville. La construction des diverses formes de représentation dans la manière de définir un quartier de Nice.

Résumé Cette thèse de doctorat se fixe comme objectif général de définir la spécificité d'un espace urbain situé en centre-ville à travers sa dynamique sociale et interethnique. Il s'agit de relever les différentes formes et manifestations ethniques, mais également de comprendre la manière dont est vécue cette socialité, tout en tenant compte du contexte urbain dans lequel elles se situent. La particularité de ce travail est le fait de proposer une démarche alliant deux disciplines à savoir la géographie et la sociologie, dans le but de proposer une recherche transdisciplinaire, avec comme conséquence la mise en oeuvre d'une démarche commune nourrie des deux disciplines. Les quartiers ethniques situés en centre-ville ne se présentent pas comme des espaces figés, puisque des changements fréquents sont observés aussi bien au niveau de leur population que de l'affectation des locaux commerciaux et des logements. Ils semblent davantage se caractériser comme des « espaces à l'avenir indéterminé » dans le sens où leur avenir aussi bien au niveau de leur structure urbaine que de leurs rapports sociaux ne peut être prévisible et qu'ils sont par ailleurs susceptibles d'être soumis à de nombreux changements, d'où l'idée également d' « espaces en mutations ». Pour aborder ces thématiques de recherche, nous avons opté pour une étude fondée sur le principe de la représentation et de l'une des formes essentielles à travers lesquelles elles se produisent et se manifestent, à savoir les images. Cette recherche privilégie ainsi une analyse dynamique du phénomène interethnique, l'inscrivant dans le courant constructiviste qui s'intéresse notamment à l'étude des pratiques des acteurs sociaux et leurs représentations, dans le sens qu'elles sont le fruit d'une construction sociale. Dans ce type de démarche analytique, il s'agit de prendre en considération à la fois l'analyse du chercheur et la réalité objective et subjective du vécu des acteurs urbains qui regroupent les riverains, les commerçants, les membres associatifs, les politiques et les médias. C'est justement l'association de ces deux approches qui permet de définir un quartier urbain. Abstract "The spatial and social stakes of an interethnic dynamics in transition in the downtown area. The construction of the various forms of representations in order to define an urban quarter in Nice" The general objective of this PHD is to define the specificity of an urban space located in the downtown area with social and interethnic dynamics. Specifically the diverse ethnic form and sign to include the understanding mariner where society has lived in a specific urban context. The particularity of this work is to advance an analytic method using two disciplines: geography and sociology, in order to bring up a transdisciplinarian research with a common method. The ethnic quarters in the inner-city area are not presented as fixed spaces, since frequent changes were observed in the population level and at the allotment level of the commercial premises and the lodgings. Defined as "undetermined future spaces" as their future in the urban structure or the social relations can not be foresee and are supposed to be subjective to many changes. Therefore speaking of "mutations spaces". The option of this type of thematic research is the study of the principle of representation and essential form which shows and produces the images. Research promotes a dynamic analyses of the interethnic phenomena, inside the constructivist which is focused the practices of social actors, their own representations, and products of a social construction. The production of analytical approach is necessary in order to consider the analyses of the research worker and the objective and subjective reality of the life of the urban actors. It gathers the residents, tradesmen, members of associations, politicians and the press. The association of these types of approach define a urban quarter.

Using a single tablet daily to treat latent tuberculosis infection in Brazil: bioequivalence of two different isoniazid formulations (300 mg and 100 mg) demonstrated by a sensitive and rapid high-performance liquid chromatography-tandem mass spectrometry method in a randomised, crossover study

The recommended treatment for latent tuberculosis (TB) infection in adults is a daily dose of isoniazid (INH) 300 mg for six months. In Brazil, INH was formulated as 100 mg tablets. The treatment duration and the high pill burden compromised patient adherence to the treatment. The Brazilian National Programme for Tuberculosis requested a new 300 mg INH formulation. The aim of our study was to compare the bioavailability of the new INH 300 mg formulation and three 100 mg tablets of the reference formulation. We conducted a randomised, single dose, open label, two-phase crossover bioequivalence study in 28 healthy human volunteers. The 90% confidence interval for the INH maximum concentration of drug observed in plasma and area under the plasma concentration vs. time curve from time zero to the last measurable concentration “time t” was 89.61-115.92 and 94.82-119.44, respectively. The main limitation of our study was that neither adherence nor the safety profile of multiple doses was evaluated. To determine the level of INH in human plasma, we developed and validated a sensitive, simple and rapid high-performance liquid chromatography-tandem mass spectrometry method. Our results showed that the new formulation was bioequivalent to the 100 mg reference product. This finding supports the use of a single 300 mg tablet daily strategy to treat latent TB. This new formulation may increase patients’ adherence to the treatment and quality of life.

Spatial relationship between the atmospheric circulation and the precipitation measured in the western Swiss Alps by means of the analogue method

An adaptation technique based on the synoptic atmospheric circulation to forecast local precipitation, namely the analogue method, has been implemented for the western Swiss Alps. During the calibration procedure, relevance maps were established for the geopotential height data. These maps highlight the locations were the synoptic circulation was found of interest for the precipitation forecasting at two rain gauge stations (Binn and Les Marécottes) that are located both in the alpine Rhône catchment, at a distance of about 100 km from each other. These two stations are sensitive to different atmospheric circulations. We have observed that the most relevant data for the analogue method can be found where specific atmospheric circulation patterns appear concomitantly with heavy precipitation events. Those skilled regions are coherent with the atmospheric flows illustrated, for example, by means of the back trajectories of air masses. Indeed, the circulation recurrently diverges from the climatology during days with strong precipitation on the southern part of the alpine Rhône catchment. We have found that for over 152 days with precipitation amount above 50 mm at the Binn station, only 3 did not show a trajectory of a southerly flow, meaning that such a circulation was present for 98% of the events. Time evolution of the relevance maps confirms that the atmospheric circulation variables have significantly better forecasting skills close to the precipitation period, and that it seems pointless for the analogue method to consider circulation information days before a precipitation event as a primary predictor. Even though the occurrence of some critical circulation patterns leading to heavy precipitation events can be detected by precursors at remote locations and 1 week ahead (Grazzini, 2007; Martius et al., 2008), time extrapolation by the analogue method seems to be rather poor. This would suggest, in accordance with previous studies (Obled et al., 2002; Bontron and Obled, 2005), that time extrapolation should be done by the Global Circulation Model, which can process atmospheric variables that can be used by the adaptation method.

Estimation of electronic coupling in π -stacked donor-bridge-acceptor systems: correction of the two-state model

Comparison of donor-acceptor electronic couplings calculated within two-state and three-state models suggests that the two-state treatment can provide unreliable estimates of Vda because of neglecting the multistate effects. We show that in most cases accurate values of the electronic coupling in a π stack, where donor and acceptor are separated by a bridging unit, can be obtained as Ṽ da = (E2 - E1) μ12 Rda + (2 E3 - E1 - E2) 2 μ13 μ23 Rda2, where E1, E2, and E3 are adiabatic energies of the ground, charge-transfer, and bridge states, respectively, μij is the transition dipole moments between the states i and j, and Rda is the distance between the planes of donor and acceptor. In this expression based on the generalized Mulliken-Hush approach, the first term corresponds to the coupling derived within a two-state model, whereas the second term is the superexchange correction accounting for the bridge effect. The formula is extended to bridges consisting of several subunits. The influence of the donor-acceptor energy mismatch on the excess charge distribution, adiabatic dipole and transition moments, and electronic couplings is examined. A diagnostic is developed to determine whether the two-state approach can be applied. Based on numerical results, we showed that the superexchange correction considerably improves estimates of the donor-acceptor coupling derived within a two-state approach. In most cases when the two-state scheme fails, the formula gives reliable results which are in good agreement (within 5%) with the data of the three-state generalized Mulliken-Hush model

Estimates of electronic coupling for excess electron transfer in DNA

Electronic coupling Vda is one of the key parameters that determine the rate of charge transfer through DNA. While there have been several computational studies of Vda for hole transfer, estimates of electronic couplings for excess electron transfer (ET) in DNA remain unavailable. In the paper, an efficient strategy is established for calculating the ET matrix elements between base pairs in a π stack. Two approaches are considered. First, we employ the diabatic-state (DS) method in which donor and acceptor are represented with radical anions of the canonical base pairs adenine-thymine (AT) and guanine-cytosine (GC). In this approach, similar values of Vda are obtained with the standard 6-31 G* and extended 6-31++ G* basis sets. Second, the electronic couplings are derived from lowest unoccupied molecular orbitals (LUMOs) of neutral systems by using the generalized Mulliken-Hush or fragment charge methods. Because the radical-anion states of AT and GC are well reproduced by LUMOs of the neutral base pairs calculated without diffuse functions, the estimated values of Vda are in good agreement with the couplings obtained for radical-anion states using the DS method. However, when the calculation of a neutral stack is carried out with diffuse functions, LUMOs of the system exhibit the dipole-bound character and cannot be used for estimating electronic couplings. Our calculations suggest that the ET matrix elements Vda for models containing intrastrand thymine and cytosine bases are essentially larger than the couplings in complexes with interstrand pyrimidine bases. The matrix elements for excess electron transfer are found to be considerably smaller than the corresponding values for hole transfer and to be very responsive to structural changes in a DNA stack

Two complementary molecular energy decomposition schemes: the Mayer and Ziegler-Rauk methods in comparison

In the present paper we discuss and compare two different energy decomposition schemes: Mayer's Hartree-Fock energy decomposition into diatomic and monoatomic contributions [Chem. Phys. Lett. 382, 265 (2003)], and the Ziegler-Rauk dissociation energy decomposition [Inorg. Chem. 18, 1558 (1979)]. The Ziegler-Rauk scheme is based on a separation of a molecule into fragments, while Mayer's scheme can be used in the cases where a fragmentation of the system in clearly separable parts is not possible. In the Mayer scheme, the density of a free atom is deformed to give the one-atom Mulliken density that subsequently interacts to give rise to the diatomic interaction energy. We give a detailed analysis of the diatomic energy contributions in the Mayer scheme and a close look onto the one-atom Mulliken densities. The Mulliken density ρA has a single large maximum around the nuclear position of the atom A, but exhibits slightly negative values in the vicinity of neighboring atoms. The main connecting point between both analysis schemes is the electrostatic energy. Both decomposition schemes utilize the same electrostatic energy expression, but differ in how fragment densities are defined. In the Mayer scheme, the electrostatic component originates from the interaction of the Mulliken densities, while in the Ziegler-Rauk scheme, the undisturbed fragment densities interact. The values of the electrostatic energy resulting from the two schemes differ significantly but typically have the same order of magnitude. Both methods are useful and complementary since Mayer's decomposition focuses on the energy of the finally formed molecule, whereas the Ziegler-Rauk scheme describes the bond formation starting from undeformed fragment densities

One- and two-center physical space partitioning of the energy in the density functional theory

A conceptually new approach is introduced for the decomposition of the molecular energy calculated at the density functional theory level of theory into sum of one- and two-atomic energy components, and is realized in the "fuzzy atoms" framework. (Fuzzy atoms mean that the three-dimensional physical space is divided into atomic regions having no sharp boundaries but exhibiting a continuous transition from one to another.) The new scheme uses the new concept of "bond order density" to calculate the diatomic exchange energy components and gives them unexpectedly close to the values calculated by the exact (Hartree-Fock) exchange for the same Kohn-Sham orbitals

Influence of configurational atomic order on the relative stability of bcc and close-packed structures in Cu-based alloys

We present a study of the influence of atomic order on the relative stability of the bcc and the 18R martensitic structures in a Cu2.96Al0.92Be0.12 crystal. Calorimetric measurements have shown that disorder increases the stability of the 18R phase, contrary to what happens in Cu-Zn-Al alloys for which it is the bcc phase that is stabilized by disordering the system. This different behavior has been explained in terms of a model recently reported. We have also proved that the entropy change at the martensitic transition is independent of the state of atomic order of the crystal, as predicted theoretically. Our results suggest that differences in the vibrational spectrum of the crystal due to different states of atomic order must be equal in the bcc and in the close-packed phases.

Vacancy-driven ordering in a two-dimensional binary alloy

Domain growth in a system with nonconserved order parameter is studied. We simulate the usual Ising model for binary alloys with concentration 0.5 on a two-dimensional square lattice by Monte Carlo techniques. Measurements of the energy, jump-acceptance ratio, and order parameters are performed. Dynamics based on the diffusion of a single vacancy in the system gives a growth law faster than the usual Allen-Cahn law. Allowing vacancy jumps to next-nearest-neighbor sites is essential to prevent vacancy trapping in the ordered regions. By measuring local order parameters we show that the vacancy prefers to be in the disordered regions (domain boundaries). This naturally concentrates the atomic jumps in the domain boundaries, accelerating the growth compared with the usual exchange mechanism that causes jumps to be homogeneously distributed on the lattice.

Pairing in a two-component ultracold Fermi gas: Phases with broken-space symmetries

We explore the phase diagram of a two-component ultracold atomic Fermi gas interacting with zero-range forces in the limit of weak coupling. We focus on the dependence of the pairing gap and the free energy on the variations in the number densities of the two species while the total density of the system is held fixed. As the density asymmetry is increased, the system exhibits a transition from a homogenous Bardeen-Cooper-Schrieffer (BCS) phase to phases with spontaneously broken global space symmetries. One such realization is the deformed Fermi surface superfluidity (DFS) which exploits the possibility of deforming the Fermi surfaces of the species into ellipsoidal form at zero total momentum of Cooper pairs. The critical asymmetries at which the transition from DFS to the unpaired state occurs are larger than those for the BCS phase. In this precritical region the DFS phase lowers the pairing energy of the asymmetric BCS state. We compare quantitatively the DFS phase to another realization of superconducting phases with broken translational symmetry: the single-plane-wave Larkin-Ovchinnikov-Fulde-Ferrell phase, which is characterized by a nonvanishing center-of-mass momentum of the Cooper pairs. The possibility of the detection of the DFS phase in the time-of-flight experiments is discussed and quantified for the case of 6Li atoms trapped in two different hyperfine states.

MERCURY QUANTIFICATION IN SOILS USING THERMAL DESORPTION AND ATOMIC ABSORPTION SPECTROMETRY: PROPOSAL FOR AN ALTERNATIVE METHOD OF ANALYSIS

Despite the considerable environmental importance of mercury (Hg), given its high toxicity and ability to contaminate large areas via atmospheric deposition, little is known about its activity in soils, especially tropical soils, in comparison with other heavy metals. This lack of information about Hg arises because analytical methods for determination of Hg are more laborious and expensive compared to methods for other heavy metals. The situation is even more precarious regarding speciation of Hg in soils since sequential extraction methods are also inefficient for this metal. The aim of this paper is to present a technique of thermal desorption associated with atomic absorption spectrometry, TDAAS, as an efficient tool for quantitative determination of Hg in soils. The method consists of the release of Hg by heating, followed by its quantification by atomic absorption spectrometry. It was developed by constructing calibration curves in different soil samples based on increasing volumes of standard Hg2+ solutions. Performance, accuracy, precision, and quantification and detection limit parameters were evaluated. No matrix interference was detected. Certified reference samples and comparison with a Direct Mercury Analyzer, DMA (another highly recognized technique), were used in validation of the method, which proved to be accurate and precise.

Experimental Method to Determine Some Physical Properties in Physics Classes

ABSTRACT Particle density, gravimetric and volumetric water contents and porosity are important basic concepts to characterize porous systems such as soils. This paper presents a proposal of an experimental method to measure these physical properties, applicable in experimental physics classes, in porous media samples consisting of spheres with the same diameter (monodisperse medium) and with different diameters (polydisperse medium). Soil samples are not used given the difficulty of working with this porous medium in laboratories dedicated to teaching basic experimental physics. The paper describes the method to be followed and results of two case studies, one in monodisperse medium and the other in polydisperse medium. The particle density results were very close to theoretical values for lead spheres, whose relative deviation (RD) was -2.9 % and +0.1 % RD for the iron spheres. The RD of porosity was also low: -3.6 % for lead spheres and -1.2 % for iron spheres, in the comparison of procedures – using particle and porous medium densities and saturated volumetric water content – and monodisperse and polydisperse media.

Impaired chemical coupling of cerebral blood flow is compatible with intact neurological outcome in neonates with perinatal risk factors.

Early detection of pathophysiological factors associated with permanent brain damage is a major issue in neonatal medicine. The aim of our study was to evaluate the significance of the CO2 reactivity of cerebral blood flow (CBF) in neonates with perinatal risk factors. Fourteen ventilated neonates with perinatal risk factors (pathological cardiotocogramm, low cord pH, postpartal encephalopathy) were enrolled into this prospective study. The study was performed 18-123 h after birth. CBF was measured using the noninvasive intravenous 133Xe method. Two measurements were taken with a minimal PaCO2-difference of 5 mm Hg. From the two CBF values the CO2 reactivity was calculated. Outcome was evaluated 1 year after birth. The CBF values at a lower PaCO2 ranged from 6.6 to 115. 2 ml/100 g brain issue/min (median = 18.2) and at a higher PaCO2 level from 7.1 to 125.7 ml/100 g brain tissue/min (median = 18.75). The calculated CO2 reactivity ranged from -9.6 to 6.6% (median 1.1%) change in CBF/mm Hg change in PaCO2. CO2 reactivity correlated with lowest pH (r2 = 0.35, p = 0.02). Two infants died, one of neonatal sepsis, the other of heart failure. Neurological outcome at the age of 1 year was normal in 11 patients, 1 had severe cerebral palsy. From the 12 surviving patients the patient with severe neurological deficit showed the highest CBF values (125.7 ml/100 g/min). Impaired chemical coupling of cerebral blood flow is compatible with intact neurological outcome in neonates with perinatal risk factors. CO2 reactivity in these newborns correlates with the lowest pH and may reflect the severity of perinatal asphyxia.