Interaccion among systems of finite size in predictive relativistic mechanics I. General framework

We derive a Hamiltonian formulation for the three-dimensional formalism of predictive relativistic mechanics. This Hamiltonian structure is used to derive a set of dynamical equations describing the interaction among systems in perturbation theory.

Interaccion among systems of finite size in predictive relativistic mechanics -II. Electromagnetic and gravitational interactions

We explicitly construct a closed system of differential equations describing the electromagnetic and gravitational interactions among bodies to first order in the coupling constants, retaining terms up to order c-2. The Breit and Barker and O'Connell Hamiltonians are recovered by means of a coordinate transformation. The method used throws light on the meaning of these coordinates.

Interactions among systems of finite size en predictive relativistic mechanics III. Short-range interaction and second-order terms

We compute up to and including all the c-2 terms in the dynamical equations for extended bodies interacting through electromagnetic, gravitational, or short-range fields. We show that these equations can be reduced to those of point particles with intrinsic angular momentum assuming spherical symmetry.

Efeito do grau de moagem, do tipo de frasco e do volume vazio sobre a variabilidade analítica do fósforo extraído pelos métodos Mehlich-1 e Mehlich-3

Na execução de uma análise química de solo são empregados diversos procedimentos que, mesmo seguindo-se o protocolo preconizado pelo método de análise utilizado, estão sujeitos a variações nos resultados analíticos, causadas por manipulação das amostras ou dos materiais utilizados. Objetivou-se neste estudo avaliar a influência do grau de moagem da amostra, do tipo de frasco e do volume vazio no frasco na execução dos métodos Mehlich-1 e Mehlich-3 para determinar o P no solo. Para tanto, foram conduzidos três experimentos. No primeiro experimento, as amostras de solo foram moídas e tamisadas em peneiras com aberturas de 2,000; 1,700; 0,850; 0,600; e 0,300 mm. No segundo, foram utilizados dois modelos de frasco (erlenmeyer e snap-cap), ambos com volume de 50 mL. Já no terceiro, foi alterado o volume vazio no frasco, mantendo-se a relação solo:solução extratora utilizando-se as quantidades dentro do frasco de 1:10; 1,5:15; 2,5:25; 3:30; e 4:40 cm³ cm-3. O grau de moagem das amostras não influenciou a capacidade de extração do Mehlich-1; entretanto, a capacidade extrativa do Mehlich-3 foi influenciada, principalmente em solos argilosos. Tanto para o Mehlich-1 quanto para o Mehlich-3, os teores de P extraído foram significativamente mais elevados com o uso de frasco tipo snap-cap em relação ao erlenmeyer. O volume vazio no frasco alterou os teores de P extraído para o Mehlich-1 e Mehlich-3 em 100 e 64 % das amostras, respectivamente. Deve-se padronizar a intensidade da moagem das amostras de solo para extração do P pela solução de Mehlich-3. Um modelo único de frasco deve ser adotado pelos laboratórios de rotina para análise do P, independentemente do método de extração, mantendo-se sempre constante no frasco o volume da amostra (cm³) para o volume de solução extratora (cm³).

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.

Fluctuations in finite systems: fluctuations in fluids in thermocapillary motion

We compute nonequilibrium correlation functions about the stationary state in which the fluid moves as a consequence of tangential stresses on the liquid surface, related to a varying surface tension (thermocapillary motion). The nature of the stationary state makes it necessary to take into account that the system is finite. We then extend a previous analysis on fluctuations about simple stationary states to include some effects related to the finite size of the sample.

Comment on "Canonical formalism for Lagrangians with nonlocality of finite extent"

The paper by Woodward [Phys. Rev. A 62, 052105 (2000)] claimed to have proved that Lagrangian theories with a nonlocality of finite extent are necessarily unstable. In this Comment we propose that this conclusion is false.

Finite-size effects in nonequilibrium correlation functions

We have shown that finite-size effects in the correlation functions away from equilibrium may be introduced through dimensionless numbers: the Nusselt numbers, accounting for both the nature of the boundaries and the size of the system. From an analysis based on fluctuating hydrodynamics, we conclude that the mean-square fluctuations satisfy scaling laws, since they depend only on the dimensionless numbers in addition to reduced variables. We focus on the case of diffusion modes and describe some physical situations in which finite-size effects may be relevant.

Adsorption kinetics in the presence of external fields

In this article we present a detailed analysis of the kinetics of a class of sequential adsorption models that take into account the effect of externally applied fields (as an electric field, or a shear rate) on the adsorption. The excluded volume interactions related to the finite size of the adsorbing particles are modified by the external fields. As a result, new adsorption mechanisms appear with respect to the ones used to describe the kinetics in a quiescent fluid. In particular, if the adsorbing particles are allowed to roll over preadsorbed ones, adsorption becomes non local even in the simplest geometry. An exact analytic theory cannot be developed, but we introduce a self-consistent theory that turns out to agree with the simulation results over all the range of the parameters.

Integrated random processes exhibiting long tails, finite moments, and power-law spectra

A dynamical model based on a continuous addition of colored shot noises is presented. The resulting process is colored and non-Gaussian. A general expression for the characteristic function of the process is obtained, which, after a scaling assumption, takes on a form that is the basis of the results derived in the rest of the paper. One of these is an expansion for the cumulants, which are all finite, subject to mild conditions on the functions defining the process. This is in contrast with the Lévy distribution¿which can be obtained from our model in certain limits¿which has no finite moments. The evaluation of the spectral density and the form of the probability density function in the tails of the distribution shows that the model exhibits a power-law spectrum and long tails in a natural way. A careful analysis of the characteristic function shows that it may be separated into a part representing a Lévy process together with another part representing the deviation of our model from the Lévy process. This

Adsorption kinetics in the presence of external fields

In this article we present a detailed analysis of the kinetics of a class of sequential adsorption models that take into account the effect of externally applied fields (as an electric field, or a shear rate) on the adsorption. The excluded volume interactions related to the finite size of the adsorbing particles are modified by the external fields. As a result, new adsorption mechanisms appear with respect to the ones used to describe the kinetics in a quiescent fluid. In particular, if the adsorbing particles are allowed to roll over preadsorbed ones, adsorption becomes non local even in the simplest geometry. An exact analytic theory cannot be developed, but we introduce a self-consistent theory that turns out to agree with the simulation results over all the range of the parameters.

Spatial analysis of erythrocyte membrane fluctuations by digital holographic microscopy.

Red blood cell (RBC) membrane fluctuations provide important insights into cell states. We present a spatial analysis of red blood cell membrane fluctuations by using digital holographic microscopy (DHM). This interferometric and dye-free technique, possessing nanometric axial and microsecond temporal sensitivities enables to measure cell membrane fluctuations (CMF) on the whole cell surface. DHM acquisition is combined with a model which allows extracting the membrane fluctuation amplitude, while taking into account cell membrane topology. Uneven distribution of CMF amplitudes over the RBC surface is observed, showing maximal values in a ring corresponding to the highest points on the RBC torus as well as in some scattered areas in the inner region of the RBC. CMF amplitudes of 35.9+/-8.9 nm and 4.7+/-0.5 nm (averaged over the cell surface) were determined for normal and ethanol-fixed RBCs, respectively.

Monte Carlo study of the finite-size effects on the magnetization of maghemite small particles

Monte Carlo simulations of a model for gamma-Fe2O3 (maghemite) single particle of spherical shape are presented aiming at the elucidation of the specific role played by the finite size and the surface on the anomalous magnetic behavior observed in small particle systems at low temperature. The influence of the finite-size effects on the equilibrium properties of extensive magnitudes, field coolings, and hysteresis loops is studied and compared to the results for periodic boundaries. It is shown that for the smallest sizes the thermal demagnetization of the surface completely dominates the magnetization while the behavior of the core is similar to that of the periodic boundary case, independently of D. The change in shape of the hysteresis loops with D demonstrates that the reversal mode is strongly influenced by the presence of broken links and disorder at the surface

Finite-size scaling investigation of the liquid-liquid critical point in ST2 water and its stability with respect to crystallization.

The liquid-liquid critical point scenario of water hypothesizes the existence of two metastable liq- uid phases low-density liquid (LDL) and high-density liquid (HDL) deep within the supercooled region. The hypothesis originates from computer simulations of the ST2 water model, but the stabil- ity of the LDL phase with respect to the crystal is still being debated. We simulate supercooled ST2 water at constant pressure, constant temperature, and constant number of molecules N for N ≤ 729 and times up to 1 μs. We observe clear differences between the two liquids, both structural and dynamical. Using several methods, including finite-size scaling, we confirm the presence of a liquid-liquid phase transition ending in a critical point. We find that the LDL is stable with respect to the crystal in 98% of our runs (we perform 372 runs for LDL or LDL-like states), and in 100% of our runs for the two largest system sizes (N = 512 and 729, for which we perform 136 runs for LDL or LDL-like states). In all these runs, tiny crystallites grow and then melt within 1 μs. Only for N ≤ 343 we observe six events (over 236 runs for LDL or LDL-like states) of spontaneous crystal- lization after crystallites reach an estimated critical size of about 70 ± 10 molecules.