Three-dimensional, fully adaptive simulations of phase-field fluid models

We present an efficient numerical methodology for the 31) computation of incompressible multi-phase flows described by conservative phase-field models We focus here on the case of density matched fluids with different viscosity (Model H) The numerical method employs adaptive mesh refinements (AMR) in concert with an efficient semi-implicit time discretization strategy and a linear, multi-level multigrid to relax high order stability constraints and to capture the flow`s disparate scales at optimal cost. Only five linear solvers are needed per time-step. Moreover, all the adaptive methodology is constructed from scratch to allow a systematic investigation of the key aspects of AMR in a conservative, phase-field setting. We validate the method and demonstrate its capabilities and efficacy with important examples of drop deformation, Kelvin-Helmholtz instability, and flow-induced drop coalescence (C) 2010 Elsevier Inc. All rights reserved

The shoving model for the glass-former LiCl center dot 6H(2)O: A molecular dynamics simulation study

Molecular dynamics (MD) simulations of LiCl center dot 6H(2)O Showed that the diffusion coefficient D, and also I lie structural relaxation time , follow a power law at high temperatures, D(-1) proportional to (T - T(0))(-mu), with the same experimental parameters for viscosity (T(0) = 207 K, mu = 2.08). Decoupling between D and occurs at T(x) similar to 1.1 T(0). High frequency acoustic excitations for the LiCl center dot 6H(2)O model were obtained by the calculation of time correlation functions of mass current fluctuations. The temperature dependence of the instantaneous shear modulus, G,(T), was considered in the shoving model for supercooled liquids [J.C. Dyre, T. Christensen, N.B. Olsen, J. Non-Cryst. Solids 352 (2006) 4635] resulting in a linear relationship log (D(-1)) vs. G root T. (C) 2009 Elsevier B.V. All rights reserved.

Modelling water adsorption on Au(210) surfaces: II. Monte Carlo simulations

Canonical Monte Carlo simulations for the Au(210)/H(2)O interface, using a force field recently proposed by us, are reported. The results exhibit the main features normally observed in simulations of water molecules in contact with different noble metal surfaces. The calculations also assess the influence of the surface topography on the structural aspects of the adsorbed water and on the distribution of the water molecules in the direction normal to the metal surface plane. The adsorption process is preferential at sites in the first layer of the metal. The analysis of the density profiles and dipole moment distributions points to two predominant orientations. Most of the molecules are adsorbed with the molecular plane parallel to surface, while others adsorb with one of the O-H bonds parallel to the surface and the other bond pointing towards the bulk liquid phase. There is also evidence of hydrogen bond formation between the first and second solvent layers at the interface. (c) 2007 Elsevier B.V. All rights reserved.

Coherence dynamics of two-mode condensates in asymmetric potentials

Detection of weak forces with an accuracy beyond the standard quantum limit holds promise both for fundamental research and for technological applications. Schemes involving ultracold atoms for such measurements are now considered to be prime candidates for increased sensitivity. In this paper we use a combination of analytical and numerical techniques to investigate the possible subshot-noise estimation of applied force fields through detection of coherence dynamics of Bose-condensed atoms in asymmetric double-well traps. Following a semiclassical description of the system dynamics and fringe visibility, we present numerical simulations of the full quantum dynamics that demonstrate the dynamical production of phase squeezing beyond the standard quantum limit. Nonlinear interactions are found to limit the achievable amount to a finite value determined by the external weak force.

Dynamics of meso and thermo citrate synthases with implicit solvation

The dynamics of hydration of meso and thermo citrate synthases has been investigated using the EEF1 methodology implemented with the CHARNM program. The native enzymes are composed of two identical subunits, each divided into a small and large domain. The dynamics behavior of both enzymes at 30 degrees C and 60 degrees C has been compared. The results of simulations show that during the hydration process, each subunit follows a different pathway of hydration, in spite of the identical sequence. The hydrated structures were compared with the crystalline structure, and the root mean square deviation (RMSD) of each residue along the trajectory was calculated. The regions with larger and smaller mobility were identified. In particular, helices belonging to the small domain are more mobile than those of the large domain. In contrast, the residues that constitute the active site show a much lower displacement compared with the crystalline structure. Hydration free energy calculations point out that Thermoplasma acidophilum citrate synthase (TCS) is more stable than chicken citrate synthase (CCS), at high temperatures. Such result has been ascribed to the higher number of superficial charges in the thermophilic homologue, which stabilizes the enzyme, while the mesophilic homologue denatures. These results are in accord with the experimental found that TCS keeps activity at temperatures farther apart from the catalysis regular temperature than the CCS. (c) 2005 Wiley Periodicals, Inc.

Population dynamics, life stage and ecological modeling in Chrysomya albiceps (Wiedemann) (Diptera: Calliphoridae)

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Spatio-temporal dynamics and transition from asymptotic equilibrium to bounded oscillations in Chrysomya albiceps (Diptera, Calliphoridae)

The sensitivity of parameters that govern the stability of population size in Chrysomya albiceps and describe its spatial dynamics was evaluated in this study. The dynamics was modeled using a density-dependent model of population growth. Our simulations show that variation in fecundity and mainly in survival has marked effect on the dynamics and indicates the possibility of transitions from one-point equilibrium to bounded oscillations. C. albiceps exhibits a two-point limit cycle, but the introduction of diffusive dispersal induces an evident qualitative shift from two-point limit cycle to a one fixed-point dynamics. Population dynamics of C. albiceps is here compared to dynamics of Cochliomyia macellaria, C. megacephala and C. putoria.

Dynamics of experimental populations of native and introduced blowflies (Diptera: Calliphoridae): mathematical modelling and the transition from asymptotic equilibrium to bounded oscillations

The equilibrium dynamics of native and introduced blowflies is modelled using a density-dependent model of population growth that takes into account important features of the life-history in these flies. A theoretical analysis indicates that the product of maximum fecundity and survival is the primary determinant of the dynamics. Cochliomyia macellaria, a blowfly native to the Americas and the introduced Chrysomya megacephala and Chrysomya putoria, differ in their dynamics in that the first species shows a damping oscillatory behavior leading to a one-point equilibrium, whereas in the last two species population numbers show a two-point limit cycle. Simulations showed that variation in fecundity has a marked effect on the dynamics and indicates the possibility of transitions from one-point equilibrium to bounded oscillations and aperiodic behavior. Variation in survival has much less influence on the dynamics.

Stochastic dynamics of Musca domestica (Diptera: Muscidae) and Chrysomya putoria (Diptera: Calliphoridae)

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

CHARGED PARTICLE DYNAMICS IN TIME DEPENDENT MAGNETIC FIELDS

In this work we study the behavior of charged particles immersed in a peculiar configuration of magnetic fields, which has a main constant field B(0) and a superimposed, transversal perturbation field B(1) sin(omega(p)t), with B(1) << B(0). By taking Cartesian coordinates and placing B(0) along the z axis and B(1) sin (omega(p)t) on the x axis, an analytical solution for y(t) may be obtained by solving an integrodifferential equation. Besides, the solution z(t) also exhibits a very interesting dynamics, and the entire system is conditioned by resonances between the particle orbit frequencies and the frequency of the magnetic transversal perturbation, omega(p). In this work we also discuss numerical simulations for the related particle trajectories, as well as potential applications in the context of separation phenomena.

Dynamics of the deconfinement transition of quarks and gluons in a finite volume

We employ finite elements methods for the approximation of solutions of the Ginzburg-Landau equations describing the deconfinement transition in quantum chromodynamics. These methods seem appropriate for situations where the deconfining transition occurs over a finite volume as in relativistic heavy ion collisions. where in addition expansion of the system and flow of matter are important. Simulation results employing finite elements are presented for a Ginzburg-Landau equation based on a model free energy describing the deconfining transition in pure gauge SU(2) theory. Results for finite and infinite system are compared. (C) 2009 Elsevier B.V. All rights reserved.

Nebular gas drag and co-orbital system dynamics

Aims. We study trajectories of planetesimals whose orbits decay due to gas drag in a primordial solar nebula and are perturbed by the gravity of the secondary body on an eccentric orbit whose mass ratio takes values from mu(2) = 10(-7) to mu(2) = 10(-3) increasing ten times at each step. Each planetesimal ultimately suffers one of the three possible fates: (1) trapping in a mean motion resonance with the secondary body; (2) collision with the secondary body and consequent increase of its mass; or (3) diffusion after crossing the orbit of the secondary body.Methods. We take the Burlirsh-Stoer numerical algorithm in order to integrate the Newtonian equations of the planar, elliptical restricted three-body problem with the secondary body and the planetesimal orbiting the primary. It is assumed that there is no interaction among planetesimals, and also that the gas does not affect the orbit of the secondary body.Results. The results show that the optimal value of the gas drag constant k for the 1: 1 resonance is between 0.9 and 1.25, representing a meter size planetesimal for each AU of orbital radius. In this study, the conditions of the gas drag are such that in theory, L4 no longer exists in the circular case for a critical value of k that defines a limit size of the planetesimal, but for a secondary body with an eccentricity larger than 0.05 when mu(2) = 10(-6), it reappears. The decrease of the cutoff collision radius increase the difusions but does not affect the distribution of trapping. The contribution to the mass accretion of the secondary body is over 40% with a collision radius 0.05R(Hill) and less than 15% with 0.005R(Hill) for mu(2) = 10(-7). The trappings no longer occur when the drag constant k reachs 30. That means that the size limit of planetesimal trapping is 0.2 m per AU of orbital radius. In most cases, this accretion occurs for a weak gas drag and small secondary eccentricity. The diffusions represent most of the simulations showing that gas drag is an efficient process in scattering planetesimals and that the trapping of planetesimals in the 1: 1 resonance is a less probable fate. These results depend on the specific drag force chosen.

Dissipative dynamics of matter-wave solitons in a nonlinear optical lattice

Dynamics and stability of solitons in two-dimensional (2D) Bose-Einstein condensates (BEC), with one-dimensional (1D) conservative plus dissipative nonlinear optical lattices, are investigated. In the case of focusing media (with attractive atomic systems), the collapse of the wave packet is arrested by the dissipative periodic nonlinearity. The adiabatic variation of the background scattering length leads to metastable matter-wave solitons. When the atom feeding mechanism is used, a dissipative soliton can exist in focusing 2D media with 1D periodic nonlinearity. In the defocusing media (repulsive BEC case) with harmonic trap in one direction and nonlinear optical lattice in the other direction, the stable soliton can exist. Variational approach simulations are confirmed by full numerical results for the 2D Gross-Pitaevskii equation.

Nonequilibrium dynamics of quantum fields

The nonequilibrium effective equation of motion for a scalar background field in a thermal bath is studied numerically. This equation emerges from a microscopic quantum field theory derivation and it is suitable to a Langevin simulation on the lattice. Results for both the symmetric and broken phases are presented.

Short-time dynamics of percolation observables

We consider the critical short-time evolution of magnetic and droplet-percolation order parameters for the Ising model in two and three dimensions, through Monte Carlo simulations with the (local) heat-bath method. We find qualitatively different dynamic behaviors for the two types of order parameters. More precisely, we find that the percolation order parameter does not have a power-law behavior as encountered for the magnetization, but develops a scale (related to the relaxation time to equilibrium) in the Monte Carlo time. We argue that this difference is due to the difficulty in forming large clusters at the early stages of the evolution. Our results show that, although the descriptions in terms of magnetic and percolation order parameters may be equivalent in the equilibrium regime, greater care must be taken to interpret percolation observables at short times. In particular, this concerns the attempts to describe the dynamics of the deconfinement phase transition in QCD using cluster observables.