O problema da matéria escura e o teorema do virial na teoria de gravidade Eddington-Born-Infeld

Dark matter is a fundamental ingredient of the modern Cosmology. It is necessary in order to explain the process of structures formation in the Universe, rotation curves of galaxies and the mass discrepancy in clusters of galaxies. However, although many efforts, in both aspects, theoretical and experimental, have been made, the nature of dark matter is still unknown and the only convincing evidence for its existence is gravitational. This rises doubts about its existence and, in turn, opens the possibility that the Einstein’s gravity needs to be modified at some scale. We study, in this work, the possibility that the Eddington-Born-Infeld (EBI) modified gravity provides en alternative explanation for the mass discrepancy in clusters of galaxies. For this purpose we derive the modified Einstein field equations and find their solutions to a spherical system of identical and collisionless point particles. Then, we took into account the collisionless relativistic Boltzmann equation and using some approximations and assumptions for weak gravitational field, we derived the generalized virial theorem in the framework of EBI gravity. In order to compare the predictions of EBI gravity with astrophysical observations we estimated the order of magnitude of the geometric mass, showing that it is compatible with present observations. Finally, considering a power law for the density of galaxies in the cluster, we derived expressions for the radial velocity dispersion of the galaxies, which can be used for testing some features of the EBI gravity.

Simulations of charged droplet collisions in shear flow

Acknowledgments This research has been enabled by the use of computing resources provided by WestGrid, the Shared Hierarchical Academic Research Computing Network (SHARCNET: www.sharcnet.ca), and Compute/Calcul Canada. O.S. thanks NSERC for an Alexander Graham Bell Canada Graduate Scholarship.

Discrete Velocity Models for Mixtures Without Nonphysical Collision Invariants

An important aspect of constructing discrete velocity models (DVMs) for the Boltzmann equation is to obtain the right number of collision invariants. It is a well-known fact that DVMs can also have extra collision invariants, so called spurious collision invariants, in plus to the physical ones. A DVM with only physical collision invariants, and so without spurious ones, is called normal. For binary mixtures also the concept of supernormal DVMs was introduced, meaning that in addition to the DVM being normal, the restriction of the DVM to any single species also is normal. Here we introduce generalizations of this concept to DVMs for multicomponent mixtures. We also present some general algorithms for constructing such models and give some concrete examples of such constructions. One of our main results is that for any given number of species, and any given rational mass ratios we can construct a supernormal DVM. The DVMs are constructed in such a way that for half-space problems, as the Milne and Kramers problems, but also nonlinear ones, we obtain similar structures as for the classical discrete Boltzmann equation for one species, and therefore we can apply obtained results for the classical Boltzmann equation.

A lattice Boltzmann method for KDV equation

We prepose a 5-bit lattice Boltzmann model for KdV equation. Using Chapman-Enskog expansion and multiscale technique, we obtained high order moments of equilibrium distribution function, and the 3rd dispersion coefficient and 4th order viscosity. The parameters of this scheme can be determined by analysing the energy dissipation.

Spectrum reconstruction from a scattering measurement using the adjoint Boltzmann transport equation for photons

Quality control of medical radiological systems is of fundamental importance, and requires efficient methods for accurately determine the X-ray source spectrum. Straightforward measurements of X-ray spectra in standard operating require the limitation of the high photon flux, and therefore the measure has to be performed in a laboratory. However, the optimal quality control requires frequent in situ measurements which can be only performed using a portable system. To reduce the photon flux by 3 magnitude orders an indirect technique based on the scattering of the X-ray source beam by a solid target is used. The measured spectrum presents a lack of information because of transport and detection effects. The solution is then unfolded by solving the matrix equation that represents formally the scattering problem. However, the algebraic system is ill-conditioned and, therefore, it is not possible to obtain a satisfactory solution. Special strategies are necessary to circumvent the ill-conditioning. Numerous attempts have been done to solve this problem by using purely mathematical methods. In this thesis, a more physical point of view is adopted. The proposed method uses both the forward and the adjoint solutions of the Boltzmann transport equation to generate a better conditioned linear algebraic system. The procedure has been tested first on numerical experiments, giving excellent results. Then, the method has been verified with experimental measurements performed at the Operational Unit of Health Physics of the University of Bologna. The reconstructed spectra have been compared with the ones obtained with straightforward measurements, showing very good agreement.

Fully Anisotropic Solution of the Three Dimensional Boltzmann Transport Equation

The development of accurate modeling techniques for nanoscale thermal transport is an active area of research. Modern day nanoscale devices have length scales of tens of nanometers and are prone to overheating, which reduces device performance and lifetime. Therefore, accurate temperature profiles are needed to predict the reliability of nanoscale devices. The majority of models that appear in the literature obtain temperature profiles through the solution of the Boltzmann transport equation (BTE). These models often make simplifying assumptions about the nature of the quantized energy carriers (phonons). Additionally, most previous work has focused on simulation of planar two dimensional structures. This thesis presents a method which captures the full anisotropy of the Brillouin zone within a three dimensional solution to the BTE. The anisotropy of the Brillouin zone is captured by solving the BTE for all vibrational modes allowed by the Born Von-Karman boundary conditions.

Multicarrier conduction and Boltzmann transport analysis of heavy hole mobility in HgCdTe near room temperature

Magnetotransport measurements in pulsed fields up to 15 T have been performed on mercury cadmium telluride (Hg1-xCdxTe, x similar to 0.2) bulk as well as liquid phase epitaxially grown samples to obtain the resistivity and conductivity tensors in the temperature range 220-300 K. Mobilities and densities of various carriers participating in conduction have been extracted using both conventional multicarrier fitting (MCF) and mobility spectrum analysis. The fits to experimental data, particularly at the highest magnetic fields, were substantially improved when MCF is applied to minimize errors simultaneously on both resistivity and conductivity tensors. The semiclassical Boltzmann transport equation has been solved without using adjustable parameters by incorporating the following scattering mechanisms to fit the mobility: ionized impurity, polar and nonpolar optical phonons, acoustic deformation potential, and alloy disorder. Compared to previous estimates based on the relaxation time approximation with outscattering only, polar optical scattering and ionized impurity scattering limited mobilities are shown to be larger due to the correct incorporation of the inscattering term taking into account the overlap integrals in the valence band.

Lattice Boltzmann study of hydrodynamic effects in lamellar ordering process of two-dimensional quenched block copolymers

By incorporating self-consistent field theory with lattice Boltzmann method, a model for polymer melts is proposed. Compared with models based on Ginzburg-Landau free energy, our model does not employ phenomenological free energies to describe systems and can consider the chain topological details of polymers. We use this model to study the effects of hydrodynamic interactions on the dynamics of microphase separation for block copolymers. In the early stage of phase separation, an exponential growth predicted by Cahn-Hilliard treatment is found. Simulation results also show that the effect of hydrodynamic interactions can be neglected in the early stage.

Calculating the equation of state parameters and predicting the spinodal curve of isotactic polypropylene/poly(ethylene-co-octene) blend by molecular dynamics simulations combined with Sanchez-Lacombe lattice fluid theory

Molecular dynamics simulations are adopted to calculate the equation of state characteristic parameters P*, rho*, and T* of isotactic polypropylene (iPP) and poly(ethylene-co-octene) (PEOC), which can be further used in the Sanchez-Lacombe lattice fluid theory (SLLFT) to describe the respective physical properties. The calculated T* is a function of the temperature, which was also found in the literature. To solve this problem, we propose a Boltzmann fitting of the data and obtain T* at the high-temperature limit. With these characteristic parameters, the pressure-volume-temperature (PVT) data of iPP and PEOC are predicted by the SLLFT equation of state. To justify the correctness of our results, we also obtain the PVT data for iPP and PEOC by experiments. Good agreement is found between the two sets of data. By integrating the Euler-Lagrange equation and the Cahn-Hilliard relation, we predict the density profiles and the surface tensions for iPP and PEOC, respectively. Furthermore, a recursive method is proposed to obtain the characteristic interaction energy parameter between iPP and PEOC. This method, which does not require fitting to the experimental phase equilibrium data, suggests an alternative way to predict the phase diagrams that are not easily obtained in experiments.

A Van der Pol-Mathieu equation for the dynamics of dust grain charge in dusty plasmas

The chaotic profile of dust grain dynamics associated with dust-acoustic oscillations in a dusty plasma is considered. The collective behaviour of the dust plasma component is described via a multi-fluid model, comprising Boltzmann distributed electrons and ions, as well as an equation of continuity possessing a source term for the dust grains, the dust momentum and Poisson's equations. A Van der Pol–Mathieu-type nonlinear ordinary differential equation for the dust grain density dynamics is derived. The dynamical system is cast into an autonomous form by employing an averaging method. Critical stability boundaries for a particular trivial solution of the governing equation with varying parameters are specified. The equation is analysed to determine the resonance region, and finally numerically solved by using a fourth-order Runge–Kutta method. The presence of chaotic limit cycles is pointed out.

pH at the micellar interface: Synthesis of pH probes derived from salicylic acid, acid-base dissociation in sodium dodecyl sulfate micelles, and Poisson-Boltzmann simulation

The study of the H+ concentration at the micellar interface is a convenient system for modeling the distribution of H+ at interfaces. We have synthesized salicylic acid derivatives to analyze the proton dissociation of both the carboxylic and phenol groups of' the probes, determining spectrophotometrically the apparent pK(a)'s (pK(ap)) in sodium dodecyl Sulfate, SDS, micelles with and without added salt. The synthesized probes were 2-hydroxy-5-(2-trimethylammoniumacetyl)benzoate; 2-hydroxy-5-(2-dimethylhexadecylammoniumacetyl)benzoate- 2-hydroxy-5-(2-dimethylhexadecylammoniumhexanoyl)benzoate-, 2-hydroxy-5-(2-diniethylhexadecylammoniumundecanoyl)betizoate; 2-hydroxy-5-acetylbenzoic acids and 2-hydroxy-5-dodecanoylbenzoic acid. Upon incorporation into SDS micelles the pK(ap)'s of both carboxylic and phenol groups increased by ca. 3 pH units and NaCl addition caused a decrease in the probe-incorporated pKap. The experimental results were fitted with a cell model Poisson-Boltzmann (P-B) equation taking in consideration the effect of salt on the aggregation number of SDS and using the distance of' the dissociating group as a parameter. The conformations of the probes were analyzed theoretically using two dielectric constants, e.g., 2 and 78. Both the P-B analysis and conformation calculations can be interpreted by assuming that the acid groups dissociate very close to, or at, the interface. Our results are consistent with the assumption that the intrinsic pK(a)'s of both carboxylic and phenol groups of the salicylic acid probes used here can be taken as those in water. Using this assumption the micellar and salt effects on the pKap's of the (trialkylammonium)benzoate probes were described accurately using a cell model P-B analysis. (c) 2005 Elsevier B.V. All rights reserved.

Generalized Metropolis dynamics with a generalized master equation: An approach for time-independent and time-dependent Monte Carlo simulations of generalized spin systems

The extension of Boltzmann-Gibbs thermostatistics, proposed by Tsallis, introduces an additional parameter q to the inverse temperature beta. Here, we show that a previously introduced generalized Metropolis dynamics to evolve spin models is not local and does not obey the detailed energy balance. In this dynamics, locality is only retrieved for q = 1, which corresponds to the standard Metropolis algorithm. Nonlocality implies very time-consuming computer calculations, since the energy of the whole system must be reevaluated when a single spin is flipped. To circumvent this costly calculation, we propose a generalized master equation, which gives rise to a local generalized Metropolis dynamics that obeys the detailed energy balance. To compare the different critical values obtained with other generalized dynamics, we perform Monte Carlo simulations in equilibrium for the Ising model. By using short-time nonequilibrium numerical simulations, we also calculate for this model the critical temperature and the static and dynamical critical exponents as functions of q. Even for q not equal 1, we show that suitable time-evolving power laws can be found for each initial condition. Our numerical experiments corroborate the literature results when we use nonlocal dynamics, showing that short-time parameter determination works also in this case. However, the dynamics governed by the new master equation leads to different results for critical temperatures and also the critical exponents affecting universality classes. We further propose a simple algorithm to optimize modeling the time evolution with a power law, considering in a log-log plot two successive refinements.

Lattice Boltzmann modeling of fluid flow and solute transport in karst aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Lattice boltzmann modeling of fluid flow to determine the permeability of a karst specimen

A limestone sample was scanned using computed tomography (CT) and the hydraulic conductivity of the 3D reconstructed sample was determined using Lattice- Boltzmann methods (LBM) at varying scales. Due to the shape and size of the original sample, it was challenging to obtain a consistent rectilinear test sample. Through visual inspection however, 91 mm and 76 mm samples were digitally cut from the original. The samples had porosities of 58% and 64% and produced hydraulic conductivity values of K= 13.5 m/s and K=34.5 m/s, respectively. Both of these samples were re-sampled to 1/8 and 1/64 of their original size to produce new virtual samples at lower resolutions of 0.542 mm/lu and 1.084 mm/lu, while still representing the same physical dimensions. The hydraulic conductivity tended to increase slightly as the resolution became coarser. In order to determine an REV, the 91 mm sample was also sub-sampled into blocks that were 1/8 and 1/64 the size of the original. The results were consistent with analytical expectations such as those produced by the Kozeny-Carman equation. A definitive REV size was not reached, however, indicating the need for a larger sample. The methods described here demonstrate the ability of LBM to test rock structures and sizes not normally attainable.