Combining Legendre's Polynomials and Genetic Algorithm in the Solution of Nonlinear Initial-Value Problems

This work develops a method for solving ordinary differential equations, that is, initial-value problems, with solutions approximated by using Legendre's polynomials. An iterative procedure for the adjustment of the polynomial coefficients is developed, based on the genetic algorithm. This procedure is applied to several examples providing comparisons between its results and the best polynomial fitting when numerical solutions by the traditional Runge-Kutta or Adams methods are available. The resulting algorithm provides reliable solutions even if the numerical solutions are not available, that is, when the mass matrix is singular or the equation produces unstable running processes.

Physical approximations for the nonlinear evolution of perturbations in inhomogeneous dark energy scenarios

The abundance and distribution of collapsed objects such as galaxy clusters will become an important tool to investigate the nature of dark energy and dark matter. Number counts of very massive objects are sensitive not only to the equation of state of dark energy, which parametrizes the smooth component of its pressure, but also to the sound speed of dark energy, which determines the amount of pressure in inhomogeneous and collapsed structures. Since the evolution of these structures must be followed well into the nonlinear regime, and a fully relativistic framework for this regime does not exist yet, we compare two approximate schemes: the widely used spherical collapse model and the pseudo-Newtonian approach. We show that both approximation schemes convey identical equations for the density contrast, when the pressure perturbation of dark energy is parametrized in terms of an effective sound speed. We also make a comparison of these approximate approaches to general relativity in the linearized regime, which lends some support to the approximations.

Nonlinear transport and oscillating magnetoresistance in double quantum wells

The nonlinear regime of low-temperature magnetoresistance of double quantum wells in the region of magnetic fields below 1 T is studied both experimentally and theoretically. The observed inversion of the magnetointersubband oscillation peaks with increasing electric current and splitting of these peaks are described by the theory based on the kinetic equation for the isotropic nonequilibrium part of electron distribution function. The inelastic-scattering time of electrons is determined from the current dependence of the inversion field.

Scalar resonances in a unitary pi pi S-wave model for D(+)->pi(+)pi(-)pi(+)

We propose a model for D(+)->pi(+)pi(-)pi(+) decays following experimental results which indicate that the two-pion interaction in the S wave is dominated by the scalar resonances f(0)(600)/sigma and f(0)(980). The weak decay amplitude for D(+)-> R pi(+), where R is a resonance that subsequently decays into pi(+)pi(-), is constructed in a factorization approach. In the S wave, we implement the strong decay R ->pi(+)pi(-) by means of a scalar form factor. This provides a unitary description of the pion-pion interaction in the entire kinematically allowed mass range m(pi pi)(2) from threshold to about 3 GeV(2). In order to reproduce the experimental Dalitz plot for D(+)->pi(+)pi(-)pi(+), we include contributions beyond the S wave. For the P wave, dominated by the rho(770)(0), we use a Breit-Wigner description. Higher waves are accounted for by using the usual isobar prescription for the f(2)(1270) and rho(1450)(0). The major achievement is a good reproduction of the experimental m(pi pi)(2) distribution, and of the partial as well as the total D(+)->pi(+)pi(-)pi(+) branching ratios. Our values are generally smaller than the experimental ones. We discuss this shortcoming and, as a by-product, we predict a value for the poorly known D ->sigma transition form factor at q(2)=m pi(2).

Nonlinear waves in a quark gluon plasma

We study the propagation of perturbations in the quark gluon plasma. This subject has been addressed in other works and in most of the theoretical descriptions of this phenomenon the hydrodynamic equations have been linearized for simplicity. We propose an alternative approach, also based on hydrodynamics but taking into account the nonlinear terms of the equations. We show that these terms may lead to localized waves or even solitons. We use a simple equation of state for the QGP and expand the hydrodynamic equations around equilibrium configurations. The resulting differential equations describe the propagation of perturbations in the energy density. We solve them numerically and find that localized perturbations can propagate for long distances in the plasma. Under certain conditions our solutions mimic the propagation of Korteweg-de Vries solitons.

Nonlinear spectra of ZnO: reverse saturable, two- and three-photon absorption

We present a broadband (460-980 nm) analysis of the nonlinear absorption processes in bulk ZnO, a large-bandgap material with potential blue-to-UV photonic device applications. Using an optical parametric amplifier we generated tunable 1-kHz repetition rate laser pulses and employed the Z-scan technique to investigate the nonlinear absorption spectrum of ZnO. For excitation wavelengths below 500 nm, we observed reverse saturable absorption due to one-photon excitation of the sample, agreeing with rate-equation modeling. Two-and three-photon absorption were observed from 540 to 980 nm. We also determined the spectral regions exhibiting mixture of nonlinear absorption mechanisms, which were confirmed by photoluminescence measurements. (C) 2010 Optical Society of America

Experimental results on the nonlinear H(infinity) control via quasi-LPV representation and game theory for wheeled mobile robots

In this paper, nonlinear dynamic equations of a wheeled mobile robot are described in the state-space form where the parameters are part of the state (angular velocities of the wheels). This representation, known as quasi-linear parameter varying, is useful for control designs based on nonlinear H(infinity) approaches. Two nonlinear H(infinity) controllers that guarantee induced L(2)-norm, between input (disturbances) and output signals, bounded by an attenuation level gamma, are used to control a wheeled mobile robot. These controllers are solved via linear matrix inequalities and algebraic Riccati equation. Experimental results are presented, with a comparative study among these robust control strategies and the standard computed torque, plus proportional-derivative, controller.

An EFG method for the nonlinear analysis of plates undergoing arbitrarily large deformations

The applicability of a meshfree approximation method, namely the EFG method, on fully geometrically exact analysis of plates is investigated. Based on a unified nonlinear theory of plates, which allows for arbitrarily large rotations and displacements, a Galerkin approximation via MLS functions is settled. A hybrid method of analysis is proposed, where the solution is obtained by the independent approximation of the generalized internal displacement fields and the generalized boundary tractions. A consistent linearization procedure is performed, resulting in a semi-definite generalized tangent stiffness matrix which, for hyperelastic materials and conservative loadings, is always symmetric (even for configurations far from the generalized equilibrium trajectory). Besides the total Lagrangian formulation, an updated version is also presented, which enables the treatment of rotations beyond the parameterization limit. An extension of the arc-length method that includes the generalized domain displacement fields, the generalized boundary tractions and the load parameter in the constraint equation of the hyper-ellipsis is proposed to solve the resulting nonlinear problem. Extending the hybrid-displacement formulation, a multi-region decomposition is proposed to handle complex geometries. A criterium for the classification of the equilibrium`s stability, based on the Bordered-Hessian matrix analysis, is suggested. Several numerical examples are presented, illustrating the effectiveness of the method. Differently from the standard finite element methods (FEM), the resulting solutions are (arbitrary) smooth generalized displacement and stress fields. (c) 2007 Elsevier Ltd. All rights reserved.

Galerkin Method and Weighting Functions Applied to Nonlinear H(infinity) Control with Output Feedback

This paper considers two aspects of the nonlinear H(infinity) control problem: the use of weighting functions for performance and robustness improvement, as in the linear case, and the development of a successive Galerkin approximation method for the solution of the Hamilton-Jacobi-Isaacs equation that arises in the output-feedback case. Design of nonlinear H(infinity) controllers obtained by the well-established Taylor approximation and by the proposed Galerkin approximation method applied to a magnetic levitation system are presented for comparison purposes.

On the convergence of successive Galerkin approximation for nonlinear output feedback H(infinity) control

Although the formulation of the nonlinear theory of H(infinity) control has been well developed, solving the Hamilton-Jacobi-Isaacs equation remains a challenge and is the major bottleneck for practical application of the theory. Several numerical methods have been proposed for its solution. In this paper, results on convergence and stability for a successive Galerkin approximation approach for nonlinear H(infinity) control via output feedback are presented. An example is presented illustrating the application of the algorithm.

Resonant Wave Interactions in the Presence of a Diurnally Varying Heat Source

Resonant interactions among equatorial waves in the presence of a diurnally varying heat source are studied in the context of the diabatic version of the equatorial beta-plane primitive equations for a motionless, hydrostatic, horizontally homogeneous and stably stratified background atmosphere. The heat source is assumed to be periodic in time and of small amplitude [i.e., O(epsilon)] and is prescribed to roughly represent the typical heating associated with deep convection in the tropical atmosphere. In this context, using the asymptotic method of multiple time scales, the free linear Rossby, Kelvin, mixed Rossby-gravity, and inertio-gravity waves, as well as their vertical structures, are obtained as leading-order solutions. These waves are shown to interact resonantly in a triad configuration at the O(e) approximation, and the dynamics of these interactions have been studied in the presence of the forcing. It is shown that for the planetary-scale wave resonant triads composed of two first baroclinic equatorially trapped waves and one barotropic Rossby mode, the spectrum of the thermal forcing is such that only one of the triad components is resonant with the heat source. As a result, to illustrate the role of the diurnal forcing in these interactions in a simplified fashion, two kinds of triads have been analyzed. The first one refers to triads composed of a k = 0 first baroclinic geostrophic mode, which is resonant with the stationary component of the diurnal heat source, and two dispersive modes, namely, a mixed Rossby-gravity wave and a barotropic Rossby mode. The other class corresponds to triads composed of two first baroclinic inertio-gravity waves in which the highest-frequency wave resonates with a transient harmonic of the forcing. The integration of the asymptotic reduced equations for these selected resonant triads shows that the stationary component of the diurnal heat source acts as an ""accelerator"" for the energy exchanges between the two dispersive waves through the excitation of the catalyst geostrophic mode. On the other hand, since in the second class of triads the mode that resonates with the forcing is the most energetically active member because of the energy constraints imposed by the triad dynamics, the results show that the convective forcing in this case is responsible for a longer time scale modulation in the resonant interactions, generating a period doubling in the energy exchanges. The results suggest that the diurnal variation of tropical convection might play an important role in generating low-frequency fluctuations in the atmospheric circulation through resonant nonlinear interactions.

Resonant Wave Interactions in the Equatorial Waveguide

Weakly nonlinear interactions among equatorial waves have been explored in this paper using the adiabatic version of the equatorial beta-plane primitive equations in isobaric coordinates. Assuming rigid lid vertical boundary conditions, the conditions imposed at the surface and at the top of the troposphere were expanded in a Taylor series around two isobaric surfaces in an approach similar to that used in the theory of surface-gravity waves in deep water and capillary-gravity waves. By adopting the asymptotic method of multiple time scales, the equatorial Rossby, mixed Rossby-gravity, inertio-gravity, and Kelvin waves, as well as their vertical structures, were obtained as leading-order solutions. These waves were shown to interact resonantly in a triad configuration at the O(epsilon) approximation. The resonant triads whose wave components satisfy a resonance condition for their vertical structures were found to have the most significant interactions, although this condition is not excluding, unlike the resonant conditions for the zonal wavenumbers and meridional modes. Thus, the analysis has focused on such resonant triads. In general, it was found that for these resonant triads satisfying the resonance condition in the vertical direction, the wave with the highest absolute frequency always acts as an energy source (or sink) for the remaining triad components, as usually occurs in several other physical problems in fluid dynamics. In addition, the zonally symmetric geostrophic modes act as catalyst modes for the energy exchanges between two dispersive waves in a resonant triad. The integration of the reduced asymptotic equations for a single resonant triad shows that, for the initial mode amplitudes characterizing realistic magnitudes of atmospheric flow perturbations, the modes in general exchange energy on low-frequency (intraseasonal and/or even longer) time scales, with the interaction period being dependent upon the initial mode amplitudes. Potential future applications of the present theory to the real atmosphere with the inclusion of diabatic forcing, dissipation, and a more realistic background state are also discussed.

Coda wave attenuation in the Parecis Basin, Amazon Craton, Brazil: sensitivity to basement depth

Small local earthquakes from two aftershock sequences in Porto dos GaA(0)chos, Amazon craton-Brazil, were used to estimate the coda wave attenuation in the frequency band of 1 to 24 Hz. The time-domain coda-decay method of a single backscattering model is employed to estimate frequency dependence of the quality factor (Q (c)) of coda waves modeled usingwhere Q (0) is the coda quality factor at frequency of 1 Hz and eta is the frequency parameter. We also used the independent frequency model approach (Morozov, Geophys J Int, 175:239-252, 2008), based in the temporal attenuation coefficient, chi(f) instead of Q(f), given by the equation for the calculation of the geometrical attenuation (gamma) and effective attenuation Q (c) values have been computed at central frequencies (and band) of 1.5 (1-2), 3.0 (2-4), 6.0 (4-8), 9.0 (6-12), 12 (8-16), and 18 (12-24) Hz for five different datasets selected according to the geotectonic environment as well as the ability to sample shallow or deeper structures, particularly the sediments of the Parecis basin and the crystalline basement of the Amazon craton. For the Parecis basin for the surrounding shield and for the whole region of Porto dos GaA(0)chos Using the independent frequency model, we found: for the cratonic zone, gamma = 0.014 s (-aEuro parts per thousand 1), nu a parts per thousand 1.12; for the basin zone with sediments of similar to 500 m, gamma = 0.031 s (-aEuro parts per thousand 1), nu a parts per thousand 1.27; and for the Parecis basin with sediments of similar to 1,000 m, gamma = 0.047 s (-aEuro parts per thousand 1), nu a parts per thousand 1.42. Analysis of the attenuation factor (Q (c)) for different values of the geometrical spreading parameter (nu) indicated that an increase of nu generally causes an increase in Q (c), both in the basin as well as in the craton. But the differences in the attenuation between different geological environments are maintained for different models of geometrical spreading. It was shown that the energy of coda waves is attenuated more strongly in the sediments, (in the deepest part of the basin), than in the basement, (in the craton). Thus, the coda wave analysis can contribute to studies of geological structures in the upper crust, as the average coda quality factor is dependent on the thickness of sedimentary layer.

Dynamics of the viscous Cahn-Hilliard equation

We study generalized viscous Cahn-Hilliard problems with nonlinearities satisfying critical growth conditions in W-0(1,p)(Omega), where Omega is a bounded smooth domain in R-n, n >= 3. In the critical growth case, we prove that the problems are locally well posed and obtain a bootstrapping procedure showing that the solutions are classical. For p = 2 and almost critical dissipative nonlinearities we prove global well posedness, existence of global attractors in H-0(1)(Omega) and, uniformly with respect to the viscosity parameter, L-infinity(Omega) bounds for the attractors. Finally, we obtain a result on continuity of regular attractors which shows that, if n = 3, 4, the attractor of the Cahn-Hilliard problem coincides (in a sense to be specified) with the attractor for the corresponding semilinear heat equation. (C) 2008 Elsevier Inc. All rights reserved.

Low-dimensional chaos and wave turbulence in plasmas

We investigated drift-wave turbulence in the plasma edge of a small tokamak by considering solutions of the Hasegawa-Mima equation involving three interacting modes in Fourier space. The resulting low-dimensional dynamics presented periodic as well as chaotic evolution of the Fourier-mode amplitudes, and we performed the control of chaotic behaviour through the application of a fourth resonant wave of small amplitude.