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.

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.

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.

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.

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.

Basic ingredients for mathematical modeling of tumor growth in vitro: Cooperative effects and search for space

Based on the literature data from HT-29 cell monolayers, we develop a model for its growth, analogous to an epidemic model, mixing local and global interactions. First, we propose and solve a deterministic equation for the progress of these colonies. Thus, we add a stochastic (local) interaction and simulate the evolution of an Eden-like aggregate by using dynamical Monte Carlo methods. The growth curves of both deterministic and stochastic models are in excellent agreement with the experimental observations. The waiting times distributions, generated via our stochastic model, allowed us to analyze the role of mesoscopic events. We obtain log-normal distributions in the initial stages of the growth and Gaussians at long times. We interpret these outcomes in the light of cellular division events: in the early stages, the phenomena are dependent each other in a multiplicative geometric-based process, and they are independent at long times. We conclude that the main ingredients for a good minimalist model of tumor growth, at mesoscopic level, are intrinsic cooperative mechanisms and competitive search for space. © 2013 Elsevier Ltd.

Modeling of piezoelectric energy harvesting considering the dependence of the rectifier circuit

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

Modeling of multiple cracks in reinforced concrete members using solid finite elements with high aspect ratio

This paper presents a new technique to model interfaces by means of degenerated solid finite elements, i.e., elements with a very high aspect ratio, with the smallest dimension corresponding to the thickness of the interfaces. It is shown that, as the aspect ratio increases, the element strains also increase, approaching the kinematics of the strong discontinuity. A tensile damage constitutive relation between strains and stresses is proposed to describe the nonlinear behavior of the interfaces associated with crack opening. To represent crack propagation, couples of triangular interface elements are introduced in between all regular (bulk) elements of the original mesh. With this technique the analyses can be performed integrally in the context of the continuum mechanics and complex crack patterns involving multiple cracks can be simulated without the need of tracking algorithms. Numerical tests are performed to show the applicability of the proposed technique, studding also aspects related to mesh objectivity.

Modeling of As2S3 raman fiber lasers operating in the mid-infrared

A numerical study of As2S3 Raman fiber lasers is carried out to show their potential for the entire coverage of the 3–4-m spectral band. Experimental results are first obtained from such a laser operated under controlled conditions in order to set the fiber parameters (i.e., gain and attenuation coefficients) to be used in the numerical model. An exhaustive numerical analysis is then performed to establish the conditions for optimum lasing performances over the entire 3–4-m spectral band.

Modeling field load tests in lateritic unsaturated soil

Many models for unsaturated soil have been developed in the last years, accompanying the development of experimental techniques to deal with such soils. The benchmark of the models for unsaturated soil can be assigned to the Barcelona Basic Model (BBM) now incorporated in some codes such as the CODE_BRIGHT. Most of those models were validated considering limited laboratory test results and not much validation is available considering real field problems. This paper presents modeling results of field plate load tests performed under known suction on a lateritic unsaturated soil. The required input data were taken from laboratory tests performed under suction control. The modeling nicely reproduces field tests allowing appreciating the influence of soil suction on the stress-settlement curve. In addition, wetting induced or collapse settlements were calculated from field tests and were nicely duplicated by the numerical analysis performed.

Graded meshes in bio-thermal problems with transmission-line modeling method

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

Numerical analysis of a modified evacuated tubes solar collector

The need for renewable energy sources, facing the consequences of Climate Change, results in growing investment for solar collectors’ use. Research in this field has accompanied this expansion and evacuated tube solar collector stands as an important study focus. Thus, several works have been published for representing the stratification of the fluid inside the tubes and the reservoir, as well as analytical modeling for the heat flow problem. Based on recent publications, this paper proposes the study of solar water heating with evacuated tubes, their operation characteristics and operating parameters. To develop this work, a computational tool will be used - in this case, the application of computational fluid dynamics (CFD) software. In possession of the implemented model, a numerical simulation will be performed to evaluate the behavior of the fluid within this solar collector and possible improvements to be applied in the model.

Failure behavior modeling of slender reinforced concrete columns subjected to eccentric load

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