Periodic time series modeling of environmental proxy records with guaranteed positive growth rate estimation

Identifying a periodic time-series model from environmental records, without imposing the positivity of the growth rate, does not necessarily respect the time order of the data observations. Consequently, subsequent observations, sampled in the environmental archive, can be inversed on the time axis, resulting in a non-physical signal model. In this paper an optimization technique with linear constraints on the signal model parameters is proposed that prevents time inversions. The activation conditions for this constrained optimization are based upon the physical constraint of the growth rate, namely, that it cannot take values smaller than zero. The actual constraints are defined for polynomials and first-order splines as basis functions for the nonlinear contribution in the distance-time relationship. The method is compared with an existing method that eliminates the time inversions, and its noise sensitivity is tested by means of Monte Carlo simulations. Finally, the usefulness of the method is demonstrated on the measurements of the vessel density, in a mangrove tree, Rhizophora mucronata, and the measurement of Mg/Ca ratios, in a bivalve, Mytilus trossulus.

(Non-)ergodicity of a degenerate diffusion modeling the fiber lay down process

We analyze the large time behavior of a stochastic model for the lay down of fibers on a moving conveyor belt in the production process of nonwovens. It is shown that under weak conditions this degenerate diffusion process has a unique invariant distribution and is even geometrically ergodic. This generalizes results from previous works [M. Grothaus and A. Klar, SIAM J. Math. Anal., 40 (2008), pp. 968–983; J. Dolbeault et al., arXiv:1201.2156] concerning the case of a stationary conveyor belt, in which the situation of a moving conveyor belt has been left open.

Sparse shape representation using the Laplace-Beltrami eigenfunctions and its application to modeling subcortical structures

We present a new sparse shape modeling framework on the Laplace-Beltrami (LB) eigenfunctions. Traditionally, the LB-eigenfunctions are used as a basis for intrinsically representing surface shapes by forming a Fourier series expansion. To reduce high frequency noise, only the first few terms are used in the expansion and higher frequency terms are simply thrown away. However, some lower frequency terms may not necessarily contribute significantly in reconstructing the surfaces. Motivated by this idea, we propose to filter out only the significant eigenfunctions by imposing l1-penalty. The new sparse framework can further avoid additional surface-based smoothing often used in the field. The proposed approach is applied in investigating the influence of age (38-79 years) and gender on amygdala and hippocampus shapes in the normal population. In addition, we show how the emotional response is related to the anatomy of the subcortical structures.

Understanding the transition to seizure by modeling the epileptiform activity of general anesthetic agents

A central difficulty in modeling epileptogenesis using biologically plausible computational and mathematical models is not the production of activity characteristic of a seizure, but rather producing it in response to specific and quantifiable physiologic change or pathologic abnormality. This is particularly problematic when it is considered that the pathophysiological genesis of most epilepsies is largely unknown. However, several volatile general anesthetic agents, whose principle targets of action are quantifiably well characterized, are also known to be proconvulsant. The authors describe recent approaches to theoretically describing the electroencephalographic effects of volatile general anesthetic agents that may be able to provide important insights into the physiologic mechanisms that underpin seizure initiation.

Modeling the summertime evolution of sea-ice melt ponds

1] We present a mathematical model describing the summer melting of sea ice. We simulate the evolution of melt ponds and determine area coverage and total surface ablation. The model predictions are tested for sensitivity to the melt rate of unponded ice, enhanced melt rate beneath the melt ponds, vertical seepage, and horizontal permeability. The model is initialized with surface topographies derived from laser altimetry corresponding to first-year sea ice and multiyear sea ice. We predict that there are large differences in the depth of melt ponds and the area of coverage between the two types of ice. We also find that the vertical seepage rate and the melt rate of unponded ice are important in determining the total surface ablation and area covered by melt ponds.

Modeling migratory grazing of zooplankton on toxic and non-toxic phytoplankton

Migratory grazing of zooplankton between non-toxic phytoplankton (NTP) and toxic phytoplankton (TPP) is a realistic phenomena unexplored so far. The present article is a first step in this direction. A mathematical model of NTP–TPP-zooplankton with constant and variable zooplankton migration is proposed and analyzed. The asymptotic dynamics of the model system around the biologically feasible equilibria is explored through local stability analysis. The dynamics of the proposed system is explored and displayed for different combination of migratory parameters and toxin inhibition parameters. Our analysis suggests that the migratory grazing of zooplankton has a significant role in determining the dynamic stability and oscillation of phytoplankton zooplankton systems.

An evaluation of three upwinding approximations for numerical modeling the flow in tubular membrane of Newtonian and non-Newtonian fluids

In this paper, the laminar fluid flow of Newtonian and non-Newtonian of aqueous solutions in a tubular membrane is numerically studied. The mathematical formulation, with associated initial and boundary conditions for cylindrical coordinates, comprises the mass conservation, momentum conservation and mass transfer equations. These equations are discretized by using the finite-difference technique on a staggered grid system. Comparisons of the three upwinding schemes for discretization of the non-linear (convective) terms are presented. The effects of several physical parameters on the concentration profile are investigated. The numerical results compare favorably with experimental data and the analytical solutions. (C) 2011 Elsevier Inc. All rights reserved.

Quantitative structure-activity relationships for a series of inhibitors of cruzain from Trypanosoma cruzi: Molecular modeling, CoMFA and CoMSIA studies

Human parasitic diseases are the foremost threat to human health and welfare around the world. Trypanosomiasis is a very serious infectious disease against which the currently available drugs are limited and not effective. Therefore, there is an urgent need for new chemotherapeutic agents. One attractive drug target is the major cysteine protease from Trypanosoma cruzi, cruzain. In the present work, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) studies were conducted on a series of thiosemicarbazone and semicarbazone derivatives as inhibitors of cruzain. Molecular modeling studies were performed in order to identify the preferred binding mode of the inhibitors into the enzyme active site, and to generate structural alignments for the three-dimensional quantitative structure-activity relationship (3D QSAR) investigations. Statistically significant models were obtained (CoMFA. r(2) = 0.96 and q(2) = 0.78; CoMSIA, r(2) = 0.91 and q(2) = 0.73), indicating their predictive ability for untested compounds. The models were externally validated employing a test set, and the predicted values were in good agreement with the experimental results. The final QSAR models and the information gathered from the 3D CoMFA and CoMSIA contour maps provided important insights into the chemical and structural basis involved in the molecular recognition process of this family of cruzain inhibitors, and should be useful for the design of new structurally related analogs with improved potency. (C) 2009 Elsevier Inc. All rights reserved.

Open Circuit Interaction of Formic Acid with Oxidized Pt Surfaces: Experiments, Modeling, and Simulations

We report in this work the study of the interaction between formic acid and an oxidized platinum surface under open circuit conditions. The investigation was carried out with the aid of in situ infrared spectroscopy, and results analyzed in terms of a mathematical model and numerical simulations. It has been found that during the first seconds of the interaction a small amount of CO(2) is produced and absolutely no adsorbed CO was observed. A sudden drop in potential then follows, which is accompanied by a steep increase first of CO(2) production and then by adsorbed CO. The steep transient was rationalized in terms of an autocatalytic production of free platinum sites which enhances the overall rate of reaction. Modeling and simulation showed nearly quantitative agreement with the experimental observations and provided further insight into some experimentally inaccessible variables such as surface free sites. Finally, based on the understanding provided from the combined experimental and theoretical approach, we discuss the general aspects influencing the open circuit transient.

Modeling upper eyelid kinematics during spontaneous and reflex blinks

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

Ecological modeling and pest population management: a possible and necessary connection in a changing world

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

Population dynamics, life stage and ecological modeling in Chrysomya albiceps (Wiedemann) (Diptera: Calliphoridae)

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

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.

Modeling and identification of fertility maps using artificial neural networks

The application of agricultural fertilizers using variable rates along the field can be made through fertility maps previously elaborated or through real-time sensors. In most of the cases applies maps previously elaborated. These maps are identified from analyzes done in soil samples collected regularly (a sample for each field cell) or irregularly along the field. At the moment, mathematical interpolation methods such as nearest neighbor, local average, weighted inverse distance, contouring and kriging are used for predicting the variables involved with elaboration of fertility maps. However, some of these methods present deficiencies that can generate different fertility maps for a same data set. Moreover, such methods can generate inprecise maps to be used in precision farming. In this paper, artificial neural networks have been applied for elaboration and identification of precise fertility maps which can reduce the production costs and environmental impacts.

Geometric modeling from the medical images in the CAD system to support prosthesis design

This paper presents an individual designing prosthesis for surgical use and proposes a methodology for such design through mathematical extrapolation of data from digital images obtained via tomography of individual patient's bones. Individually tailored prosthesis designed to fit particular patient requirements as accurately as possible should result in more successful reconstruction, enable better planning before surgery and consequently fewer complications during surgery. Fast and accurate design and manufacture of personalized prosthesis for surgical use in bone replacement or reconstruction is potentially feasible through the application and integration of several different existing technologies, which are each at different stages of maturity. Initial case study experiments have been undertaken to validate the research concepts by making dimensional comparisons between a bone and a virtual model produced using the proposed methodology and a future research directions are discussed.