An exact equation for the free surface of a fluid in a porous medium

We study the problem of the evolution of the free surface of a fluid in a saturated porous medium, bounded from below by a. at impermeable bottom, and described by the Laplace equation with moving-boundary conditions. By making use of a convenient conformal transformation, we show that the solution to this problem is equivalent to the solution of the Laplace equation on a fixed domain, with new variable coefficients, the boundary conditions. We use a kernel of the Laplace equation which allows us to write the Dirichlet-to-Neumann operator, and in this way we are able to find an exact differential-integral equation for the evolution of the free surface in one space dimension. Although not amenable to direct analytical solutions, this equation turns out to allow an easy numerical implementation. We give an explicit illustrative case at the end of the article.

The rigid-flexible nonlinear robotic manipulator: Modeling and control

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

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.

Host-guest complexation of a nitroheterocyclic compound with cyclodextrins: a spectrofluorimetric and molecular modeling study

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

Turbulence modeling in engineering applications

A brief review is given of turbulence models in use today for engineering applications. The main categories covered are simple eddy-viscosity models, the k-ε two-equation model and Reynolds-stress-equation models as well as their algebraic stress derivatives. Calculation examples are presented for a variety of 2D and 3D flows.

Object oriented modeling of piston engines

This work reports a conception phase of a piston engine global model. The model objective is forecast the motor performance (power, torque and specific consumption as a function of rotation and environmental conditions). Global model or Zero-dimensional is based on flux balance through each engine component. The resulting differential equations represents a compressive unsteady flow, in which, all dimensional variables are areas or volumes. A review is presented first. The ordinary differential equation system is presented and a Runge-Kutta method is proposed to solve it numerically. The model includes the momentum conservation equation to link the gas dynamics with the engine moving parts rigid body mechanics. As an oriented to objects model the documentation follows the UML standard. A discussion about the class diagrams is presented, relating the classes with physical model related. The OOP approach allows evolution from simple models to most complex ones without total code rewrite. Copyright © 2001 Society of Automotive Engineers, Inc.

Modeling approach based on experimental results for prediction of measurement errors in energy meters

This paper presents a general modeling approach to investigate and to predict measurement errors in active energy meters both induction and electronic types. The measurement error modeling is based on Generalized Additive Model (GAM), Ridge Regression method and experimental results of meter provided by a measurement system. The measurement system provides a database of 26 pairs of test waveforms captured in a real electrical distribution system, with different load characteristics (industrial, commercial, agricultural, and residential), covering different harmonic distortions, and balanced and unbalanced voltage conditions. In order to illustrate the proposed approach, the measurement error models are discussed and several results, which are derived from experimental tests, are presented in the form of three-dimensional graphs, and generalized as error equations. © 2009 IEEE.

Modeling the growth of LT and TL-oriented fatigue cracks in longitudinally and transversely pre-strained Al 2524-T3 alloy

The aluminum alloy 2524 (Al-Cu-Mg) was developed during the 90s mainly to be employed in aircraft fuselage panels, replacing the standard Al 2024. In the present analysis the fatigue crack growth (FCG) behavior of 2524-T3 was investigated, regarding the influence of three parameters: load ratio, pre strain and crack plane orientation of the material. The pre strain of aluminum alloys is usually performed in order to obtain a more homogeneous precipitates distribution, accompanied by an increase in the yield strength. In this work, it was evaluated the resistance of Al 2524-T3 sheet samples to the fatigue crack growth, having L-T and T-L crack orientations. FCG tests were performed under constant amplitude loading at three distinct positive load ratios. The three material conditions were tested: as received(AR), pre strained longitudinally (SL) and transversally (ST) in relation to rolling direction. In order to describe FCG behavior, two-parameter kinetic equations were compared: a Paris-type potential model and a new exponential equation introduced in a previous work conducted by our research group. It was observed that the exponential model, which takes into account the deviations from linearity presented by da/dN versus AK data, describes more adequately the FCG behavior of Al 224-T3 in relation to load ratio, pre strain effects and crack plane orientation. © 2011 Published by Elsevier Ltd.

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.

ECOLOGICAL MODELING OF TALITROIDES TOPITOTUM (CRUSTACEA: AMPHIPODA)

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

About the minimum time transference in a fractional differential equation

The use of fractional calculus when modeling phenomena allows new queries concerning the deepest parts of the physical laws involved in. Here we will be dealing with an apparent paradox in which the time of transference from zero in a system with fractional derivatives can be strictly shortened relatively to the minimal time transference done in an equivalent system in the frame of the entire derivatives.

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

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

Equation of Motion Governing the Dynamics of Vertically Collapsing Buildings

The present paper aims at contributing to a discussion, opened by several authors, on the proper equation of motion that governs the vertical collapse of buildings. The most striking and tragic example is that of the World Trade Center Twin Towers, in New York City, about 10 years ago. This is a very complex problem and, besides dynamics, the analysis involves several areas of knowledge in mechanics, such as structural engineering, materials sciences, and thermodynamics, among others. Therefore, the goal of this work is far from claiming to deal with the problem in its completeness, leaving aside discussions about the modeling of the resistive load to collapse, for example. However, the following analysis, restricted to the study of motion, shows that the problem in question holds great similarity to the classic falling-chain problem, very much addressed in a number of different versions as the pioneering one, by von Buquoy or the one by Cayley. Following previous works, a simple single-degree-of-freedom model was readdressed and conceptually discussed. The form of Lagrange's equation, which leads to a proper equation of motion for the collapsing building, is a general and extended dissipative form, which is proper for systems with mass varying explicitly with position. The additional dissipative generalized force term, which was present in the extended form of the Lagrange equation, was shown to be derivable from a Rayleigh-like energy function. DOI: 10.1061/(ASCE)EM.1943-7889.0000453. (C) 2012 American Society of Civil Engineers.

Modeling the exergy behavior of human body

Exergy analysis is applied to assess the energy conversion processes that take place in the human body, aiming at developing indicators of health and performance based on the concepts of exergy destroyed rate and exergy efficiency. The thermal behavior of the human body is simulated by a model composed of 15 cylinders with elliptical cross section representing: head, neck, trunk, arms, forearms, hands, thighs, legs, and feet. For each, a combination of tissues is considered. The energy equation is solved for each cylinder, being possible to obtain transitory response from the body due to a variation in environmental conditions. With this model, it is possible to obtain heat and mass flow rates to the environment due to radiation, convection, evaporation and respiration. The exergy balances provide the exergy variation due to heat and mass exchange over the body, and the exergy variation over time for each compartments tissue and blood, the sum of which leads to the total variation of the body. Results indicate that exergy destroyed and exergy efficiency decrease over lifespan and the human body is more efficient and destroys less exergy in lower relative humidities and higher temperatures. (C) 2012 Elsevier Ltd. All rights reserved.

Generalized Metropolis dynamics with a generalized master equation: An approach for time-independent and time-dependent Monte Carlo simulations of generalized spin systems

The extension of Boltzmann-Gibbs thermostatistics, proposed by Tsallis, introduces an additional parameter q to the inverse temperature beta. Here, we show that a previously introduced generalized Metropolis dynamics to evolve spin models is not local and does not obey the detailed energy balance. In this dynamics, locality is only retrieved for q = 1, which corresponds to the standard Metropolis algorithm. Nonlocality implies very time-consuming computer calculations, since the energy of the whole system must be reevaluated when a single spin is flipped. To circumvent this costly calculation, we propose a generalized master equation, which gives rise to a local generalized Metropolis dynamics that obeys the detailed energy balance. To compare the different critical values obtained with other generalized dynamics, we perform Monte Carlo simulations in equilibrium for the Ising model. By using short-time nonequilibrium numerical simulations, we also calculate for this model the critical temperature and the static and dynamical critical exponents as functions of q. Even for q not equal 1, we show that suitable time-evolving power laws can be found for each initial condition. Our numerical experiments corroborate the literature results when we use nonlocal dynamics, showing that short-time parameter determination works also in this case. However, the dynamics governed by the new master equation leads to different results for critical temperatures and also the critical exponents affecting universality classes. We further propose a simple algorithm to optimize modeling the time evolution with a power law, considering in a log-log plot two successive refinements.