Application of a bounded upwinding scheme to complex fluid dynamics problems

This paper is concerned with an overview of upwinding schemes, and further nonlinear applications of a recently introduced high resolution upwind differencing scheme, namely the ADBQUICKEST [V.G. Ferreira, F.A. Kurokawa, R.A.B. Queiroz, M.K. Kaibara, C.M. Oishi, J.A.Cuminato, A.F. Castelo, M.F. Tomé, S. McKee, assessment of a high-order finite difference upwind scheme for the simulation of convection-diffusion problems, International Journal for Numerical Methods in Fluids 60 (2009) 1-26]. The ADBQUICKEST scheme is a new TVD version of the QUICKEST [B.P. Leonard, A stable and accurate convective modeling procedure based on quadratic upstream interpolation, Computer Methods in Applied Mechanics and Engineering 19 (1979) 59-98] for solving nonlinear balance laws. The scheme is based on the concept of NV and TVD formalisms and satisfies a convective boundedness criterion. The accuracy of the scheme is compared with other popularly used convective upwinding schemes (see, for example, Roe (1985) [19], Van Leer (1974) [18] and Arora & Roe (1997) [17]) for solving nonlinear conservation laws (for example, Buckley-Leverett, shallow water and Euler equations). The ADBQUICKEST scheme is then used to solve six types of fluid flow problems of increasing complexity: namely, 2D aerosol filtration by fibrous filters; axisymmetric flow in a tubular membrane; 2D two-phase flow in a fluidized bed; 2D compressible Orszag-Tang MHD vortex; axisymmetric jet onto a flat surface at low Reynolds number and full 3D incompressible flows involving moving free surfaces. The numerical simulations indicate that this convective upwinding scheme is a good generic alternative for solving complex fluid dynamics problems. © 2012.

A comparison of two nonlinear damping mechanisms in a vibration isolator

The vibration transmissibility characteristics of a single-degree-of- freedom (SDOF) passive vibration isolation system with different nonlinear dampers are investigated in this paper. In one configuration, the damper is assumed to be linear and viscous, and is connected to the mass so that it is perpendicular to the spring (horizontal damper). The vibration is in the direction of the spring. The second configuration is one in which the damper is in parallel with the spring but the damping force is proportional to the cube of the relative velocity across the damper (cubic damping). Both configurations are studied for small amplitudes of excitation, when some analysis can be conducted based on analytical expressions, and for large amplitudes of excitation, where the analysis is based on numerical simulations. It is found that the two nonlinear systems can outperform the linear system when force transmissibility is considered. However, for displacement transmissibility, the system with the horizontal damper exhibits some desirable properties, but the system with cubic damping does not. © 2012 Elsevier Ltd.

A compound model for the origin of Earth's water

One of the most important subjects of debate in the formation of the solar system is the origin of Earth's water. Comets have long been considered as the most likely source of the delivery of water to Earth. However, elemental and isotopic arguments suggest a very small contribution from these objects. Other sources have also been proposed, among which local adsorption of water vapor onto dust grains in the primordial nebula and delivery through planetesimals and planetary embryos have become more prominent. However, no sole source of water provides a satisfactory explanation for Earth's water as a whole. In view of that, using numerical simulations, we have developed a compound model incorporating both the principal endogenous and exogenous theories, and investigating their implications for terrestrial planet formation and water delivery. Comets are also considered in the final analysis, as it is likely that at least some of Earth's water has cometary origin. We analyze our results comparing two different water distribution models, and complement our study using the D/H ratio, finding possible relative contributions from each source and focusing on planets formed in the habitable zone. We find that the compound model plays an important role by showing greater advantage in the amount and time of water delivery in Earth-like planets. © 2013. The American Astronomical Society. All rights reserved.

Study of the ℤ3 symmetry in QCD at finite temperature and chemical potential using a worm algorithm

Traditional Monte Carlo simulations of QCD in the presence of a baryon chemical potential are plagued by the complex phase problem and new numerical approaches are necessary for studying the phase diagram of the theory. In this work we consider a ℤ3 Polyakov loop model for the deconfining phase transition in QCD and discuss how a flux representation of the model in terms of dimer and monomer variable solves the complex action problem. We present results of numerical simulations using a worm algorithm for the specific heat and two-point correlation function of Polyakov loops. Evidences of a first order deconfinement phase transition are discussed. © 2013 American Institute of Physics.

An analysis of the Hygiea asteroid family orbital region

(10) Hygiea is the fourth largest asteroid of the main belt, by volume and mass, and it is the largest member of its family, that is made mostly by low-albedo, C-type asteroids, typical of the outer main belt. Like many other large families, it is associated with a 'halo' of objects, that extends far beyond the boundary of the core family, as detected by traditional hierarchical clustering methods (HCM) in proper element domains. Numerical simulations of the orbital evolution of family members may help in estimating the family and halo family age, and the original ejection velocity field. But, in order to minimize the errors associated with including too many interlopers, it is important to have good estimates of family membership that include available data on local asteroid taxonomy, geometrical albedo and local dynamics. For this purpose, we obtained synthetic proper elements and frequencies of asteroids in the Hygiea orbital region, with their errors. We revised the current knowledge on asteroid taxonomy, including Sloan Digital Sky Survey-Moving Object Catalog 4th release (SDSS-MOC 4) data, and geometric albedo data from Wide-field Infrared Survey Explorer (WISE) and Near-Earth Object WISE (NEOWISE). We identified asteroid family members using HCM in the domain of proper elements (a, e, sin (i)) and in the domains of proper frequencies most appropriate to study diffusion in the local web of secular resonances, and eliminated possible interlopers based on taxonomic and geometrical albedo considerations. To identify the family halo, we devised a new hierarchical clustering method in an extended domain that includes proper elements, principal components PC1, PC2 obtained based on SDSS photometric data and, for the first time, WISE and NEOWISE geometric albedo. Data on asteroid size distribution, light curves and rotations were also revised for the Hygiea family. The Hygiea family is the largest group in its region, with two smaller families in proper element domain and 18 families in various frequencies domains identified in this work for the first time. Frequency groups tend to extend vertically in the (a, sin (i)) plane and cross not only the Hygiea family but also the near C-type families of Themis and Veritas, causing a mixture of objects all of relatively low albedo in the Hygiea family area. A few high-albedo asteroids, most likely associated with the Eos family, are also present in the region. Finally, the new multidomains hierarchical clustering method allowed us to obtain a good and robust estimate of the membership of the Hygiea family halo, quite separated from other asteroids families halo in the region, and with a very limited (about 3 per cent) presence of likely interlopers. © 2013 The Author Published by Oxford University Press on behalf of the Royal Astronomical Society.

Small particles in Pluto's environment: Effects of the solar radiation pressure

Impacts of micrometeoroids on the surfaces of the plutonian small satellites Nix and Hydra can generate dust particles. Even in this region so far from the Sun these tiny ejected particles are under the effects of the solar radiation pressure. In this work, we investigate the orbital evolution of the escaping ejecta from both the small satellites under the effects of the radiation pressure combined with the gravitational effects of Pluto, Charon, Nix and Hydra. The mass production rate of micron-sized dust particles generated by micrometeoroids hitting the satellites is obtained, and numerical simulations are performed to derive the lifetime of the ejecta. These pieces of information allow us to estimate the optical depth of a putative ring, which extends from the orbits of Nix to Hydra. The ejected particles, between the orbits of Nix and Hydra, form a wide ring of about 16 000 km. Collisions with the massive bodies and escape from the system are mainly determined by the effects of the solar radiation pressure. This is an important loss mechanism, removing 30 per cent of the initial set of 1 μm-sized particles in 1 yr. The surviving particles form a ring too faint to be detectable with the derived maximum optical depth of 4 × 10-11. © 2013 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society.

Dynamics characterization of modified Gross-Pitaevskii equation

The dynamics of dissipative and coherent N-body systems, such as a Bose-Einstein condensate, which can be described by an extended Gross-Pitaevskii formalism, is investigated. In order to analyze chaotic and unstable regimes, two approaches are considered: a metric one, based on calculations of Lyapunov exponents, and an algorithmic one, based on the Lempel-Ziv criterion. The consistency of both approaches is established, with the Lempel-Ziv algorithmic found as an efficient complementary approach to the metric one for the fast characterization of dynamical behaviors obtained from finite sequences. © 2013 Elsevier B.V. All rights reserved.

On a model for the growth of an invasive avascular tumor

This paper is devoted to study the 1D model of invasive avascular tumor growth, which takes into account cell division, death, and motility, proposed by Kolobov and collaborators in 2009. First, we examine the existence and uniqueness of the solution to this model. Second, we studied qualitatively and numerically the traveling wave solutions. Finally, we show some numerical simulations for the cell density and nutrient concentration. © 2013 NSP Natural Sciences Publishing Cor.

Otimização vetorial e técnicas de mistura de herbicidas aplicadas ao controle de plantas daninhas

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

Regiões de estabilidade no sistema Plutão-Caronte

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

Análise teórica e experimental de um método interferométrico de detecção de fase óptica, auto-consistente e com elevada faixa dinâmica, aplicado à caracterização de atuadores piezoelétricos flextensionais

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

Fenômenos não-lineares, incluindo-se os não-ideias, em captura de energia utilizando-se dispositivos piezoelétricos

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

Comportamento dinâmico não linear e controle de sistemas eletromecânicos em macro e micro escalas

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

Modelagem matemática em câncer: dinâmica angiogênica e quimioterapia anti-neoplásica

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

Modelagem matemática da doença do caranguejo letárgico via ondas viajantes

Pós-graduação em Biometria - IBB