A marker-and-cell approach to free surface 2-D multiphase flows

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

Formation of bright solitons and soliton trains in a fermion-fermion mixture by modulational instability

We employ a time- dependent mean- field- hydrodynamic model to study the generation of bright solitons in a degenerate fermion - fermion mixture in a cigar- shaped geometry using variational and numerical methods. Due to a strong Pauli- blocking repulsion among identical spin- polarized fermions at short distances there cannot be bright solitons for repulsive interspecies interactions. Employing a linear stability analysis we demonstrate the formation of stable solitons due to modulational instability of a constant-amplitude solution of the model equations for a sufficiently attractive interspecies interaction. We perform a numerical stability analysis of these solitons and also demonstrate the formation of soliton trains by jumping the effective interspecies interaction from repulsive to attractive. These fermionic solitons can be formed and studied in laboratory with present technology.

Tightly bound gap solitons in a Fermi gas

Within the framework of the mean-field hydrodynamic model of a degenerate Fermi gas ( DFG), we study, by means of numerical methods and variational approximation ( VA), the formation of fundamental gap solitons ( FGSs) in a DFG ( or in a BCS superfluid generated by weak interaction between spin- up and spin- down fermions), which is trapped in a periodic optical- lattice ( OL) potential. An effectively one- dimensional ( 1D) con. guration is considered, assuming strong transverse confinement; in parallel, a proper 1D model of the DFG ( which amounts to the known quintic equation for the Tonks- Girardeau gas in the OL) is considered too. The FGSs found in the first two bandgaps of the OL- induced spectrum ( unless they are very close to edges of the gaps) feature a ( tightly bound) shape, being essentially confined to a single cell of the OL. In the second bandgap, we also find antisymmetric tightly bound subfundamental solitons ( SFSs), with zero at the midpoint. The SFSs are also confined to a single cell of the OL, but, unlike the FGSs, they are unstable. The predicted solitons, consisting of similar to 10(4) - 10(5) atoms, can be created by available experimental techniques in the DFG of Li-6 atoms.

Bright solitons and soliton trains in a fermion-fermion mixture

We use a time-dependent dynamical mean-field-hydrodynamic model to predict and study bright solitons in a degenerate fermion-fermion mixture in a quasi-one-dimensional cigar-shaped geometry using variational and numerical methods. Due to a strong Pauli-blocking repulsion among identical spin-polarized fermions at short distances there cannot be bright solitons for repulsive interspecies fermion-fermion interactions. However, stable bright solitons can be formed for a sufficiently attractive interspecies interaction. We perform a numerical stability analysis of these solitons and also demonstrate the formation of soliton trains. These fermionic solitons can be formed and studied in laboratory with present technology.

Cubic interactions and quantum criticality in dimerized antiferromagnets

In certain Mott-insulating dimerized antiferromagnets, triplet excitations of the paramagnetic phase display both three-particle and four-particle interactions. When such a magnet undergoes a quantum phase transition into a magnetically ordered state, the three-particle interaction becomes part of the critical theory provided that the lattice ordering wave vector is zero. One microscopic example is the staggered-dimer antiferromagnet on the square lattice, for which deviations from O(3) universality have been reported in numerical studies. Using both symmetry arguments and microscopic calculations, we show that a nontrivial cubic term arises in the relevant order-parameter quantum field theory, and we assess its consequences using a combination of analytical and numerical methods. We also present finite-temperature quantum Monte Carlo data for the staggered-dimer antiferromagnet which complement recently published results. The data can be consistently interpreted in terms of critical exponents identical to that of the standard O(3) universality class, but with anomalously large corrections to scaling. We argue that the cubic interaction of critical triplons, although irrelevant in two spatial dimensions, is responsible for the leading corrections to scaling due to its small scaling dimension.

High-order upwinding and the hydraulic jump

This work is concerned with the computation of incompressible axisymmetric and fall three-dimensional free-surface flows. In particular, the circular-hydraulic jump is simulated and compared with approximate analytic solutions. However, the principal thrust of this paper is to provide a real problem as a test bed for comparing the many existing convective approximations. Their performance is compared; SMART, HLPA and VONOS emerge as acceptable upwinding methods for this problem. Copyright (C) 2002 John Wiley Sons, Ltd.

Escape and transport for an open bouncer: Stretched exponential decays

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

On non-linear dynamics and an optimal control synthesis of the action potential of membranes (ideal and non-ideal cases) of the Hodgkin-Huxley (HH) mathematical model

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

Modelagem matemática e simulação numérica do resfriamento rápido de morango com ar forçado

Este trabalho abordou o resfriamento rápido com ar forçado de morango via simulação numérica. Para tanto, foi empregado o modelo matemático que descreve o processo de transferência de calor, com base na lei de Fourier, escrito em coordenadas esféricas e simplificado para descrever o processo unidimensional. A resolução da equação expressa pelo modelo matemático deu-se por meio da implementação de um algoritmo, fundamentado no esquema explícito do método numérico das diferenças finitas, executado no ambiente de computação científica MATLAB 6.1. A validação do modelo matemático foi realizada a partir da comparação de dados teóricos com dados obtidos num experimento, no qual morangos foram resfriados com ar forçado. Os resultados mostraram que esse tipo de investigação para a determinação do coeficiente de transferência de calor por convecção é promissora como ferramenta no suporte à decisão do uso ou desenvolvimento de equipamentos na área de resfriamento rápido de frutos esféricos com ar forçado.

Solitons in Bose-Einstein condensates trapped in a double-well potential

We investigate, analytically and numerically, families of bright solitons in a system of two linearly coupled nonlinear Schrodinger/Gross-Pitaevskii equations, describing two Bose-Einstein condensates trapped in an asymmetric double-well potential, in particular, when the scattering lengths in the condensates have arbitrary magnitudes and opposite signs. The solitons are found to exist everywhere where they are permitted by the dispersion law. Using the Vakhitov-Kolokolov criterion and numerical methods, we show that, except for small regions in the parameter space, the solitons are stable to small perturbations. Some of them feature self-trapping of almost all the atoms in the condensate with no atomic interaction or weak repulsion is coupled to the self-attractive condensate. An unusual bifurcation is found, when the soliton bifurcates from the zero solution with vanishing amplitude and width simultaneously diverging but at a finite number of atoms in the soliton. By means of numerical simulations, it is found that, depending on values of the parameters and the initial perturbation, unstable solitons either give rise to breathers or completely break down into incoherent waves (radiation). A version of the model with the self-attraction in both components, which applies to the description of dual-core fibers in nonlinear optics, is considered too, and new results are obtained for this much studied system. (C) 2003 Elsevier B.V. All rights reserved.

An application of the Element Free Galerkin Method to the analysis of quantum well structures

In order to obtain the quantum-mechanical properties of layered semicondutor structures (quantum well and superlattice structures, for instance), solutions of the Schrodinger equation should be obtained for arbitrary potential profiles. In this paper, it is shown that such problems may be also studied by the Element Free Galerkin Method.

Gap solitons in superfluid boson-fermion mixtures

Using coupled equations for the bosonic and fermionic order parameters, we construct families of gap solitons (GSs) in a nearly one-dimensional Bose-Fermi mixture trapped in a periodic optical-lattice (OL) potential, the boson and fermion components being in the states of the Bose-Einstein condensation and Bardeen-Cooper-Schrieffer superfluid, respectively. Fundamental GSs are compact states trapped, essentially, in a single cell of the lattice. Full families of such solutions are constructed in the first two band gaps of the OL-induced spectrum, by means of variational and numerical methods, which are found to be in good agreement. The families include both intragap and intergap solitons, with the chemical potentials of the boson and fermion components falling in the same or different band gaps, respectively. Nonfundamental states, extended over several lattice cells, are constructed too. The GSs are stable against strong perturbations.

On non-linear dynamics and a periodic control design applied to the potential of membrane action

The Fitzhugh-Nagumo (fn) mathematical model characterizes the action potential of the membrane. The dynamics of the Fitzhugh-Nagumo model have been extensively studied both with a view to their biological implications and as a test bed for numerical methods, which can be applied to more complex models. This paper deals with the dynamics in the (FH) model. Here, the dynamics are analyzed, qualitatively, through the stability diagrams to the action potential of the membrane. Furthermore, we also analyze quantitatively the problem through the evaluation of Floquet multipliers. Finally, the nonlinear periodic problem is controlled, based on the Chebyshev polynomial expansion, the Picard iterative method and on Lyapunov-Floquet transformation (L-F transformation).

Charge-orbital ordering and phase separation in the two-orbital model for manganites: Roles of Jahn-Teller phononic and Coulombic interactions

The main properties of realistic models for manganites are studied using analytic mean-field approximations and computational numerical methods, focusing on the two-orbital model with electrons interacting through Jahn-Teller (JT) phonons and/or Coulombic repulsions. Analyzing the model including both interactions by the combination of the mean-field approximation and the exact diagonalization method, it is argued that the spin-charge-orbital structure in the insulating phase of the purely JT-phononic model with a large Hund couphng J(H) is not qualitatively changed by the inclusion of the Coulomb interactions. As an important application of the present mean-held approximation, the CE-type antiferromagnetic state, the charge-stacked structure along the z axis, and (3x(2) - r(2))/(3y(2) - r(2))-type orbital ordering are successfully reproduced based on the JT-phononic model with large JH for the half-doped manganite, in agreement with recent Monte Carlo simulation results. Topological arguments and the relevance of the Heisenberg exchange among localized t(2g) spins explains why the inclusion of the nearest-neighbor Coulomb interaction does not destroy the charge stacking structure. It is also verified that the phase-separation tendency is observed both in purely JT-phononic (large JH) and purely Coulombic models in the vicinity of the hole undoped region, as long as realistic hopping matrices are used. This highlights the qualitative similarities of both approaches and the relevance of mixed-phase tendencies in the context of manganites. In addition, the rich and complex phase diagram of the two-orbital Coulombic model in one dimension is presented. Our results provide robust evidence that Coulombic and JT-phononic approaches to manganites are not qualitatively different ways to carry out theoretical calculations, but they share a variety of common features.

Bayesian outlier analysis in binary regression

We propose alternative approaches to analyze residuals in binary regression models based on random effect components. Our preferred model does not depend upon any tuning parameter, being completely automatic. Although the focus is mainly on accommodation of outliers, the proposed methodology is also able to detect them. Our approach consists of evaluating the posterior distribution of random effects included in the linear predictor. The evaluation of the posterior distributions of interest involves cumbersome integration, which is easily dealt with through stochastic simulation methods. We also discuss different specifications of prior distributions for the random effects. The potential of these strategies is compared in a real data set. The main finding is that the inclusion of extra variability accommodates the outliers, improving the adjustment of the model substantially, besides correctly indicating the possible outliers.