967 resultados para Differential equations, Nonlinear -- Numerical solutions -- Computer programs
Resumo:
In this paper, we study the travelling wave reductions for certain (2 + 1)- and (3 + 1)-dimensional physically important nonlinear evolutionary equations by using the recently proposed Homogenous Balance Method (HBM). Through this analysis we explore certain new solutions for the equations we have studied. (C) 2001 Published by Elsevier B.V.
Resumo:
In this paper, we investigate the invariance and integrability properties of an integrable two-component reaction-diffusion equation. We perform Painleve analysis for both the reaction-diffusion equation modelled by a coupled nonlinear partial differential equations and its general similarity reduced ordinary differential equation and confirm its integrability. Further, we perform Lie symmetry analysis for this model. Interestingly our investigations reveals a rich variety of particular solutions, which have not been reported in the literature, for this model. (C) 2000 Elsevier B.V. Ltd. All rights reserved.
Resumo:
We examine two-component Gross-Pitaevskii equations with nonlinear and linear couplings, assuming self-attraction in one species and self-repulsion in the other, while the nonlinear inter-species coupling is also repulsive. For initial states with the condensate placed in the self-attractive component, a sufficiently strong linear coupling switches the collapse into decay (in the free space). Setting the linear-coupling coefficient to be time-periodic (alternating between positive and negative values, with zero mean value) can make localized states quasi-stable for the parameter ranges considered herein, but they slowly decay. The 2D states can then be completely stabilized by a weak trapping potential. In the case of the high-frequency modulation of the coupling constant, averaged equations are derived, which demonstrate good agreement with numerical solutions of the full equations. (C) 2007 Elsevier B.V. All rights reserved.
Resumo:
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.
Resumo:
The time evolution of the matter produced in high energy heavy-ion collisions seems to be well described by relativistic viscous hydrodynamics. In addition to the hydrodynamic degrees of freedom related to energy-momentum conservation, degrees of freedom associated with order parameters of broken continuous symmetries must be considered because they are all coupled to each other. of particular interest is the coupling of degrees of freedom associated with the chiral symmetry of QCD. Quantum and thermal fluctuations of the chiral fields act as noise sources in the classical equations of motion, turning them into stochastic differential equations in the form of Ginzburg-Landau-Langevin (GLL) equations. Analytic solutions of GLL equations are attainable only in very special circumstances and extensive numerical simulations are necessary, usually by discretizing the equations on a spatial lattice. However, a not much appreciated issue in the numerical simulations of GLL equations is that ultraviolet divergences in the form of lattice-spacing dependence plague the solutions. The divergences are related to the well-known Rayleigh-Jeans catastrophe in classical field theory. In the present communication we present a systematic lattice renormalization method to control the catastrophe. We discuss the implementation of the method for a GLL equation derived in the context of a model for the QCD chiral phase transition and consider the nonequilibrium evolution of the chiral condensate during the hydrodynamic flow of the quark-gluon plasma.
Resumo:
The investigation of the behavior of a nonlinear system consists in the analysis of different stages of its motion, where the complexity varies with the proximity of a resonance region. Near this region the stability domain of the system undergoes sudden changes due basically to competition and interaction between periodic and saddle solutions inside the phase portrait, leading to the occurrence of the most different phenomena. Depending of the domain of the chosen control parameter, these events can reveal interesting geometric features of the system so that the phase portrait is not capable to express all them, since the projection of these solutions on the two-dimensional surface can hide some aspects of these events. In this work we will investigate the numerical solutions of a particular pendulum system close to a secondary resonance region, where we vary the control parameter in a restrict domain in order to draw a preliminary identification about what happens with this system. This domain includes the appearance of non-hyperbolic solutions where the basin of attraction in the center of the phase portrait diminishes considerably, almost disappearing, and afterwards its size increases with the direction of motion inverted. This phenomenon delimits a boundary between low and high frequency of the external excitation.
Resumo:
In this paper, a nonideal mechanical system with the LuGre friction damping model is considered. The mechanical model of the system is an oscillator not necessarily linear connected with an unbalanced motor of excitation with limited power supply. The control of motion and the attenuation of the Sommerfeld effect of the considered nonideal system are analyzed in this paper The mathematical model of the system is represented by coupled non-linear differential equations. The identification of some interesting nonlinear phenomenon in the transient and steady state motion of the system during the passage through resonance (using applied voltages at dc motor as control parameter) is investigated in detail using numerical simulation. [DOI: 10.1115/1.3124783]
Resumo:
Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)
Resumo:
A simple mathematical model is developed to explain the appearance of oscillations in the dispersal of larvae from the food source in experimental populations of certain species of blowflies. The life history of the immature stage in these flies, and in a number of other insects, is a system with two populations, one of larvae dispersing on the soil and the other of larvae that burrow in the soil to pupate. The observed oscillations in the horizontal distribution of buried pupae at the end of the dispersal process are hypothesized to be a consequence of larval crowding at a given point in the pupation substrate. It is assumed that dispersing larvae are capable of perceiving variations in density of larvae buried at a given point in the substrate of pupation, and that pupal density may influence pupation of dispersing larvae. The assumed interaction between dispersing larvae and the larvae that are burrowing to pupate is modeled using the concept of non-local effects. Numerical solutions of integro-partial differential equations developed to model density-dependent immature dispersal demonstrate that variation in the parameter that governs the non-local interaction between dispersing and buried larvae induces oscillations in the final horizontal distribution of pupae. (C) 1997 Academic Press Limited.
Resumo:
Immobilized cell utilization in tower-type bioreactor is one of the main alternatives being studied to improve the industrial bioprocess. Other alternatives for the production of beta -lactam antibiotics, such as a cephalosporin C fed-batch p recess in an aerated stirred-tank bioreactor with free cells of Cepha-losporium acremonium or a tower-type bioreactor with immobilized cells of this fungus, have proven to be more efficient than the batch profess. In the fed-batch process, it is possible to minimize the catabolite repression exerted by the rapidly utilization of carbon sources (such as glucose) in the synthesis of antibiotics by utilizing a suitable flow rate of supplementary medium. In this study, several runs for cephalosporin C production, each lasting 200 h, were conducted in a fed-batch tower-type bioreactor using different hydrolyzed sucrose concentrations, For this study's model, modifications were introduced to take into account the influence of supplementary medium flow rate. The balance equations considered the effect of oxygen limitation inside the bioparticles. In the Monod-type rate equations, eel concentrations, substrate concentrations, and dissolved oxygen were included as reactants affecting the bioreaction rate. The set of differential equations was solved by the numerical method, and the values of the parameters were estimated by the classic nonlinear regression method following Marquardt's procedure with a 95% confidence interval. The simulation results showed that the proposed model fit well with the experimental data,and based on the experimental data and the mathematical model an optimal mass flow rate to maximize the bioprocess productivity could be proposed.
Resumo:
The spatial distribution of water and sugars in half-fresh apples dehydrated in sucrose solutions (30% and 50% w/w, 27 degrees C) for 2, 4 and 8 h, was determined. Each half was sliced as from the exposed surface. The density, water and sugar contents were determined for each piece. A mathematical model was fitted to the experimental data of the water and sucrose contents considering the overall flux and tissue shrinkage. A numerical method of finite differences permitted the calculation of the effective diffusion coefficients as a function of concentration, using material coordinates and integrating the two differential equations (for water and sucrose) simultaneously. The coefficients obtained were one or even two orders of magnitude lower than those for pure solutions and presented unusual concentration dependence. The behaviour of the apple tissue was also studied using light microscopy techniques to obtain images of the osmotically treated pieces (20%, 30% and 50% w/w sucrose solutions for 2, 4 and 8 h). (c) 2006 Elsevier Ltd. All rights reserved.
Resumo:
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.
Resumo:
A formulation used to determine the time-optimal geomagnetic attitude maneuvers subject to dynamic and geometric constraints is proposed in this paper. This was obtained by a direct search procedure based on a control function parametrization method, using linear programming to obtain numerical suboptimal solutions by linear perturbation. Due to its characteristics it can be used in small computers and to generate computer programs of general application. The dynamic modeling, the magnetic torque model and the suboptimal control procedure are presented. Simulation runs have verified the feasibility of the formulation thus derived and have shown a notable improvement in performance.
Resumo:
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.
Resumo:
A fourth-order numerical method for solving the Navier-Stokes equations in streamfunction/vorticity formulation on a two-dimensional non-uniform orthogonal grid has been tested on the fluid flow in a constricted symmetric channel. The family of grids is generated algebraically using a conformal transformation followed by a non-uniform stretching of the mesh cells in which the shape of the channel boundary can vary from a smooth constriction to one which one possesses a very sharp but smooth corner. The generality of the grids allows the use of long channels upstream and downstream as well as having a refined grid near the sharp corner. Derivatives in the governing equations are replaced by fourth-order central differences and the vorticity is eliminated, either before or after the discretization, to form a wide difference molecule for the streamfunction. Extra boundary conditions, necessary for wide-molecule methods, are supplied by a procedure proposed by Henshaw et al. The ensuing set of non-linear equations is solved using Newton iteration. Results have been obtained for Reynolds numbers up to 250 for three constrictions, the first being smooth, the second having a moderately sharp corner and the third with a very sharp corner. Estimates of the error incurred show that the results are very accurate and substantially better than those of the corresponding second-order method. The observed order of the method has been shown to be close to four, demonstrating that the method is genuinely fourth-order. © 1977 John Wiley & Sons, Ltd.