MIXMASTER NUMERICAL BEHAVIOR AND GENERALIZATIONS

The numerical simulation of the mixmaster universe serves the purpose of suggesting two kinds of results. The intrinsic time evolution, during contraction, will be seen to be nonchaotic. This is a necessary feature of relativistic cosmological models undergoing this kind of motion. The mixmaster model also provides a clue on how to define chaoticity for systems described by nonautonomous sets of differential equations.

Spatial distribution of solutes and water in sucrose solution dehydrated apples

Half-fresh apples were immersed in sucrose solution (50% w/w, 27 degrees C) during different times of exposition (2, 4, and 8 h). Then each fruit was sliced from the transversal exposed surface. Density, water, and sugar content were determined for each slice. A mathematical model was fitted to experimental data of water and sucrose content considering the global flux and the tissue shrinkage. By numerical analysis, the binary effective diffusion coefficients as a function of concentration were calculated, using material coordinates and integrating simultaneously two differential equations (for water and sucrose). The coefficients obtained are one or even two orders of magnitude lower than the ones for pure solutions and present an unusual concentration dependence. This comparison shows the influence of the tissue resistance to the diffusion.

Centers on center manifolds in the Lu system

We confirm a conjecture of Mello and Coelho [Phys. Lett. A 373 (2009) 1116] concerning the existence of centers on local center manifolds at equilibria of the Lu system of differential equations on R(3). Our proof shows that the local center manifolds are algebraic ruled surfaces, and are unique. (C) 2011 Elsevier B.V. All rights reserved.

GENERALIZED HEAVISIDE FUNCTIONS IN THE COLOMBEAU THEORY CONTEXT

We defined generalized Heaviside functions for a variable x in R-n, and for variables (x, t) in R-n x R-m. Then study properties such as: composition, invertibility, and association relation (the weak equality). This work is developed in the Colombeau generalized functions context.

Two-body Dirac equation approach to the deuteron

The two-body Dirac(Breit) equation with potentials associated to one-boson-exchanges with cutoff masses is solved for the deuteron and its observables calculated. The 16-component wave-function for the Jπ = 1+ state contains four independent radial functions which satisfy a system of four coupled differential equations of first order. This system is numerically integrated, from infinity towards the origin, by fixing the value of the deuteron binding energy and imposing appropriate boundary conditions at infinity. For the exchange potential of the pion, a mixture of direct plus derivative couplings to the nucleon is considered. We varied the pion-nucleon coupling constant, and the best results of our calculations agree with the lower values recently determined for this constant.

On the integrable perturbations of the Camassa-Holm equation

We present an investigation of the nonlinear partial differential equations (PDE) which are asymptotically representable as a linear combination of the equations from the Camassa-Holm hierarchy. For this purpose we use the infinitesimal transformations of dependent and independent variables of the original PDE. This approach is helpful for the analysis of the systems of the PDE which can be asymptotically represented as the evolution equations of polynomial structure. © 2000 American Institute of Physics.

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.

A two phase flow nonlinear parameter estimation for capillary tube-suction line heat exchanger

Here we present two-phase flow nonlinear parameter estimation for HFC's flow through capillary tube-suction line heat exchangers, commonly used as expansion devices in small refrigeration systems. The simplifying assumptions adopted are: steady state, pure refrigerant, one-dimensional flow, negligible axial heat conduction in the fluid, capillary tube and suction line walls. Additionally, it is considered that the refrigerant is free from oil and both phases are assumed to be at the same pressure, that is, surface tension effects are neglected. Metastable flow effects are also disregarded, and the vapor is assumed to be saturated at the local pressure. The so-called homogeneous model, involving three, first order, ordinary differential equations is applied to analyze the two-phase flow region. Comparison is done with experimental measurements of the mass flow rate and temperature distribution along capillary tubes working with refrigerant HFC-134a in different operating conditions.

On non-ideal and non-linear portal frame dynamics analysis using Bogoliubov averaging method

We apply the Bogoliubov Averaging Method to the study of the vibrations of an elastic foundation, forced by a Non-ideal energy source. The considered model consists of a portal plane frame with quadratic nonlinearities, with internal resonance 1:2, supporting a direct current motor with limited power. The non-ideal excitation is in primary resonance in the order of one-half with the second mode frequency. The results of the averaging method, plotted in time evolution curve and phase diagrams are compared to those obtained by numerically integrating of the original differential equations. The presence of the saturation phenomenon is verified by analytical procedures.

Controlling collapse in Bose-Einstein condensates by temporal modulation of the scattering length

A study was conducted on the dynamics of 2D and 3D Bose-Einstein condensates in the case when the scattering length in the Gross-Pitaevskii (GP) equation which contains constant (dc) and time-variable (ac) parts. Using the variational approximation (VA), simulating the GP equation directly, and applying the averaging procedure to the GP equation without the use of the VA, it was demonstrated that the ac component of the nonlinearity makes it possible to maintain the condensate in a stable self-confined state without external traps.

Stable two-dimensional dispersion-managed soliton

The existence of a dispersion-managed soliton in two-dimensional nonlinear Schrodinger equation with periodically varying dispersion has been explored. The averaged equations for the soliton width and chirp are obtained which successfully describe the long time evolution of the soliton. The slow dynamics of the soliton around the fixed points for the width and chirp are investigated and the corresponding frequencies are calculated. Analytical predictions are confirmed by direct partial differential equation (PDE) and ordinary differential equation (ODE) simulations. Application to a Bose-Einstein condensate in optical lattice is discussed. The existence of a dispersion-managed matter-wave soliton in such system is shown.

Analytical/numerical hybrid solution for one-dimensional ablation problem

Ablation is a thermal protection process with several applications in engineering, mainly in the field of airspace industry. The use of conventional materials must be quite restricted, because they would suffer catastrophic flaws due to thermal degradation of their structures. However, the same materials can be quite suitable once being protected by well-known ablative materials. The process that involves the ablative phenomena is complex, could involve the whole or partial loss of material that is sacrificed for absorption of energy. The analysis of the ablative process in a blunt body with revolution geometry will be made on the stagnation point area that can be simplified as a one-dimensional plane plate problem, hi this work the Generalized Integral Transform Technique (GITT) is employed for the solution of the non-linear system of coupled partial differential equations that model the phenomena. The solution of the problem is obtained by transforming the non-linear partial differential equation system to a system of coupled first order ordinary differential equations and then solving it by using well-established numerical routines. The results of interest such as the temperature field, the depth and the rate of removal of the ablative material are presented and compared with those ones available in the open literature.

Modal Analysis of a Beam with a Tip Rotor by Using a Fundamental Response

The shape modes of a damped-free beam model with a tip rotor are determined by using a dynamical basis that is generated by a fundamental spatial free response. This is a non-classical distributed model for the displacements in the transverse directions of the beam which turns out to be coupled through boundary conditions due to rotation. Numerical calculations are performed by using the Ritz-Rayleigh method with several approximating basis.

Behavioural Analysis of a Nonlinear Mechanical System Using Transient Trajectories

This work aims at a better comprehension of the features of the solution surface of a dynamical system presenting a numerical procedure based on transient trajectories. For a given set of initial conditions an analysis is made, similar to that of a return map, looking for the new configuration of this set in the first Poincaré sections. The mentioned set of I.C. will result in a curve that can be fitted by a polynomial, i.e. an analytical expression that will be called initial function in the undamped case and transient function in the damped situation. Thus, it is possible to identify using analytical methods the main stable regions of the phase portrait without a long computational time, making easier a global comprehension of the nonlinear dynamics and the corresponding stability analysis of its solutions. This strategy allows foreseeing the dynamic behavior of the system close to the region of fundamental resonance, providing a better visualization of the structure of its phase portrait. The application chosen to present this methodology is a mechanical pendulum driven through a crankshaft that moves horizontally its suspension point.

Euclidean 4d exact solitons in a Skyrme type model

We introduce a Skyrme type, four-dimensional Euclidean field theory made of a triplet of scalar fields n→, taking values on the sphere S2, and an additional real scalar field φ, which is dynamical only on a three-dimensional surface embedded in R4. Using a special ansatz we reduce the 4d non-linear equations of motion into linear ordinary differential equations, which lead to the construction of an infinite number of exact soliton solutions with vanishing Euclidean action. The theory possesses a mass scale which fixes the size of the solitons in way which differs from Derrick's scaling arguments. The model may be relevant to the study of the low energy limit of pure SU(2) Yang-Mills theory. © 2004 Elsevier B.V. All rights reserved.