3D calculation of Tucson-Melbourne 3NF effect in triton binding energy

As an application of the new realistic three-dimensional (3D) formalism reported recently for three-nucleon (3N) bound states, an attempt is made to study the effect of three-nucleon forces (3NFs) in triton binding energy in a non partial wave (PW) approach. The spin-isospin dependent 3N Faddeev integral equations with the inclusion of 3NFs, which are formulated as function of vector Jacobi momenta, specifically the magnitudes of the momenta and the angle between them, are solved with Bonn-B and Tucson-Melbourne NN and 3N forces in operator forms which can be incorporated in our 3D formalism. The comparison with numerical results in both, novel 3D and standard PW schemes, shows that non PW calculations avoid the very involved angular momentum algebra occurring for the permutations and transformations and it is more efficient and less cumbersome for considering the 3NF.

Superfluid Bose-Fermi mixture from weak coupling to unitarity

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

Permutability of Backlund transformations for N=2 supersymmetric sine-Gordon

The permutability of two Backlund transformations is employed to construct a nonlinear superposition formula and to generate a class of solutions for the N=2 super sine-Gordon model. We present explicitly the one and two soliton solutions.

Toward the Application of Three-Dimensional Approach to Few-body Atomic Bound States

The first step toward the application of an effective non partial wave (PW) numerical approach to few-body atomic bound states has been taken. The two-body transition amplitude which appears in the kernel of three-dimensional Faddeev-Yakubovsky integral equations is calculated as function of two-body Jacobi momentum vectors, i.e. as a function of the magnitude of initial and final momentum vectors and the angle between them. For numerical calculation the realistic interatomic interactions HFDHE2, HFD-B, LM2M2 and TTY are used. The angular and momentum dependence of the fully off-shell transition amplitude is studied at negative energies. It has been numerically shown that, similar to the nuclear case, the transition amplitude exhibits a characteristic angular behavior in the vicinity of He-4 dimer pole.

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 implicit BEM formulation to model strong discontinuities in solids

A boundary element method (BEM) formulation to predict the behavior of solids exhibiting displacement (strong) discontinuity is presented. In this formulation, the effects of the displacement jump of a discontinuity interface embedded in an internal cell are reproduced by an equivalent strain field over the cell. To compute the stresses, this equivalent strain field is assumed as the inelastic part of the total strain. As a consequence, the non-linear BEM integral equations that result from the proposed approach are similar to those of the implicit BEM based on initial strains. Since discontinuity interfaces can be introduced inside the cell independently on the cell boundaries, the proposed BEM formulation, combined with a tracking scheme to trace the discontinuity path during the analysis, allows for arbitrary discontinuity propagation using a fixed mesh. A simple technique to track the crack path is outlined. This technique is based on the construction of a polygonal line formed by segments inside the cells, in which the assumed failure criterion is reached. Two experimental concrete fracture tests were analyzed to assess the performance of the proposed formulation.

MODEL INDEPENDENCE OF SCATTERING OF 3 IDENTICAL BOSONS IN 2 DIMENSIONS

Within the framework of scattering integral equations in momentum space, we present numerical results of scattering of three identical bosons at low energies in two dimensions for short-range separable potentials. An analysis of the present numerical results reveals the three-particle scattering observables to be independent of potential shape provided the low-energy two-particle binding energy and scattering length are held fixed throughout the investigation. We think that the present conclusion of model independence will be valid for any potential, local or nonlocal, whose range is much smaller than the size of the two-particle bound state.

Linear analysis of building floor structures by a BEM formulation based on Reissner's theory

In this work, the plate bending formulation of the boundary element method (BEM) based on the Reissner's hypothesis is extended to the analysis of zoned plates in order to model a building floor structure. In the proposed formulation each sub-region defines a beam or a slab and depending on the way the sub-regions are represented, one can have two different types of analysis. In the simple bending problem all sub-regions are defined by their middle surface. on the other hand, for the coupled stretching-bending problem all sub-regions are referred to a chosen reference surface, therefore eccentricity effects are taken into account. Equilibrium and compatibility conditions are automatically imposed by the integral equations, which treat this composed structure as a single body. The bending and stretching values defined on the interfaces are approximated along the beam width, reducing therefore the number of degrees of freedom. Then, in the proposed model the set of equations is written in terms of the problem values on the beam axis and on the external boundary without beams. Finally some numerical examples are presented to show the accuracy of the proposed model.

A BEM formulation based on Reissner's theory to perform simple bending analysis of plates reinforced by rectangular beams

In this work, the plate bending formulation of the boundary element method (BEM), based on the Reissner's hypothesis, is extended to the analysis of plates reinforced by rectangular beams. This composed structure is modelled by a zoned plate, being the beams represented by narrow sub-regions with larger thickness. The integral equations are derived by applying the weighted residual method to each sub-region, and summing them to get the equation for the whole plate. Equilibrium and compatibility conditions are automatically imposed by the integral equations, which treat this composed structure as a single body. In order to decrease the number of degrees of freedom, some approximations are considered for both displacements and tractions along the beam width. The accuracy of the proposed model is illustrated by simple examples whose exact solution are known as well as by more complex examples whose numerical results are compared with a well-known finite element code.

Wide localized solitons in systems with time- and space-modulated nonlinearities

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

Theoretical Dynamics of Experimental Populations of Introduced and Native Blowflies (Diptera: Calliphoridae)

Equilibrium dynamics in experimental populations of Chrysomya megacephala (F.) and C. putoria (Wiedemann), which have recently invaded the Americas, and the native species Cochliomyia macellaria (F.), were investigated using nonlinear difference equations. A theoretical analysis of the mathematical model using bifurcation theory established the combination of demographic parameters responsible for producing shifts in blowfly population dynamics from stable equilibria to bounded cycles and aperiodic behavior. Mathematical modeling shows that the populations of the 2 introduced Chrysomya species will form stable oscillations with numbers fluctuating 3-4 times in successive generations. However, in the native species C. macellaria, the dynamics is characterized by damping oscillations in population size, leading to a stable population level.

Approximate solution for non-Darcy transient film condensation in a porous medium

The problem of non-darcian transient film condensation adjacent to a vertical flat plate embedded in a porous medium has been considered. The governing equation for the boundary layer thickness was obtained by an integral method and solved approximately by the method of integral relations. It is shown that the results are in good agreement with those obtained exactly by the method of characteristics.

Scaling limit of virtual states of triatomic systems

Natural scales determine the physics of quantum few-body systems with short-range interactions. Thus, the scaling limit is found when the ratio between the scattering length and the interaction range tends to infinity, while the ratio between the physical scales are kept fixed. From the formal point of view, the relation of the scaling limit and the renormalization aspects of a few-body model with a zero-range interaction, through the derivation of subtracted three-body T-matrix equations that are renormalization-group invariant.

Thermally developing forced convection of non-Newtonian fluids inside elliptical ducts

Laminar-forced convection inside tubes of various cross-section shapes is of interest in the design of a low Reynolds number heat exchanger apparatus. Heat transfer to thermally developing, hydrodynamically developed forced convection inside tubes of simple geometries such as a circular tube, parallel plate, or annular duct has been well studied in the literature and documented in various books, but for elliptical duct there are not much work done. The main assumptions used in this work are a non-Newtonian fluid, laminar flow, constant physical properties, and negligible axial heat diffusion (high Peclet number). Most of the previous research in elliptical ducts deal mainly with aspects of fully developed laminar flow forced convection, such as velocity profile, maximum velocity, pressure drop, and heat transfer quantities. In this work, we examine heat transfer in a hydrodynamically developed, thermally developing laminar forced convection flow of fluid inside an elliptical tube under a second kind of a boundary condition. To solve the thermally developing problem, we use the generalized integral transform technique (GITT), also known as Sturm-Liouville transform. Actually, such an integral transform is a generalization of the finite Fourier transform, where the sine and cosine functions are replaced by more general sets of orthogonal functions. The axes are algebraically transformed from the Cartesian coordinate system to the elliptical coordinate system in order to avoid the irregular shape of the elliptical duct wall. The GITT is then applied to transform and solve the problem and to obtain the once unknown temperature field. Afterward, it is possible to compute and present the quantities of practical interest, such as the bulk fluid temperature, the local Nusselt number, and the average Nusselt number for various cross-section aspect ratios.

Evaluating residues and integrals through negative dimensional integration method (NDIM)

The standard way of evaluating residues and some real integrals through the residue theorem (Cauchy's theorem) is well-known and widely applied in many branches of Physics. Herein we present an alternative technique based on the negative dimensional integration method (NDIM) originally developed to handle Feynman integrals. The advantage of this new technique is that we need only to apply Gaussian integration and solve systems of linear algebraic equations, with no need to determine the poles themselves or their residues, as well as obtaining a whole class of results for differing orders of poles simultaneously.