Nanocomposites: the contribution of multiscale modelling to improve the dispersion levels

Tese de Doutoramento em Ciência e Engenharia de Polímeros e Compósitos

Predicting optimal lengths of random knots

In a thermally fluctuating long linear polymeric chain in a solution, the ends, from time to time, approach each other. At such an instance, the chain can be regarded as closed and thus will form a knot or rather a virtual knot. Several earlier studies of random knotting demonstrated that simpler knots show a higher occurrence for shorter random walks than do more complex knots. However, up to now there have been no rules that could be used to predict the optimal length of a random walk, i.e. the length for which a given knot reaches its highest occurrence. Using numerical simulations, we show here that a power law accurately describes the relation between the optimal lengths of random walks leading to the formation of different knots and the previously characterized lengths of ideal knots of a corresponding type.

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.

Gravity-driven slug motion in capillary tubes

The velocity of a liquid slug falling in a capillary tube is lower than predicted for Poiseuille flow due to presence of menisci, whose shapes are determined by the complex interplay of capillary, viscous, and gravitational forces. Due to the presence of menisci, a capillary pressure proportional to surface curvature acts on the slug and streamlines are bent close to the interface, resulting in enhanced viscous dissipation at the wedges. To determine the origin of drag-force increase relative to Poiseuille flow, we compute the force resultant acting on the slug by integrating Navier-Stokes equations over the liquid volume. Invoking relationships from differential geometry we demonstrate that the additional drag is due to viscous forces only and that no capillary drag of hydrodynamic origin exists (i.e., due to hydrodynamic deformation of the interface). Requiring that the force resultant is zero, we derive scaling laws for the steady velocity in the limit of small capillary numbers by estimating the leading order viscous dissipation in the different regions of the slug (i.e., the unperturbed Poiseuille-like bulk, the static menisci close to the tube axis and the dynamic regions close to the contact lines). Considering both partial and complete wetting, we find that the relationship between dimensionless velocity and weight is, in general, nonlinear. Whereas the relationship obtained for complete-wetting conditions is found in agreement with the experimental data of Bico and Quere [J. Bico and D. Quere, J. Colloid Interface Sci. 243, 262 (2001)], the scaling law under partial-wetting conditions is validated by numerical simulations performed with the Volume of Fluid method. The simulated steady velocities agree with the behavior predicted by the theoretical scaling laws in presence and in absence of static contact angle hysteresis. The numerical simulations suggest that wedge-flow dissipation alone cannot account for the entire additional drag and that the non-Poiseuille dissipation in the static menisci (not considered in previous studies) has to be considered for large contact angles.

Stirring and mixing in two-dimensional divergent flow

While stirring and mixing properties in the stratosphere are reasonably well understood in the context of balanced (slow) dynamics, as is evidenced in numerous studies of chaotic advection, the strongly enhanced presence of high-frequency gravity waves in the mesosphere gives rise to a significant unbalanced (fast) component to the flow. The present investigation analyses result from two idealized shallow-water numerical simulations representative of stratospheric and mesospheric dynamics on a quasi-horizontal isentropic surface. A generalization of the Hua–Klein Eulerian diagnostic to divergent flow reveals that velocity gradients are strongly influenced by the unbalanced component of the flow. The Lagrangian diagnostic of patchiness nevertheless demonstrates the persistence of coherent features in the zonal component of the flow, in contrast to the destruction of coherent features in the meridional component. Single-particle statistics demonstrate t2 scaling for both the stratospheric and mesospheric regimes in the case of zonal dispersion, and distinctive scaling laws for the two regimes in the case of meridional dispersion. This is in contrast to two-particle statistics, which in the mesospheric (unbalanced) regime demonstrate a more rapid approach to Richardson’s t3 law in the case of zonal dispersion and is evidence of enhanced meridional dispersion.

Regimes of conformational transitions of a diblock polyampholyte

The conformational properties of symmetric flexible diblock polyampholytes are investigated by scaling theory and molecular dynamics simulations. The electrostatically driven coil-globule transition of a symmetric diblock polyampholyte is found to consist of three regimes identified with increasing electrostatic interaction strength. In the first (folding) regime the electrostatic attraction causes the chain to fold through the overlap of the two blocks, while each block is slightly stretched by self-repulsion. The second (weak association or scrambled egg) regime is the classical collapse of the chain into a globule dominated by the fluctuation-induced attractions between oppositely charged sections of the chain. The structure of the formed globule can be represented as a dense packing of the charged chain sections (electrostatic attraction blobs). The third (strong association or ion binding) regime starts with direct binding of oppositely charged monomers (dipole formation), followed by a cascade of multipole formation (quadrupole, hexapole, octupole, etc.), leading to multiplets analogous to those found in ionomers. The existence of the multiplet cascade has also been confirmed in the simulations of solutions of short polymers with only one single charge (either positive or negative) in the middle of each chain. We use scaling theory to estimate the average chain size and the electrostatic correlation length as functions of the chain length, strength of electrostatic interactions, charge fraction, and solvent quality. The theoretically predicted scaling laws of these conformational properties are in very good agreement with our simulation results.

A scenario for torus T(2) destruction via a global bifurcation

We show a scenario of a two-frequeney torus breakdown, in which a global bifurcation occurs due to the collision of a quasi-periodic torus T(2) with saddle points, creating a heteroclinic saddle connection. We analyze the geometry of this torus-saddle collision by showing the local dynamics and the invariant manifolds (global dynamics) of the saddle points. Moreover, we present detailed evidences of a heteroclinic saddle-focus orbit responsible for the type-if intermittency induced by this global bifurcation. We also characterize this transition to chaos by measuring the Lyapunov exponents and the scaling laws. (C) 2007 Elsevier Ltd. All rights reserved.

Lie algebra methods for the applications to the statistical theory of turbulence

Approximate Lie symmetries of the Navier-Stokes equations are used for the applications to scaling phenomenon arising in turbulence. In particular, we show that the Lie symmetries of the Euler equations are inherited by the Navier-Stokes equations in the form of approximate symmetries that allows to involve the Reynolds number dependence into scaling laws. Moreover, the optimal systems of all finite-dimensional Lie subalgebras of the approximate symmetry transformations of the Navier-Stokes are constructed. We show how the scaling groups obtained can be used to introduce the Reynolds number dependence into scaling laws explicitly for stationary parallel turbulent shear flows. This is demonstrated in the framework of a new approach to derive scaling laws based on symmetry analysis [11]-[13].

Simulação de fluxo de fluidos em meios porosos desordenados uma análise de efeito de escala na estimativa da permeabilidade e do coeficiente de arrasto

The present study provides a methodology that gives a predictive character the computer simulations based on detailed models of the geometry of a porous medium. We using the software FLUENT to investigate the flow of a viscous Newtonian fluid through a random fractal medium which simplifies a two-dimensional disordered porous medium representing a petroleum reservoir. This fractal model is formed by obstacles of various sizes, whose size distribution function follows a power law where exponent is defined as the fractal dimension of fractionation Dff of the model characterizing the process of fragmentation these obstacles. They are randomly disposed in a rectangular channel. The modeling process incorporates modern concepts, scaling laws, to analyze the influence of heterogeneity found in the fields of the porosity and of the permeability in such a way as to characterize the medium in terms of their fractal properties. This procedure allows numerically analyze the measurements of permeability k and the drag coefficient Cd proposed relationships, like power law, for these properties on various modeling schemes. The purpose of this research is to study the variability provided by these heterogeneities where the velocity field and other details of viscous fluid dynamics are obtained by solving numerically the continuity and Navier-Stokes equations at pore level and observe how the fractal dimension of fractionation of the model can affect their hydrodynamic properties. This study were considered two classes of models, models with constant porosity, MPC, and models with varying porosity, MPV. The results have allowed us to find numerical relationship between the permeability, drag coefficient and the fractal dimension of fractionation of the medium. Based on these numerical results we have proposed scaling relations and algebraic expressions involving the relevant parameters of the phenomenon. In this study analytical equations were determined for Dff depending on the geometrical parameters of the models. We also found a relation between the permeability and the drag coefficient which is inversely proportional to one another. As for the difference in behavior it is most striking in the classes of models MPV. That is, the fact that the porosity vary in these models is an additional factor that plays a significant role in flow analysis. Finally, the results proved satisfactory and consistent, which demonstrates the effectiveness of the referred methodology for all applications analyzed in this study.

Não-extensividade em percolação por ligações de longo - alcance

A linear chain do not present phase transition at any finite temperature in a one dimensional system considering only first neighbors interaction. An example is the Ising ferromagnet in which his critical temperature lies at zero degree. Analogously, in percolation like disordered geometrical systems, the critical point is given by the critical probability equals to one. However, this situation can be drastically changed if we consider long-range bonds, replacing the probability distribution by a function like . In this kind of distribution the limit α → ∞ corresponds to the usual first neighbor bond case. In the other hand α = 0 corresponds to the well know "molecular field" situation. In this thesis we studied the behavior of Pc as a function of a to the bond percolation specially in d = 1. Our goal was to check a conjecture proposed by Tsallis in the context of his Generalized Statistics (a generalization to the Boltzmann-Gibbs statistics). By this conjecture, the scaling laws that depend with the size of the system N, vary in fact with the quantitie

Excitações elementares em super-redes quasiperiódicas com simetria espelho

In this thesis, we investigated the magnonic and photonic structures that exhibit the so-called deterministic disorder. Speci cally, we studied the effects of the quasiperiodicity, associated with an internal structural symmetry, called mirror symmetry, on the spectra of photonics and magnonics multilayer. The quasiperiodicity is introduced when stacked layers following the so-called substitutional sequences. The three sequences used here were the Fibonacci sequence, Thue-Morse and double-period, all with mirror symmetry. Aiming to study the propagation of light waves in multilayer photonic, and spin waves propagation in multilayer magnonic, we use a theoretical model based on transfer matrix treatment. For the propagation of light waves, we present numerical results that show that the quasiperiodicity associated with a mirror symmetry greatly increases the intensity of transmission and the transmission spectra exhibit a pro le self-similar. The return map plotted for this system show that the presence of internal symmetry does not alter the pattern of Fibonacci maps when compared with the case without symmetry. But when comparing the maps of Thue-Morse and double-time sequences with their case without the symmetry mirror, is evident the change in the pro le of the maps. For magnetic multilayers, we work with two di erent systems, multilayer composed of a metamagnetic material and a non-magnetic material, and multilayers composed of two cubic Heisenberg ferromagnets. In the rst case, our calculations are carried out in the magnetostatic regime and calculate the dispersion relation of spin waves for the metamgnetic material considered FeBr2. We show the e ect of mirror symmetry in the spectra of spin waves, and made the analysis of the location of bulk bands and the scaling laws between the full width of the bands allowed and the number of layers of unit cell. Finally, we calculate the transmission spectra of spin waves in quasiperiodic multilayers consisting of Heisenberg ferromagnets. The transmission spectra exhibit self-similar patterns, with regions of scaling well-de ned in frequency and the return maps indicates only dependence of the particular sequence used in the construction of the multilayer

Universal aspects of light halo nuclei

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Phase Transition in Dynamical Systems: Defining Classes of Universality for Two-Dimensional Hamiltonian Mappings via Critical Exponents

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Suppressing Fermi acceleration in two-dimensional driven billiards

We consider a dissipative oval-like shaped billiard with a periodically moving boundary. The dissipation considered is proportional to a power of the velocity V of the particle. The three specific types of power laws used are: (i) F proportional to-V; (ii) F proportional to-V-2 and (iii) F proportional to-V-delta with 1 < delta < 2. In the course of the dynamics of the particle, if a large initial velocity is considered, case (i) shows that the decay of the particle's velocity is a linear function of the number of collisions with the boundary. For case (ii), an exponential decay is observed, and for 1 < delta < 2, an powerlike decay is observed. Scaling laws were used to characterize a phase transition from limited to unlimited energy gain for cases (ii) and (iii). The critical exponents obtained for the phase transition in the case (ii) are the same as those obtained for the dissipative bouncer model. Therefore near this phase transition, these two rather different models belong to the same class of universality. For all types of dissipation, the results obtained allow us to conclude that suppression of the unlimited energy growth is indeed observed.