Ballistic and diffusive corrections to front propagation in the presence of multiplicative noise

We study the dynamics of reaction-diffusion fronts under the influence of multiplicative noise. An approximate theoretical scheme is introduced to compute the velocity of the front and its diffusive wandering due to the presence of noise. The theoretical approach is based on a multiple scale analysis rather than on a small noise expansion and is confirmed with numerical simulations for a wide range of the noise intensity. We report on the possibility of noise sustained solutions with a continuum of possible velocities, in situations where only a single velocity is allowed without noise.

Front dynamics in turbulent media

A study of a stable front propagating in a turbulent medium is presented. The front is generated through a reaction-diffusion equation, and the turbulent medium is statistically modeled using a Langevin equation. Numerical simulations indicate the presence of two different dynamical regimes. These regimes appear when the turbulent flow either wrinkles a still rather sharp propagating interfase or broadens it. Specific dependences of the propagating velocities on stirring intensities appropriate to each case are found and fitted when possible according to theoretically predicted laws. Different turbulent spectra are considered.

Physicochemical and biological characterization of single-walled and double-walled carbon nanotubes in biological media

To study the toxicity of nanoparticles under relevant conditions, it is important to reproducibly disperse nanoparticles in biological media in in vitro and in vivo studies. Here, single-walled nanotubes (SWNTs) and double-walled nanotubes (DWNTs) were physicochemically and biologically characterized when dispersed in phosphate-buffered saline (PBS) and bovine serum albumin (BSA). BSA-SWNT/DWNT interaction resulted in a reduction of aggregation and an increase in particle stabilization. Based on the protein sequence coverage and protein binding results, DWNTs exhibited higher protein binding than SWNTs. SWNT and DWNT suspensions in the presence of BSA increased interleukin-6 (IL-6) levels and reduced tumor necrosis factor-alpha (TNF-α) levels in A549 cells as compared to corresponding samples in the absence of BSA. We next determined the effects of SWNTs and DWNTs on pulmonary protein modification using bronchoalveolar lavage fluid (BALF) as a surrogate collected form BALB/c mice. The BALF proteins bound to SWNTs (13 proteins) and DWNTs (11 proteins), suggesting that these proteins were associated with blood coagulation pathways. Lastly, we demonstrated the importance of physicochemical and biological alterations of SWNTs and DWNTs when dispersed in biological media, since protein binding may result in the misinterpretation of in vitro results and the activation of protein-regulated biological responses.

Propagation in a thermal graviton background

It is well known that radiative corrections evaluated in nontrivial backgrounds lead to effective dispersion relations which are not Lorentz invariant. Since gravitational interactions increase with energy, gravity-induced radiative corrections could be relevant for the trans-Planckian problem. As a first step to explore this possibility, we compute the one-loop radiative corrections to the self-energy of a scalar particle propagating in a thermal bath of gravitons in Minkowski spacetime. We obtain terms which originate from the thermal bath and which indeed break the Lorentz invariance that possessed the propagator in the vacuum. Rather unexpectedly, however, the terms which break Lorentz invariance vanish in the high three-momentum limit. We also found that the imaginary part, which gives the rate of approach to thermal equilibrium, vanishes at one loop.

Multiscale descriptions of density-driven flow instabilities in porous media

Les instabilités engendrées par des gradients de densité interviennent dans une variété d'écoulements. Un exemple est celui de la séquestration géologique du dioxyde de carbone en milieux poreux. Ce gaz est injecté à haute pression dans des aquifères salines et profondes. La différence de densité entre la saumure saturée en CO2 dissous et la saumure environnante induit des courants favorables qui le transportent vers les couches géologiques profondes. Les gradients de densité peuvent aussi être la cause du transport indésirable de matières toxiques, ce qui peut éventuellement conduire à la pollution des sols et des eaux. La gamme d'échelles intervenant dans ce type de phénomènes est très large. Elle s'étend de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères à laquelle interviennent les phénomènes à temps long. Une reproduction fiable de la physique par la simulation numérique demeure donc un défi en raison du caractère multi-échelles aussi bien au niveau spatial et temporel de ces phénomènes. Il requiert donc le développement d'algorithmes performants et l'utilisation d'outils de calculs modernes. En conjugaison avec les méthodes de résolution itératives, les méthodes multi-échelles permettent de résoudre les grands systèmes d'équations algébriques de manière efficace. Ces méthodes ont été introduites comme méthodes d'upscaling et de downscaling pour la simulation d'écoulements en milieux poreux afin de traiter de fortes hétérogénéités du champ de perméabilité. Le principe repose sur l'utilisation parallèle de deux maillages, le premier est choisi en fonction de la résolution du champ de perméabilité (grille fine), alors que le second (grille grossière) est utilisé pour approximer le problème fin à moindre coût. La qualité de la solution multi-échelles peut être améliorée de manière itérative pour empêcher des erreurs trop importantes si le champ de perméabilité est complexe. Les méthodes adaptatives qui restreignent les procédures de mise à jour aux régions à forts gradients permettent de limiter les coûts de calculs additionnels. Dans le cas d'instabilités induites par des gradients de densité, l'échelle des phénomènes varie au cours du temps. En conséquence, des méthodes multi-échelles adaptatives sont requises pour tenir compte de cette dynamique. L'objectif de cette thèse est de développer des algorithmes multi-échelles adaptatifs et efficaces pour la simulation des instabilités induites par des gradients de densité. Pour cela, nous nous basons sur la méthode des volumes finis multi-échelles (MsFV) qui offre l'avantage de résoudre les phénomènes de transport tout en conservant la masse de manière exacte. Dans la première partie, nous pouvons démontrer que les approximations de la méthode MsFV engendrent des phénomènes de digitation non-physiques dont la suppression requiert des opérations de correction itératives. Les coûts de calculs additionnels de ces opérations peuvent toutefois être compensés par des méthodes adaptatives. Nous proposons aussi l'utilisation de la méthode MsFV comme méthode de downscaling: la grille grossière étant utilisée dans les zones où l'écoulement est relativement homogène alors que la grille plus fine est utilisée pour résoudre les forts gradients. Dans la seconde partie, la méthode multi-échelle est étendue à un nombre arbitraire de niveaux. Nous prouvons que la méthode généralisée est performante pour la résolution de grands systèmes d'équations algébriques. Dans la dernière partie, nous focalisons notre étude sur les échelles qui déterminent l'évolution des instabilités engendrées par des gradients de densité. L'identification de la structure locale ainsi que globale de l'écoulement permet de procéder à un upscaling des instabilités à temps long alors que les structures à petite échelle sont conservées lors du déclenchement de l'instabilité. Les résultats présentés dans ce travail permettent d'étendre les connaissances des méthodes MsFV et offrent des formulations multi-échelles efficaces pour la simulation des instabilités engendrées par des gradients de densité. - Density-driven instabilities in porous media are of interest for a wide range of applications, for instance, for geological sequestration of CO2, during which CO2 is injected at high pressure into deep saline aquifers. Due to the density difference between the C02-saturated brine and the surrounding brine, a downward migration of CO2 into deeper regions, where the risk of leakage is reduced, takes place. Similarly, undesired spontaneous mobilization of potentially hazardous substances that might endanger groundwater quality can be triggered by density differences. Over the last years, these effects have been investigated with the help of numerical groundwater models. Major challenges in simulating density-driven instabilities arise from the different scales of interest involved, i.e., the scale at which instabilities are triggered and the aquifer scale over which long-term processes take place. An accurate numerical reproduction is possible, only if the finest scale is captured. For large aquifers, this leads to problems with a large number of unknowns. Advanced numerical methods are required to efficiently solve these problems with today's available computational resources. Beside efficient iterative solvers, multiscale methods are available to solve large numerical systems. Originally, multiscale methods have been developed as upscaling-downscaling techniques to resolve strong permeability contrasts. In this case, two static grids are used: one is chosen with respect to the resolution of the permeability field (fine grid); the other (coarse grid) is used to approximate the fine-scale problem at low computational costs. The quality of the multiscale solution can be iteratively improved to avoid large errors in case of complex permeability structures. Adaptive formulations, which restrict the iterative update to domains with large gradients, enable limiting the additional computational costs of the iterations. In case of density-driven instabilities, additional spatial scales appear which change with time. Flexible adaptive methods are required to account for these emerging dynamic scales. The objective of this work is to develop an adaptive multiscale formulation for the efficient and accurate simulation of density-driven instabilities. We consider the Multiscale Finite-Volume (MsFV) method, which is well suited for simulations including the solution of transport problems as it guarantees a conservative velocity field. In the first part of this thesis, we investigate the applicability of the standard MsFV method to density- driven flow problems. We demonstrate that approximations in MsFV may trigger unphysical fingers and iterative corrections are necessary. Adaptive formulations (e.g., limiting a refined solution to domains with large concentration gradients where fingers form) can be used to balance the extra costs. We also propose to use the MsFV method as downscaling technique: the coarse discretization is used in areas without significant change in the flow field whereas the problem is refined in the zones of interest. This enables accounting for the dynamic change in scales of density-driven instabilities. In the second part of the thesis the MsFV algorithm, which originally employs one coarse level, is extended to an arbitrary number of coarse levels. We prove that this keeps the MsFV method efficient for problems with a large number of unknowns. In the last part of this thesis, we focus on the scales that control the evolution of density fingers. The identification of local and global flow patterns allows a coarse description at late times while conserving fine-scale details during onset stage. Results presented in this work advance the understanding of the Multiscale Finite-Volume method and offer efficient dynamic multiscale formulations to simulate density-driven instabilities. - Les nappes phréatiques caractérisées par des structures poreuses et des fractures très perméables représentent un intérêt particulier pour les hydrogéologues et ingénieurs environnementaux. Dans ces milieux, une large variété d'écoulements peut être observée. Les plus communs sont le transport de contaminants par les eaux souterraines, le transport réactif ou l'écoulement simultané de plusieurs phases non miscibles, comme le pétrole et l'eau. L'échelle qui caractérise ces écoulements est définie par l'interaction de l'hétérogénéité géologique et des processus physiques. Un fluide au repos dans l'espace interstitiel d'un milieu poreux peut être déstabilisé par des gradients de densité. Ils peuvent être induits par des changements locaux de température ou par dissolution d'un composé chimique. Les instabilités engendrées par des gradients de densité revêtent un intérêt particulier puisque qu'elles peuvent éventuellement compromettre la qualité des eaux. Un exemple frappant est la salinisation de l'eau douce dans les nappes phréatiques par pénétration d'eau salée plus dense dans les régions profondes. Dans le cas des écoulements gouvernés par les gradients de densité, les échelles caractéristiques de l'écoulement s'étendent de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères sur laquelle interviennent les phénomènes à temps long. Etant donné que les investigations in-situ sont pratiquement impossibles, les modèles numériques sont utilisés pour prédire et évaluer les risques liés aux instabilités engendrées par les gradients de densité. Une description correcte de ces phénomènes repose sur la description de toutes les échelles de l'écoulement dont la gamme peut s'étendre sur huit à dix ordres de grandeur dans le cas de grands aquifères. Il en résulte des problèmes numériques de grande taille qui sont très couteux à résoudre. Des schémas numériques sophistiqués sont donc nécessaires pour effectuer des simulations précises d'instabilités hydro-dynamiques à grande échelle. Dans ce travail, nous présentons différentes méthodes numériques qui permettent de simuler efficacement et avec précision les instabilités dues aux gradients de densité. Ces nouvelles méthodes sont basées sur les volumes finis multi-échelles. L'idée est de projeter le problème original à une échelle plus grande où il est moins coûteux à résoudre puis de relever la solution grossière vers l'échelle de départ. Cette technique est particulièrement adaptée pour résoudre des problèmes où une large gamme d'échelle intervient et évolue de manière spatio-temporelle. Ceci permet de réduire les coûts de calculs en limitant la description détaillée du problème aux régions qui contiennent un front de concentration mobile. Les aboutissements sont illustrés par la simulation de phénomènes tels que l'intrusion d'eau salée ou la séquestration de dioxyde de carbone.

Generalized percolation in random directed networks

We develop a general theory for percolation in directed random networks with arbitrary two-point correlations and bidirectional edgesthat is, edges pointing in both directions simultaneously. These two ingredients alter the previously known scenario and open new views and perspectives on percolation phenomena. Equations for the percolation threshold and the sizes of the giant components are derived in the most general case. We also present simulation results for a particular example of uncorrelated network with bidirectional edges confirming the theoretical predictions.

Tuning clustering in random networks with arbitrary degree distributions

We present a generator of random networks where both the degree-dependent clustering coefficient and the degree distribution are tunable. Following the same philosophy as in the configuration model, the degree distribution and the clustering coefficient for each class of nodes of degree k are fixed ad hoc and a priori. The algorithm generates corresponding topologies by applying first a closure of triangles and second the classical closure of remaining free stubs. The procedure unveils an universal relation among clustering and degree-degree correlations for all networks, where the level of assortativity establishes an upper limit to the level of clustering. Maximum assortativity ensures no restriction on the decay of the clustering coefficient whereas disassortativity sets a stronger constraint on its behavior. Correlation measures in real networks are seen to observe this structural bound.

Quantum relaxation in random magnets

We present a comprehensive study of the low-temperature magnetic relaxation in random magnets. The first part of the paper contains theoretical analysis of the expected features of the relaxation, based upon current theories of quantum tunneling of magnetization. Models of tunneling, dissipation, the crossover from the thermal to the quantum regime, and the effect of barrier distribution on the relaxation rate are discussed. It is argued that relaxation-type experiments are ideally suited for the observation of magnetic tunneling, since they automatically provide the condition of very low barriers. The second part of the paper contains experimental results on transition-metal¿rare-earth amorphous magnets. Structural and magnetic characterization of materials is presented. The temperature and field dependence of the magnetic relaxation is studied. Our key observation is a nonthermal character of the relaxation below a few kelvin. The observed features are in agreement with theoretical suggestions on quantum tunneling of magnetization.

Metastability and front propagation in the first-order optical Fréedericksz transition

A general dynamical model for the first-order optical Fréedericksz transition incorporating spatial transverse inhomogeneities and hydrodynamic effects is discussed in the framework of a time-dependent Ginzburg-Landau model. The motion of an interface between two coexisting states with different director orientations is considered. A uniformly translating front solution of the dynamical equations for the motion of that interface is described.

Enhanced pulse propagation in non-linear arrays of oscillators

The propagation of a pulse in a nonlinear array of oscillators is influenced by the nature of the array and by its coupling to a thermal environment. For example, in some arrays a pulse can be speeded up while in others a pulse can be slowed down by raising the temperature. We begin by showing that an energy pulse (one dimension) or energy front (two dimensions) travels more rapidly and remains more localized over greater distances in an isolated array (microcanonical) of hard springs than in a harmonic array or in a soft-springed array. Increasing the pulse amplitude causes it to speed up in a hard chain, leaves the pulse speed unchanged in a harmonic system, and slows down the pulse in a soft chain. Connection of each site to a thermal environment (canonical) affects these results very differently in each type of array. In a hard chain the dissipative forces slow down the pulse while raising the temperature speeds it up. In a soft chain the opposite occurs: the dissipative forces actually speed up the pulse, while raising the temperature slows it down. In a harmonic chain neither dissipation nor temperature changes affect the pulse speed. These and other results are explained on the basis of the frequency vs energy relations in the various arrays

Front dynamics in turbulent media

A study of a stable front propagating in a turbulent medium is presented. The front is generated through a reaction-diffusion equation, and the turbulent medium is statistically modeled using a Langevin equation. Numerical simulations indicate the presence of two different dynamical regimes. These regimes appear when the turbulent flow either wrinkles a still rather sharp propagating interfase or broadens it. Specific dependences of the propagating velocities on stirring intensities appropriate to each case are found and fitted when possible according to theoretically predicted laws. Different turbulent spectra are considered.

Changes in Aortic Pulse Wave Velocity in Hypertensive Postmenopausal Women: Comparison Between a Calcium Channel Blocker vs Angiotensin Receptor Blocker Regimen.

J Clin Hypertens (Greenwich). 2012;14:773-778. ©2012 Wiley Periodicals, Inc. Postmenopausal women are at greater risk for hypertension-related cardiovascular disease. Antihypertensive therapy may help alleviate arterial stiffness that represents a potential modifiable risk factor of hypertension. This randomized controlled study investigated the difference between an angiotensin receptor blocker and a calcium channel blocker in reducing arterial stiffness. Overall, 125 postmenopausal hypertensive women (age, 61.4±6 years; systolic blood pressure/diastolic blood pressure [SBP/DBP], 158±11/92±9 mm Hg) were randomized to valsartan 320 mg±hydrochlorothiazide (HCTZ) (n=63) or amlodipine 10 mg±HCTZ (n=62). The primary outcome was carotid-to-femoral pulse wave velocity (PWV) changes after 38 weeks of treatment. Both treatments lowered peripheral blood pressure (BP) (-22.9/-10.9 mm Hg for valsartan and -25.2/-11.7 mm Hg for amlodipine, P=not significant) and central BP (-15.7/-7.6 mm Hg for valsartan and -19.2/-10.3 mm Hg for amlodipine, P<.05 for central DBP). Both treatments similarly reduced the carotid-femoral PWV (-1.9 vs -1.7 m/s; P=not significant). Amlodipine was associated with a higher incidence of peripheral edema compared with the valsartan group (77% vs 14%, P<.001). BP lowering in postmenopausal women led to a reduction in arterial stiffness as assessed by PWV measurement. Both regimens reduced PWV to a similar degree after 38 weeks of treatment despite differences in central BP lowering, suggesting that the effect of valsartan on PWV is mediated through nonhemodynamic effects.

Changes in Aortic Pulse Wave Velocity in Hypertension Post-Menopausal Women: Comparison Between Calcium Channel Blocker (CCB) Vs Angiotensin Receptor Blocker (ARB) Regimen

Objective: To compare effects of a non-renin-angiotensin system (RAS) blocker, using a CCB, or a RAS blocker, using an ARB regimen on the arterial stiffness reduction in postmenopausal hypertensive women. Methods: In this prospective study, a total of 125 hypertensive women (age: 61.4_6 yrs; 98% Caucasian; BW: 71.9_14 kg; BMI: 27.3_5 kg/m2; SBP/ DBP: 158_11/92_9 mmHg) were randomized between ARB (valsartan 320mg_HCTZ) and CCB (amlodipine 10mg _ HCTZ). The primary outcome was carotid-femoral pulse wave velocity (PWV) changes after 38 weeks of treatment. Results: There were no significant differences in baseline demographic data between the two groups. Both treatments effectively lowered BP at the end of the study with similar (p>0.05) reductions in the valsartan (_22.9/_10.9 mmHg) and amlodipine based (_25.2/_11.7 mmHg) treatment groups. Despite a lower (p<0.05 for DBP) central SBP/DBP in the CCB group (_19.2/_10.3 mmHg) compared to the valsartan group (_15.7/_7.6 mmHg) at week 38, a similar reduction in carotid-femoral PWV (_1.7 vs _1.9 m/sec; p>0.05) was observed between both groups. The numerically larger BP reduction observed in the CCB group was associated with a much higher incidence of peripheral edema (77% vs 14%) than the valsartan group. Conclusion: In summary, BP lowering in postmenopausal women led to a reduction in arterial stiffness assessed by PWV measurement. Both regimens reduced PWV at 38 weeks of treatment to a similar degree, despite differences in BP lowering suggesting that the effect of RAS blockade to influence PWV may partly be independent of BP.

A framework for hybrid simulations of two-phase flow in porous media

Rough a global coarse problem. Although these techniques are usually employed for problems in which the fine-scale processes are described by Darcy's law, they can also be applied to pore-scale simulations and used as a mathematical framework for hybrid methods that couples a Darcy and pore scales. In this work, we consider a pore-scale description of fine-scale processes. The Navier-Stokes equations are numerically solved in the pore geometry to compute the velocity field and obtain generalized permeabilities. In the case of two-phase flow, the dynamics of the phase interface is described by the volume of fluid method with the continuum surface force model. The MsFV method is employed to construct an algorithm that couples a Darcy macro-scale description with a pore-scale description at the fine scale. The hybrid simulations results presented are in good agreement with the fine-scale reference solutions. As the reconstruction of the fine-scale details can be done adaptively, the presented method offers a flexible framework for hybrid modeling.