Mathematical modeling and control of population systems: Applications in biological pest control

The aim of this paper is to apply methods from optimal control theory, and from the theory of dynamic systems to the mathematical modeling of biological pest control. The linear feedback control problem for nonlinear systems has been formulated in order to obtain the optimal pest control strategy only through the introduction of natural enemies. Asymptotic stability of the closed-loop nonlinear Kolmogorov system is guaranteed by means of a Lyapunov function which can clearly be seen to be the solution of the Hamilton-Jacobi-Bellman equation, thus guaranteeing both stability and optimality. Numerical simulations for three possible scenarios of biological pest control based on the Lotka-Volterra models are provided to show the effectiveness of this method. (c) 2007 Elsevier B.V. All rights reserved.

Remarks on Parametric Surface Waves in a Nonlinear and Non-Ideally Excited Tank

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

Surface effects on vortex dynamics in thin type-II superconducting stripes with random pinning

We analyze the dynamics of a driven vortex lattice moving in a thin Superconducting stripe. The two dimensional stripe is assumed to be finite in the longitudinal direction, where we take into account the Surface effects, and infinite in the transversal direction. The numerical simulations are performed using the Langevin dynamics, including the vortex-vortex interaction, interaction of vortices with the surface current, vortex images, transport current and randomly distributed pinning centers. We show results for the differential resistivity and the vortex trajectories as a function of the external force. (C) 2004 Elsevier B.V. All rights reserved.

Clustering of rare earth in glasses, aluminum effect: experiments and modeling

Luminescent spectra of Eu3+-doped sol-gel glasses have been analyzed during the densification process and compared according to the presence or not of aluminum as a codoping ion. A transition temperature from hydrated to dehydroxyled environments has been found different for doped and codoped samples. However, only slight modifications have been displayed from luminescence measurements beyond this transition. To support the experimental analysis, molecular dynamics simulations have been performed to model the doped and codoped glass structures. Despite no evidence of rare earth clustering reduction due to aluminum has been found, the modeled structures have shown that the luminescent ions are mainly located in aluminum-rich domains. The synthesis of both experimental and numerical analyses has lead us to interpret the aluminum effect as responsible for differences in structure of the luminescent sites rather than for an effective dispersion of the rare earth ions. (C) 2004 Elsevier B.V. All rights reserved.

Numerical and experimental analysis of natural convection in a cavity heated from below

In this work, an analysis of the natural convection flow caused by heat sources dissipating energy at a constant rate simulating electronic components mounted at the bottom surface of a cavity symmetrically cooled from the sides and insulated at the top is performed. This problem was studied numerically and experimentally for several aspect ratios (height/width), for different levels of dissipation in the sources, and for different side wall temperatures. Temperature and velocity fields were determined as well as the temperature variation along the surface where the sources are mounted and the average Nusselt number in the source surfaces. Numerical and experimental results were found to agree.

Modeling and simulation of cephalosporin C production in a fed-batch tower-type bioreactor

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.

Surface effects on moving vortices in superconducting stripes

Size and surface dynamical effects are investigated in thin superconducting stripes with variable width. We perform numerical simulations of the vortex dynamics, with the inclusion of the surface confining potential and a random distribution of pinning centers. To fully characterize the vortex flow, we calculate the differential resistance, the transverse diffusion coefficient, the structure factor and the intensity of the Bragg peaks, as functions of the transport force. We found that surface effects induce a premature ordering of the flux line lattice, and the system displays plastic and smectic behavior only in a very narrow range of forces. (c) 2007 Elsevier B.V. All rights reserved.

Modeling hourly diffuse solar-radiation in the city of São Paulo using a neural-network technique

dIn this work, a perceptron neural-network technique is applied to estimate hourly values of the diffuse solar-radiation at the surface in São Paulo City, Brazil, using as input the global solar-radiation and other meteorological parameters measured from 1998 to 2001. The neural-network verification was performed using the hourly measurements of diffuse solar-radiation obtained during the year 2002. The neural network was developed based on both feature determination and pattern selection techniques. It was found that the inclusion of the atmospheric long-wave radiation as input improves the neural-network performance. on the other hand traditional meteorological parameters, like air temperature and atmospheric pressure, are not as important as long-wave radiation which acts as a surrogate for cloud-cover information on the regional scale. An objective evaluation has shown that the diffuse solar-radiation is better reproduced by neural network synthetic series than by a correlation model. (C) 2004 Elsevier Ltd. All rights reserved.

Numerical study of circulation on the inner Amazon Shelf

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

Cluster identification and characterization in the riser of a circulating fluidized bed from numerical simulation results

A methodology of identification and characterization of coherent structures mostly known as clusters is applied to hydrodynamic results of numerical simulation generated for the riser of a circulating fluidized bed. The numerical simulation is performed using the MICEFLOW code, which includes the two-fluids IIT's hydrodynamic model B. The methodology for cluster characterization that is used is based in the determination of four characteristics, related to average life time, average volumetric fraction of solid, existing time fraction and frequency of occurrence. The identification of clusters is performed by applying a criterion related to the time average value of the volumetric solid fraction. A qualitative rather than quantitative analysis is performed mainly owing to the unavailability of operational data used in the considered experiments. Concerning qualitative analysis, the simulation results are in good agreement with literature. Some quantitative comparisons between predictions and experiment were also presented to emphasize the capability of the modeling procedure regarding the analysis of macroscopic scale coherent structures. (c) 2007 Elsevier B.V. All rights reserved.

A case study of the 9 August 1988 South Atlantic storm: Numerical simulations of the wave activity

During 9-11 August 1988, a cyclone developed over Uruguay in the lee of the Andes Mountains and moved over the South Atlantic Ocean, where it redeveloped into an intense storm. This storm was responsible for unusual wave activity along the Brazilian shoreline from 22° to 32°S. The Brazilian news media reported the loss of at least one life, waves of 3 m and higher, and the disappearance of a drainage pipe, which weighed 8000 kg, off the shores of Rio de Janeiro. In this paper, the evolution of this intense storm and the associated ocean wave response is studied through European Centre for Medium-Range Weather Forecasts analyses, a hydrostatic limited-area meteorological model, and a second-generation prognostic wave model. The atmospheric model results indicated the presence of a long-lived and large fetch with surface wind velocities higher than 12 m s -1 directed toward the coast. Some areas with velocities of 20 m s -1 were embedded in the fetch. The wave model forced by this wind field was able to simulate waves with a significant height of 8 m far from the coast and about 4 m in regions very close to the Brazilian coast in agreement with the occurrence reported at Rio de Janeiro. The swell propagation toward the coast of Rio de Janeiro was obstructed by a northeastward 10-m wind during the first 24-h period of the model's integration. During the second 24-h period, the fetch was still large and strong, but the obstacle was removed by a counterclockwise rotation of wind direction favoring the swell and windsea propagation toward the Rio de Janeiro coast.

Fractional polynomial response surface models

Second-order polynomial models have been used extensively to approximate the relationship between a response variable and several continuous factors. However, sometimes polynomial models do not adequately describe the important features of the response surface. This article describes the use of fractional polynomial models. It is shown how the models can be fitted, an appropriate model selected, and inference conducted. Polynomial and fractional polynomial models are fitted to two published datasets, illustrating that sometimes the fractional polynomial can give as good a fit to the data and much more plausible behavior between the design points than the polynomial model. © 2005 American Statistical Association and the International Biometric Society.

On nonlinear dynamics in a free fluid surface of a tank excited by a non-ideal power source

We investigate the nonlinear oscillations in a free surface of a fluid in a cylinder tank excited by non-ideal power source, an electric motor with limited power supply. We study the possibility of parametric resonance in this system, showing that the excitation mechanism can generate chaotic response. Additionally, the dynamics of parametrically excited surface waves in the tank can reveal new characteristics of the system. The fluid-dynamic system is modeled in such way as to obtain a nonlinear differential equation system. Numerical experiments are carried out to find the regions of chaotic solutions. Simulation results are presented as phase-portrait diagrams characterizing the resonant vibrations of free fluid surface and the existence of several types of regular and chaotic attractors. We also describe the energy transfer in the interaction process between the hydrodynamic system and the electric motor. Copyright © 2011 by ASME.

Simulation of dynamic pinion course using Runge-Kutta's method and impact modeling

Once defined the relationship between the Starter Motor components and their functions, it is possible to develop a mathematical model capable to predict the Starter behavior during operation. One important aspect is the engagement system behavior. The development of a mathematical tool capable of predicting it is a valuable step in order to reduce the design time, cost and engineering efforts. A mathematical model, represented by differential equations, can be developed using physics laws, evaluating force balance and energy flow through the systems degrees of freedom. Another important physical aspect to be considered in this modeling is the impact conditions (particularly on the pinion and ring-gear contact). This work is a report of those equations application on available mathematical software and the resolution of those equations by Runge-Kutta's numerical integration method, in order to build an accessible engineering tool. Copyright © 2011 SAE International.

Mutual coupling modeling in transmission lines directly in the phase domain

This paper describes a computational model based on lumped elements for the mutual coupling between phases in three-phase transmission lines without the explicit use of modal transformation matrices. The self and mutual parameters and the coupling between phases are modeled using modal transformation techniques. The modal representation is developed from the intrinsic consideration of the modal transformation matrix and the resulting system of time-domain differential equations is described as state equations. Thus, a detailed profile of the currents and the voltages through the line can be easily calculated using numerical or analytical integration methods. However, the original contribution of the article is the proposal of a time-domain model without the successive phase/mode transformations and a practical implementation based on conventional electrical circuits, without the use of electromagnetic theory to model the coupling between phases. © 2011 IEEE.