Non-linear growth models for bullfrog tadpoles

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

Iterated deferred correction for linear two-point boundary value problems

An analysis of iterated deferred correction based on various classes of implicit Runge-Kutta formulae is given. Out of different possibilities considered, it is shown that those based purely on Lobatto formulae have the best stability. The enhanced stability of Lobatto schemes is very important for the efficient integration of excessively stiff boundary value problems and this is demonstrated by means of some numerical results.

Effect of tillage on fractal indices describing soil surface microrelief of a Brazilian Alfisol

Soil surface roughness is known to influence water infiltration, runoff and erosion. Soil surface roughness changes with management and weather and its mathematical description still remains an important issue. The main objective of this study was to investigate the effect of tillage on the two fractal indices, fractal dimension, D, and crossover length, 1, currently used in characterizing soil surface microrelief. The statistical index random roughness, RR, was also assessed. Field experiments were done on an Alfisol located at Rio Grande do Sul State (Brazil). Two tillage treatments (conventional versus direct drilling) were tested. The soil surface microrelief was assessed by point elevation measurements in 16 plots for each treatment. The sampling scheme was a square grid with 20 x 20 mm between point spacing and the plot size was 280 x 280 mm, so that each data set consisted of 225 individual elevation points. All indices were calculated after trend removal, both by slope correction, i.e., oriented microrelief, and by slope plus tillage marks correction, i.e., random microrelief. The implemented algorithm for estimating D and 1 consisted in evaluating the roughness around the local root mean square deviation (RMS) of the point elevation values. Irrespective of tillage treatment and detrending procedure, fractal behavior extended only over a bounded range of scales, from 40 to 100 mm, due to the experimental setup. In these conditions, assessing fractal indices was not always straightforward. The statistical index RR and the fractal index I were significantly different between tillage treatments for oriented and random surface conditions. D values of random soil surfaces were not affected by tillage treatment, whereas D values of oriented microrelief were significantly lower in the direct drilled plots. Removal of tillage marks trend resulted in a significant increase in D values. Within each tillage treatment, 1 and D were significantly correlated. (c) 2006 Elsevier B.V. All rights reserved.

Model for electrodes' filaments of hot cathode fluorescent lamps, during preheating with constant rms current

This paper presents a new model for the representation of electrodes' filaments of hot-cathode fluorescent lamps, during preheating processes based on the injection of currents with constant root mean square (rms) values. The main improvement obtained with this model is the prediction of the R-h/R-c ratio during the preheating process, as a function of the preheating time and of the rms current injected in the electrodes. Using the proposed model, it is possible to obtain an estimate of the time interval and the current that should be provided by the electronic ballast, in order to ensure a suitable preheating process. is estimate of time and current can be used as input data in the design of electronic ballasts with programmed lamp start, permitting the prediction of the R-h/R-c ratio during the initial steps of the design (theoretical analysis and digital simulation). Therefore, the use of the proposed model permits to reduce the necessity of several empirical adjustments in the prototype, in order to set the operation of electronic ballasts during the preheating process. This fact reduces time and costs associated to the global design procedure of electronic ballasts.

Viscoelastic behavior of persimmons dried at constant air temperature

The viscoelastic behavior of dried persimmons at different air-drying temperatures and velocities was evaluated. Air temperatures and velocities were varied according to a second-order central composite design, with temperature ranging from 40degreesC to 70degreesC and air velocity from 0.8 to 2.0 m/s. After drying, persimmons were equilibrated at four different water activities: 0.432, 0.576, 0.625 and 0.751. The rheological behavior of dried and conditioned persimmons was studied under uniaxial compression-relaxation tests. Three different rheological models were fitted to the experimental relaxation curves: Maxwell, Generalized Maxwell and Peleg and Normand. Based on the root mean square of residuals, the Generalized Maxwell model showed the best fit and a regression analysis was applied to obtain response surfaces for the model parameters. The dependence of the rheological properties on water activity was also analysed. Results showed that only the linear effect of air temperature was significant at a 5% level on the equilibrium stress and relaxation times. In a general way, these parameters increased with increasing air temperature and decreasing water activity. (C) 2004 Swiss Society of Food Science and Technology. Published by Elsevier Ltd. All rights reserved.

Validity and reliability of the Portuguese version of Mandibular Function Impairment Questionnaire

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

INFLUENCE of HIP EXTERNAL ROTATION on HIP ADDUCTOR and RECTUS FEMORIS MYOELECTRIC ACTIVITY DURING A DYNAMIC PARALLEL SQUAT

Pereira, GR, Leporace, G, Chagas, DV, Furtado, LFL, Praxedes, J, and Batista, LA. Influence of hip external rotation on hip adductor and rectus femoris myoelectric activity during a dynamic parallel squat. J Strength Cond Res 24(10): 27492754, 2010-This study sought to compare the myoelectric activity of the hip adductors (HAs) and rectus femoris (RF) when the hip was in a neutral position or externally rotated by 30 degrees or 50 degrees (H0, H30, and H50, respectively) during a parallel squat. Ten healthy subjects performed 10 repetitions of squats in each of the 3 hip positions and the myoelectric activities of the HAs and RF were recorded. The signal was then divided into categories representing concentric (C) and eccentric (E) contractions in the following ranges of motion: 0-30 degrees (C1 and E1), 30-60 degrees (C2 and E2), and 60-90 degrees (C3 and E3) of knee flexion. From those signals, an root mean square (RMS) value for each range of motion in each hip position was obtained. All values were normalized to those obtained during maximum voluntary isometric contraction. We found that HAs showed a significant increase in myoelectric activity during C3 and E3 in the H30 and H50 positions, as compared with H0. Meanwhile, RF activity did not significantly differ between hip positions. Both muscles showed higher activation during 60-90 degrees (C3 and E3) of knee flexion, as compared with 0-30 degrees (C1 and E1) and 30-60 degrees (C2 and E2). The results suggest that if the aim is to increase HA activity despite the low percentage of muscle activation, squats should be performed with 30 degrees of external rotation and at least 90 degrees of knee flexion.

Tm and Tm-Tb-doped germanate glasses for S-band amplifiers

The mechanism involved in the Tm(3+)((3)F(4)) -> Tb(3+)((7)F(0,1,2)) energy transfer as a function of the Tb concentration was investigated in Tm:Tb-doped germanate (GLKZ) glass. The experimental transfer rate was determined from the best fit of the (3)F(4) luminescence decay due to the Tm -> Tb energy transfer using the Burshtein model. The result showed that the 1700 nm emission from (3)F(4) can be completely quenched by 0.8 mol% of Tb(3+). As a consequence, the (7)F(3) state of Tb(3+) interacts with the (3)H(4) upper excited state of TM(3+) slighting decreasing its population. The effective amplification coefficient beta(cm(-1)) that depends on the population density difference Delta n = n((3)H(4))-n((3)F(4)) involved in the optical transition of Tm(3+) (S-band) was calculated by solving the rate equations of the system for continuous pumping with laser at 792 nm, using the Runge-Kutta numerical method including terms of fourth order. The population density inversion An as a function of Tb(3+) concentration was calculated by computational simulation for three pumping intensities, 0.2, 2.2 and 4.4 kWcm(-2). These calculations were performed using the experimental Tm -> Tb transfer rates and the optical constants of the Tm (0.1 mol%) system. It was demonstrated that 0.2 mol% of Tb(3+) propitiates best population density inversion of Tin(3+) maximizing the amplification coefficient of Tm-doped (0.1 mol%) GLKZ glass when operating as laser intensity amplification at 1.47 mu m. (C) 2007 Elsevier B.V. All rights reserved.

Integration Methods Used in Numerical Simulations of Electromagnetic Transients

This paper made an analysis of some numerical integration methods that can be used in electromagnetic transient simulations. Among the existing methods, we analyzed the trapezoidal integration method (or Heun formula), Simpson's Rule and Runge-Kutta. These methods were used in simulations of electromagnetic transients in power systems, resulting from switching operations and maneuvers that occur in transmission lines. Analyzed the characteristics such as accuracy, computation time and robustness of the methods of integration.

An evaluation of ocean color model estimates of marine primary productivity in coastal and pelagic regions across the globe

Nearly half of the earth's photosynthetically fixed carbon derives from the oceans. To determine global and region specific rates, we rely on models that estimate marine net primary productivity (NPP) thus it is essential that these models are evaluated to determine their accuracy. Here we assessed the skill of 21 ocean color models by comparing their estimates of depth-integrated NPP to 1156 in situ C-14 measurements encompassing ten marine regions including the Sargasso Sea, pelagic North Atlantic, coastal Northeast Atlantic, Black Sea, Mediterranean Sea, Arabian Sea, subtropical North Pacific, Ross Sea, West Antarctic Peninsula, and the Antarctic Polar Frontal Zone. Average model skill, as determined by root-mean square difference calculations, was lowest in the Black and Mediterranean Seas, highest in the pelagic North Atlantic and the Antarctic Polar Frontal Zone, and intermediate in the other six regions. The maximum fraction of model skill that may be attributable to uncertainties in both the input variables and in situ NPP measurements was nearly 72%. on average, the simplest depth/wavelength integrated models performed no worse than the more complex depth/wavelength resolved models. Ocean color models were not highly challenged in extreme conditions of surface chlorophyll-a and sea surface temperature, nor in high-nitrate low-chlorophyll waters. Water column depth was the primary influence on ocean color model performance such that average skill was significantly higher at depths greater than 250 m, suggesting that ocean color models are more challenged in Case-2 waters (coastal) than in Case-1 (pelagic) waters. Given that in situ chlorophyll-a data was used as input data, algorithm improvement is required to eliminate the poor performance of ocean color NPP models in Case-2 waters that are close to coastlines. Finally, ocean color chlorophyll-a algorithms are challenged by optically complex Case-2 waters, thus using satellite-derived chlorophyll-a to estimate NPP in coastal areas would likely further reduce the skill of ocean color models.

A diquark model for baryon spectroscopy

A quark-diquark approximation is used to investigate the mass spectroscopy of the spin-1/2 baryons belonging to the SU(3)-flavor group in a nonrelativistic potential approach. The baryon spectra obtained are confronted with relativistic results and experimental data. Root-mean-square radii are also calculated. © 1993 Springer-Verlag.

Another prefiltering algorithm for adaptive IIR filtering

An algorithm for adaptive IIR filtering that uses prefiltering structure in direct form is presented. This structure has an estimation error that is a linear function of the coefficients. This property greatly simplifies the derivation of gradient-based algorithms. Computer simulations show that the proposed structure improves convergence speed.

Phonons in binary BCC lithium-sodium alloys

The dispersion relations along the principal symmetry directions in BCC lithium-sodium alloys are calculated using second-order perturbation theory. The local modified Hoshino-Youngmodel potential was used for the lithium and the local Harrison model potential for sodium. The phonon density of states, the root mean square displacements and (Θ-T) curves are also calculated. In the absence of experimental data, just the theoretical predictions are presented here.

Improved numerical approach for the time-independent Gross-Pitaevskii nonlinear Schrödinger equation

In the present work, we improve a numerical method, developed to solve the Gross-Pitaevkii nonlinear Schrödinger equation. A particular scaling is used in the equation, which permits us to evaluate the wave-function normalization after the numerical solution. We have a two-point boundary value problem, where the second point is taken at infinity. The differential equation is solved using the shooting method and Runge-Kutta integration method, requiring that the asymptotic constants, for the function and its derivative, be equal for large distances. In order to obtain fast convergence, the secant method is used. © 1999 The American Physical Society.

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.