Propriedades dinâmicas de fluidos por simulação computacional: métodos híbridos atomístico-contínuo

Computational methods for the calculation of dynamical properties of fluids might consider the system as a continuum or as an assembly of molecules. Molecular dynamics (MD) simulation includes molecular resolution, whereas computational fluid dynamics (CFD) considers the fluid as a continuum. This work provides a review of hybrid methods MD/CFD recently proposed in the literature. Theoretical foundations, basic approaches of computational methods, and dynamical properties typically calculated by MD and CFD are first presented in order to appreciate the similarities and differences between these two methods. Then, methods for coupling MD and CFD, and applications of hybrid simulations MD/CFD, are presented.

Turbulence in collisionless plasmas: statistical analysis from numerical simulations with pressure anisotropy

In recent years, we have experienced increasing interest in the understanding of the physical properties of collisionless plasmas, mostly because of the large number of astrophysical environments (e. g. the intracluster medium (ICM)) containing magnetic fields that are strong enough to be coupled with the ionized gas and characterized by densities sufficiently low to prevent the pressure isotropization with respect to the magnetic line direction. Under these conditions, a new class of kinetic instabilities arises, such as firehose and mirror instabilities, which have been studied extensively in the literature. Their role in the turbulence evolution and cascade process in the presence of pressure anisotropy, however, is still unclear. In this work, we present the first statistical analysis of turbulence in collisionless plasmas using three-dimensional numerical simulations and solving double-isothermal magnetohydrodynamic equations with the Chew-Goldberger-Low laws closure (CGL-MHD). We study models with different initial conditions to account for the firehose and mirror instabilities and to obtain different turbulent regimes. We found that the CGL-MHD subsonic and supersonic turbulences show small differences compared to the MHD models in most cases. However, in the regimes of strong kinetic instabilities, the statistics, i.e. the probability distribution functions (PDFs) of density and velocity, are very different. In subsonic models, the instabilities cause an increase in the dispersion of density, while the dispersion of velocity is increased by a large factor in some cases. Moreover, the spectra of density and velocity show increased power at small scales explained by the high growth rate of the instabilities. Finally, we calculated the structure functions of velocity and density fluctuations in the local reference frame defined by the direction of magnetic lines. The results indicate that in some cases the instabilities significantly increase the anisotropy of fluctuations. These results, even though preliminary and restricted to very specific conditions, show that the physical properties of turbulence in collisionless plasmas, as those found in the ICM, may be very different from what has been largely believed.

Unconstrained Finite Element for Geometrical Nonlinear Dynamics of Shells

This paper presents a positional FEM formulation to deal with geometrical nonlinear dynamics of shells. The main objective is to develop a new FEM methodology based on the minimum potential energy theorem written regarding nodal positions and generalized unconstrained vectors not displacements and rotations. These characteristics are the novelty of the present work and avoid the use of large rotation approximations. A nondimensional auxiliary coordinate system is created, and the change of configuration function is written following two independent mappings from which the strain energy function is derived. This methodology is called positional and, as far as the authors' knowledge goes, is a new procedure to approximated geometrical nonlinear structures. In this paper a proof for the linear and angular momentum conservation property of the Newmark beta algorithm is provided for total Lagrangian description. The proposed shell element is locking free for elastic stress-strain relations due to the presence of linear strain variation along the shell thickness. The curved, high-order element together with an implicit procedure to solve nonlinear equations guarantees precision in calculations. The momentum conserving, the locking free behavior, and the frame invariance of the adopted mapping are numerically confirmed by examples. Copyright (C) 2009 H. B. Coda and R. R. Paccola.

A positional FEM Formulation for geometrical non-linear analysis of shells

This work presents a fully non-linear finite element formulation for shell analysis comprising linear strain variation along the thickness of the shell and geometrically exact description for curved triangular elements. The developed formulation assumes positions and generalized unconstrained vectors as the variables of the problem, not displacements and finite rotations. The full 3D Saint-Venant-Kirchhoff constitutive relation is adopted and, to avoid locking, the rate of thickness variation enhancement is introduced. As a consequence, the second Piola-Kirchhoff stress tensor and the Green strain measure are employed to derive the specific strain energy potential. Curved triangular elements with cubic approximation are adopted using simple notation. Selected numerical simulations illustrate and confirm the objectivity, accuracy, path independence and applicability of the proposed technique.

Pulse train fluorescence technique for measuring triplet state dynamics

We report on a method to study the dynamics of triplet formation based on the fluorescence signal produced by a pulse train. Basically, the pulse train acts as sequential pump-probe pulses that precisely map the excited-state dynamics in the long time scale. This allows characterizing those processes that affect the population evolution of the first excited singlet state, whose decay gives rise to the fluorescence. The technique was proven to be valuable to measure parameters of triplet formation in organic molecules. Additionally, this single beam technique has the advantages of simplicity, low noise and background-free signal detection. (C) 2011 Optical Society of America

Brazilian territorial dynamics and the inversion of the ""frontier thesis"" in the southern new world

It is a generally acknowledged fact that the dynamics of frontier advance deeply influenced the broad experience of American post colonial societies. The colonization, which started most from the east boundaries of the continent, appropriated and gradually transformed the American territories from east to west. The advance, initially represented by the arrival of the European settlers, went on to become an important trace of that society which did not come to know any physical limits of a restricted territory. However, despite the common identity granted by these territorial dynamics, the later developments and consequences seem to have shaped differently the Northern representatives from their Southern counterparts. In addition, the interpretation of these facts bore in each of these regions different meanings and traits.

Notes on Population Dynamics of Amblyomma Ticks (Acari: Ixodidae) in Brazil

Previous population dynamics data, generated for Amblyomma parvum Aragao and Aniblyomma cajennense (Fabricius) in Argentina and southeastern Brazil, have indicated that these ticks complete I generation per year, with larvae predominating in autumn, nymphs in winter, and mostly adults during spring and summer. The present study reports population dynamics data for free-living Amblyomma spp. ticks in northern Brazil (Amazon forest, latitude 10 degrees S, 63 degrees W), and for Amblyomma spp. ticks collected oil birds in Southeastern Brazil (latitude 23 degrees S, 45 degrees W). In northern Brazil, adult ticks predominated from mid-spring to mid-autumn, larvae predominated in early winter, and nymphs from mid-winter to mid-spring. Seven Amblyomma spp. were identified, although A. cajannense predominated in I of the 2 sites sampled. In southeastern Brazil, larval infestations on birds peaked in autumn, followed by a nymphal infestation peak in late winter. At least 32% and 75% of these larvae and nymphs, respectively, were identified as Amblyomma longirostre (Koch). Similar to previous work, the present study showed that Amblyomma spp. larvae and nymphs predominated during autumn-winter months, and mostly adults during spring-summer months, a pattern compatible with 1 genration/yr, even at latitude 10 degrees S in northern Brazil.

The structure and dynamics of Abell 1942

Aims. We present a dynamical analysis of the galaxy cluster Abell 1942 based on a set of 128 velocities obtained at the European Southern Observatory. Methods. Data on individual galaxies are presented and the accuracy of the determined velocities as some properties of the cluster are discussed. We have also made use of publicly available Chandra X-ray data. Results. We obtained an improved mean redshift value z = 0.22513 +/- 0.0008 and velocity dispersion sigma = 908(139)(+147) km s(-1). Our analysis indicates that inside a radius of similar to 1.5 h(70)(-1) Mpc (similar to 7 arcmin) the cluster is well relaxed, without any remarkable features and the X-ray emission traces the galaxy distribution fairly well. Two possible optical substructures are seen at similar to 5 arcmin from the centre in the northwest and the southwest directions, but are not confirmed by the velocity field. These clumps are, however, kinematically bound to the main structure of Abell 1942. X-ray spectroscopic analysis of Chandra data resulted in a temperature kT = 5.5+/-0.5 keV and metal abundance Z = 0.33 +/- 0.15 Z(circle dot). The velocity dispersion corresponding to this temperature using the T(X-sigma) scaling relation is in good agreement with the measured galaxy velocities. Our photometric redshift analysis suggests that the weak lensing signal observed to the south of the cluster and previously attributed to a ""dark clump"" is produced by background sources, possibly distributed as a filamentary structure.

Rotation curve bifurcations as indicators of close recent galaxy encounters

Context. Rotation curves of interacting galaxies often show that velocities are either rising or falling in the direction of the companion galaxy. Aims. We seek to reproduce and analyse these features in the rotation curves of simulated equal-mass galaxies suffering a one-to-one encounter as possible indicators of close encounters. Methods. Using simulations of major mergers in 3D, we study the time evolution of these asymmetries in a pair of galaxies during the first passage. Results. Our main results are: (a) the rotation curve asymmetries appear right at pericentre of the first passage, (b) the significant disturbed rotation velocities occur within a small time interval, of similar to 0.5 Gyr h(-1), and, therefore, the presence of bifurcation in the velocity curve could be used as an indicator of the pericentre occurrence. These results are in qualitative agreement with previous findings for minor mergers and flybys.

Dynamics and constraints of the massive graviton dark matter flat cosmologies

We discuss the dynamics of the Universe within the framework of the massive graviton cold dark matter scenario (MGCDM) in which gravitons are geometrically treated as massive particles. In this modified gravity theory, the main effect of the gravitons is to alter the density evolution of the cold dark matter component in such a way that the Universe evolves to an accelerating expanding regime, as presently observed. Tight constraints on the main cosmological parameters of the MGCDM model are derived by performing a joint likelihood analysis involving the recent supernovae type Ia data, the cosmic microwave background shift parameter, and the baryonic acoustic oscillations as traced by the Sloan Digital Sky Survey red luminous galaxies. The linear evolution of small density fluctuations is also analyzed in detail. It is found that the growth factor of the MGCDM model is slightly different (similar to 1-4%) from the one provided by the conventional flat Lambda CDM cosmology. The growth rate of clustering predicted by MGCDM and Lambda CDM models are confronted to the observations and the corresponding best fit values of the growth index (gamma) are also determined. By using the expectations of realistic future x-ray and Sunyaev-Zeldovich cluster surveys we derive the dark matter halo mass function and the corresponding redshift distribution of cluster-size halos for the MGCDM model. Finally, we also show that the Hubble flow differences between the MGCDM and the Lambda CDM models provide a halo redshift distribution departing significantly from the those predicted by other dark energy models. These results suggest that the MGCDM model can observationally be distinguished from Lambda CDM and also from a large number of dark energy models recently proposed in the literature.

Insights on the Milky Way bulge formation from the correlations between kinematics and metallicity

Context. Two main scenarios for the formation of the Galactic bulge are invoked, the first one through gravitational collapse or hierarchical merging of subclumps, the second through secular evolution of the Galactic disc. Aims. We aim to constrain the formation of the Galactic bulge through studies of the correlation between kinematics and metallicities in Baade's Window (l = 1 degrees, b = -4 degrees) and two other fields along the bulge minor axis (l = 0 degrees, b = -6 degrees and b = -12 degrees). Methods. We combine the radial velocity and the [Fe/H] measurements obtained with FLAMES/GIRAFFE at the VLT with a spectral resolution of R = 20 000, plus for the Baade's Window field the OGLE-II proper motions, and compare these with published N-body simulations of the Galactic bulge. Results. We confirm the presence of two distinct populations in Baade's Window found in Hill et al. (2010, A&A, submitted): the metal-rich population presents bar-like kinematics while the metal-poor population shows kinematics corresponding to an old spheroid or a thick disc. In this context the metallicity gradient along the bulge minor axis observed by Zoccali et al. (2008, A&A, 486, 177), visible also in the kinematics, can be related to a varying mix of these two populations as one moves away from the Galactic plane, alleviating the apparent contradiction between the kinematic evidence of a bar and the existence of a metallicity gradient. Conclusions. We show evidence that the two main scenarios for the bulge formation co-exist within the Milky Way bulge.

Social Networks Shape the Transmission Dynamics of Hepatitis C Virus

Hepatitis C virus (HCV) infects 170 million people worldwide, and is a major public health problem in Brazil, where over 1% of the population may be infected and where multiple viral genotypes co-circulate. Chronically infected individuals are both the source of transmission to others and are at risk for HCV-related diseases, such as liver cancer and cirrhosis. Before the adoption of anti-HCV control measures in blood banks, this virus was mainly transmitted via blood transfusion. Today, needle sharing among injecting drug users is the most common form of HCV transmission. Of particular importance is that HCV prevalence is growing in non-risk groups. Since there is no vaccine against HCV, it is important to determine the factors that control viral transmission in order to develop more efficient control measures. However, despite the health costs associated with HCV, the factors that determine the spread of virus at the epidemiological scale are often poorly understood. Here, we sequenced partial NS5b gene sequences sampled from blood samples collected from 591 patients in Sao Paulo state, Brazil. We show that different viral genotypes entered Sao Paulo at different times, grew at different rates, and are associated with different age groups and risk behaviors. In particular, subtype 1b is older and grew more slowly than subtypes 1a and 3a, and is associated with multiple age classes. In contrast, subtypes 1a and 3b are associated with younger people infected more recently, possibly with higher rates of sexual transmission. The transmission dynamics of HCV in Sao Paulo therefore vary by subtype and are determined by a combination of age, risk exposure and underlying social network. We conclude that social factors may play a key role in determining the rate and pattern of HCV spread, and should influence future intervention policies.

Particle competition for complex network community detection

In many real situations, randomness is considered to be uncertainty or even confusion which impedes human beings from making a correct decision. Here we study the combined role of randomness and determinism in particle dynamics for complex network community detection. In the proposed model, particles walk in the network and compete with each other in such a way that each of them tries to possess as many nodes as possible. Moreover, we introduce a rule to adjust the level of randomness of particle walking in the network, and we have found that a portion of randomness can largely improve the community detection rate. Computer simulations show that the model has good community detection performance and at the same time presents low computational complexity. (C) 2008 American Institute of Physics.

Multistability and Self-Similarity in the Parameter-Space of a Vibro-Impact System

The dynamics of a dissipative vibro-impact system called impact-pair is investigated. This system is similar to Fermi-Ulam accelerator model and consists of an oscillating one-dimensional box containing a point mass moving freely between successive inelastic collisions with the rigid walls of the box. In our numerical simulations, we observed multistable regimes, for which the corresponding basins of attraction present a quite complicated structure with smooth boundary. In addition, we characterize the system in a two-dimensional parameter space by using the largest Lyapunov exponents, identifying self-similar periodic sets. Copyright (C) 2009 Silvio L.T. de Souza et al.

Two-dimensional supersonic nonlinear Schrodinger flow past an extended obstacle

Supersonic flow of a superfluid past a slender impenetrable macroscopic obstacle is studied in the framework of the two-dimensional (2D) defocusing nonlinear Schroumldinger (NLS) equation. This problem is of fundamental importance as a dispersive analog of the corresponding classical gas-dynamics problem. Assuming the oncoming flow speed is sufficiently high, we asymptotically reduce the original boundary-value problem for a steady flow past a slender body to the one-dimensional dispersive piston problem described by the nonstationary NLS equation, in which the role of time is played by the stretched x coordinate and the piston motion curve is defined by the spatial body profile. Two steady oblique spatial dispersive shock waves (DSWs) spreading from the pointed ends of the body are generated in both half planes. These are described analytically by constructing appropriate exact solutions of the Whitham modulation equations for the front DSW and by using a generalized Bohr-Sommerfeld quantization rule for the oblique dark soliton fan in the rear DSW. We propose an extension of the traditional modulation description of DSWs to include the linear ""ship-wave"" pattern forming outside the nonlinear modulation region of the front DSW. Our analytic results are supported by direct 2D unsteady numerical simulations and are relevant to recent experiments on Bose-Einstein condensates freely expanding past obstacles.