Assessment of a high-order finite difference upwind scheme for the simulation of convection-diffusion problems

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

Object oriented modeling of piston engines

This work reports a conception phase of a piston engine global model. The model objective is forecast the motor performance (power, torque and specific consumption as a function of rotation and environmental conditions). Global model or Zero-dimensional is based on flux balance through each engine component. The resulting differential equations represents a compressive unsteady flow, in which, all dimensional variables are areas or volumes. A review is presented first. The ordinary differential equation system is presented and a Runge-Kutta method is proposed to solve it numerically. The model includes the momentum conservation equation to link the gas dynamics with the engine moving parts rigid body mechanics. As an oriented to objects model the documentation follows the UML standard. A discussion about the class diagrams is presented, relating the classes with physical model related. The OOP approach allows evolution from simple models to most complex ones without total code rewrite. Copyright © 2001 Society of Automotive Engineers, Inc.

Array of Bose-Einstein condensates under time-periodic Feshbach-resonance management

The investigation of the dynamics of a discrete soliton in an array of Bose-Einstein condensates under the action of a periodically time-modulated atomic scattering length [Feshbach-resonance management (FRM)] was discussed. The slow and rapid modulations, in comparison with the tunneling frequency were considered. An averaged equation, which was a generalized discrete nonlinear Schrödinger equation, including higher-order effective nonlinearities and intersite nonlinear interactions was derived in the case of the rapid modulation. It was demonstrated that the modulations of sufficient strength results in splitting of the soliton by direct simulations.

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.

Gravitational Capture of Asteroids by Gas Drag

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

Nonlinear Schrodinger equation for a superfluid Fermi gas in the BCS-BEC crossover

We introduce a quasianalytic nonlinear Schrodinger equation with beyond mean-field corrections to describe the dynamics of a zero-temperature dilute superfluid Fermi gas in the crossover from the weak-coupling Bardeen-Cooper-Schrieffer (BCS) regime, where k(F)parallel to a parallel to << 1 with a the s-wave scattering length and k(F) the Fermi momentum, through the unitarity limit k(F)a ->+/-infinity to the Bose-Einstein condensation (BEC) regime where k(F)a > 0. The energy of our model is parametrized using the known asymptotic behavior in the BCS, BEC, and the unitarity limits and is in excellent agreement with accurate Green's-function Monte Carlo calculations. The model generates good results for frequencies of collective breathing oscillations of a trapped Fermi superfluid.

Nonlinear Schrodinger equation for a superfluid Bose gas from weak coupling to unitarity: Study of vortices

We introduce a nonlinear Schrodinger equation to describe the dynamics of a superfluid Bose gas in the crossover from the weak-coupling regime, where an(1/3)<<1 with a the interatomic s-wave scattering length and n the bosonic density, to the unitarity limit, where a ->+infinity. We call this equation the unitarity Schrodinger equation (USE). The zero-temperature bulk equation of state of this USE is parametrized by the Lee-Yang-Huang low-density expansion and Jastrow calculations at unitarity. With the help of the USE we study the profiles of quantized vortices and vortex-core radius in a uniform Bose gas. We also consider quantized vortices in a Bose gas under cylindrically symmetric harmonic confinement and study their profile and chemical potential using the USE and compare the results with those obtained from the Gross-Pitaevskii-type equations valid in the weak-coupling limit. Finally, the USE is applied to calculate the breathing modes of the confined Bose gas as a function of the scattering length.

Fermi acceleration and its suppression in a time-dependent Lorentz gas

Some dynamical properties for a Lorentz gas were studied considering both static and time-dependent boundaries. For the static case, it was confirmed that the system has a chaotic component characterized with a positive Lyapunov exponent. For the time-dependent perturbation, the model was described using a four-dimensional nonlinear map. The behaviour of the average velocity is considered in two different situations: (i) non-dissipative and (ii) dissipative dynamics. Our results confirm that unlimited energy growth is observed for the non-dissipative case. However, and totally new for this model, when dissipation via inelastic collisions is introduced, the scenario changes and the unlimited energy growth is suppressed, thus leading to a phase transition from unlimited to limited energy growth. The behaviour of the average velocity is described using scaling arguments. (C) 2010 Elsevier B.V. All rights reserved.

In-flight and collisional dissipation as a mechanism to suppress Fermi acceleration in a breathing Lorentz gas

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

Temporal dynamics of stomatal conductance of plants under water deficit: can homeostasis be improved by more complex dynamics?

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

Fermentation and epiphytic microflora dynamics in sugar cane silage

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

Fluid dynamics of air in a packed bed: velocity profiles and the continuum model assumption

Air flow through packed beds was analyzed experimentally under conditions ranging from those that reinforce the effect of the wall on the void fraction to those that minimize it. The packing was spherical particles, with a tube-to-particle diameter ratio (D/dp) between 3 and 60. Air flow rates were maintained between 1.3 and 4.44 m3/min, and gas velocity was measured with a Pitot tube positioned above the bed exit. Measurements were made at various radial and angular coordinate values, allowing the distribution of air flow across the bed to be described in detail. Comparison of the experimentally observed radial profiles with those derived from published equations revealed that at high D/dp ratios the measured and calculated velocity profiles behaved similarly. At low ratios, oscillations in the velocity profiles agreed with those in the voidage profiles, signifying that treating the porous medium as a continuum medium is questionable in these cases.

Lattice dynamics of metallic glass Ca70Mg30 on the model of Bhatia and Singh

The lattice dynamical studies of the metallic glass Ca70Mg30 by Bhatia and Singh on their model contained two shortcomings, firstly the electron-ion interaction matrix was wrong and secondly, the numerical value of the bulk modulus of the electron gas was accepted arbitrarily. By modifying the electron-ion dynamical matrix and determining all the model parameters from the experimental data, we made a fresh study of the lattice dynamics of Ca70Mg30 and compared it to the earlier studies of Bhatia and Singh and also with experimental phonons.

Chaos in collapsing Bose-condensed gas

We reinvestigate the dynamics of the grow and collapse of Bose-Einstein condensates in a system of trapped ultracold atoms with negative scattering lengths, and found a new behavior in the long time scale evolution: the number of atoms can go far beyond the static stability limit. The condensed state is described by the solution of the time-dependent nonlinear Schrödinger equation, in a model that includes atomic feeding and three-body dissipation. Our results for the model show that, by changing the feeding parameter and when a substantial depletion of the ground-state exists, a chaotic behavior is found. We consider a criterion proposed by Deissler and Kaneko [Phys. Lett. A 119, 397 (1987)] to diagnose spatiotemporal chaos. ©2000 The American Physical Society.

Dynamics of a collapsing and exploding Bose-Einstein condensed vortex state

The dynamics of small repulsive Bose-Einstein condensed vortex states of 85Rb atoms in a cylindrical traps with low angular momentum was studied. The time-dependent mean-field Gross-Pitaevskii equation was used for the study. The condensates collapsed and atoms ejected via explosion and a remnant condensate with a smaller number of atoms emerges that survived for a long time.