Mathematical models of generalized diffusion

We discuss in this paper equations describing processes involving non-linear and higher-order diffusion. We focus on a particular case (u(t) = 2 lambda (2)(uu(x))(x) + lambda (2)u(xxxx)), which is put into analogy with the KdV equation. A balance of nonlinearity and higher-order diffusion enables the existence of self-similar solutions, describing diffusive shocks. These shocks are continuous solutions with a discontinuous higher-order derivative at the shock front. We argue that they play a role analogous to the soliton solutions in the dispersive case. We also discuss several physical instances where such equations are relevant.

Higher order Painleve equations and their symmetries via reductions of a class of integrable models

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

Path integral quantization of generalized quantum electrodynamics

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

The canonical structure of Podolsky's generalized electrodynamics on the null-plane

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

A generalized alternating-direction implicit scheme for incompressible magnetohydrodynamic viscous flows at low magnetic Reynolds number

This paper presents numerical simulations of incompressible fluid flows in the presence of a magnetic field at low magnetic Reynolds number. The equations governing the flow are the Navier-Stokes equations of fluid motion coupled with Maxwell's equations of electromagnetics. The study of fluid flows under the influence of a magnetic field and with no free electric charges or electric fields is known as magnetohydrodynamics. The magnetohydrodynamics approximation is considered for the formulation of the non-dimensional problem and for the characterization of similarity parameters. A finite-difference technique is used to discretize the equations. In particular, an extension of the generalized Peaceman and Rachford alternating-direction implicit (ADI) scheme for simulating two-dimensional fluid flows is presented. The discretized conservation equations are solved in stream function-vorticity formulation. We compare the ADI and generalized ADI schemes, and show that the latter is more efficient in simulating low Reynolds number and magnetic Reynolds number problems. Numerical results demonstrating the applicability of this technique are also presented. The simulation of incompressible magneto hydrodynamic fluid flows is illustrated by numerical solution for two-dimensional cases. (c) 2007 Elsevier B.V. All rights reserved.

The thermohydrodynamical picture of Brownian motion via a generalized Smoluchowsky equation

Via an operator continued fraction scheme, we expand Kramers equation in the high friction limit. Then all distribution moments are expressed in terms of the first momemt (particle density). The latter satisfies a generalized Smoluchowsky equation. As an application, we present the nonequilibrium thermodynamics and hydrodynamical picture for the one-dimensional Brownian motion. (C) 2000 Elsevier B.V. B.V. All rights reserved.

Estimating the aboveground biomass and wood volume of savanna woodlands in Brazil's Pantanal wetlands based on allometric correlations

The quantity and distribution of vegetal biomass are important aspects to consider in ecosystem studies. However, little information is available about Brazil's Pantanal woodland savannas. This work involved the development of regression equations of the aerial biomass and wood volume of native tree species in a region of woodland savanna on Rio Negro farm in the Pantanal of Nhecolandia, Brazil. Samples were taken from 10 trees of each of five species: Protium heptaphyllum (Aub1.) Marchand, Magonia pubescens A. St.-Hil., Diptychandra aurantiaca Tul., Terminalia argentea Mart. and Zucc. and Licania minutiflora (Sagot) Fritsch and from a miscellaneous group of I I different species. Linear and nonlinear regression analyses were developed relating the diameter at breast height to the dry weight of wood, branches and leaves, wood volume and total aerial biomass. All the regressions showed a significance of P < 0.05 and an R-2 close to or above 0.8. The biomass curve predicted by linear regression analysis of the studied species was similar to the nonlinear regression, with the exception of L. minutiflora and the miscellaneous group. The breast height diameter proved a good choice for estimating biomass and wood volume. The estimated wood volume and biomass of the Pantanal woodland savanna is crucial information for understanding the carbon cycle and for ensuring the region's conservation and sustainable use. (c) 2006 Elsevier B.V. All rights reserved.

GENERALIZED HEAVISIDE FUNCTIONS IN THE COLOMBEAU THEORY CONTEXT

We defined generalized Heaviside functions for a variable x in R-n, and for variables (x, t) in R-n x R-m. Then study properties such as: composition, invertibility, and association relation (the weak equality). This work is developed in the Colombeau generalized functions context.

Tau-functions and dressing transformations for zero-curvature affine integrable equations

The solutions of a large class of hierarchies of zero-curvature equations that includes Toda- and KdV-type hierarchies are investigated. All these hierarchies are constructed from affine (twisted or untwisted) Kac-Moody algebras g. Their common feature is that they have some special vacuum solutions corresponding to Lax operators lying in some Abelian (up to the central term) subalgebra of g; in some interesting cases such subalgebras are of the Heisenberg type. Using the dressing transformation method, the solutions in the orbit of those vacuum solutions are constructed in a uniform way. Then, the generalized tau-functions for those hierarchies are defined as an alternative set of variables corresponding to certain matrix elements evaluated in the integrable highest-weight representations of g. Such definition of tau-functions applies for any level of the representation, and it is independent of its realization (vertex operator or not). The particular important cases of generalized mKdV and KdV hierarchies as well as the Abelian and non-Abelian affine Toda theories are discussed in detail. © 1997 American Institute of Physics.

Generalized ladder operators for shape-invariant potentials

A general form for ladder operators is used to construct a method to solve bound-state Schrödinger equations. The characteristics of supersymmetry and shape invariance of the system are the start point of the approach. To show the elegance and the utility of the method we use it to obtain energy spectra and eigenfunctions for the one-dimensional harmonic oscillator and Morse potentials and for the radial harmonic oscillator and Coulomb potentials.

Charged Brownian particles: Kramers and Smoluchowski equations and the hydrothermodynamical picture

We consider a charged Brownian gas under the influence of external and non-uniform electric, magnetic and mechanical fields, immersed in a non-uniform bath temperature. With the collision time as an expansion parameter, we study the solution to the associated Kramers equation, including a linear reactive term. To the first order we obtain the asymptotic (overdamped) regime, governed by transport equations, namely: for the particle density, a Smoluchowski- reactive like equation; for the particle's momentum density, a generalized Ohm's-like equation; and for the particle's energy density, a MaxwellCattaneo-like equation. Defining a nonequilibrium temperature as the mean kinetic energy density, and introducing Boltzmann's entropy density via the one particle distribution function, we present a complete thermohydrodynamical picture for a charged Brownian gas. We probe the validity of the local equilibrium approximation, Onsager relations, variational principles associated to the entropy production, and apply our results to: carrier transport in semiconductors, hot carriers and Brownian motors. Finally, we outline a method to incorporate non-linear reactive kinetics and a mean field approach to interacting Brownian particles. © 2011 Elsevier B.V. All rights reserved.

Development of sustainable precision farming systems for swine: Estimating realtime individual amino acid requirements in growing-finishing pigs

The objective of this study was to develop and evaluate a mathematical model used to estimate the daily amino acid requirements of individual growing-finishing pigs. The model includes empirical and mechanistic model components. The empirical component estimates daily feed intake (DFI), BW, and daily gain (DG) based on individual pig information collected in real time. Based on DFI, BW, and DG estimates, the mechanistic component uses classic factorial equations to estimate the optimal concentration of amino acids that must be offered to each pig to meet its requirements. The model was evaluated with data from a study that investigated the effect of feeding pigs with a 3-phase or daily multiphase system. The DFI and BW values measured in this study were compared with those estimated by the empirical component of the model. The coherence of the values estimated by the mechanistic component was evaluated by analyzing if it followed a normal pattern of requirements. Lastly, the proposed model was evaluated by comparing its estimates with those generated by the existing growth model (InraPorc). The precision of the proposed model and InraPorc in estimating DFI and BW was evaluated through the mean absolute error. The empirical component results indicated that the DFI and BW trajectories of individual pigs fed ad libitum could be predicted 1 d (DFI) or 7 d (BW) ahead with the average mean absolute error of 12.45 and 1.85%, respectively. The average mean absolute error obtained with the InraPorc for the average individual of the population was 14.72% for DFI and 5.38% for BW. Major differences were observed when estimates from InraPorc were compared with individual observations. The proposed model, however, was effective in tracking the change in DFI and BW for each individual pig. The mechanistic model component estimated the optimal standardized ileal digestible Lys to NE ratio with reasonable between animal (average CV = 7%) and overtime (average CV = 14%) variation. Thus, the amino acid requirements estimated by model are animal- and time-dependent and follow, in real time, the individual DFI and BW growth patterns. The proposed model can follow the average feed intake and feed weight trajectory of each individual pig in real time with good accuracy. Based on these trajectories and using classical factorial equations, the model makes it possible to estimate dynamically the AA requirements of each animal, taking into account the intake and growth changes of the animal. © 2012 American Society of Animal Science. All rights reserved.

Development of equations to estimate microbial contamination in ruminal incubation residues of forage produced under tropical conditions using 15N as a label1

The objective of this study was to use 15N to label microbial cells to allow development of equations for estimating the microbial contamination in ruminal in situ incubation residues of forage produced under tropical conditions. A total of 24 tropical forages were ruminal incubated in 3 steers at 3 separate times. To determine microbial contamination of the incubated residues, ruminal bacteria were labeled with 15N by continuous intraruminal infusion 60 h before the first incubation and continued until the last day of incubation. Ruminal digesta was collected for the isolation of bacteria before the first infusion of 15N on adaptation period and after the infusion of 15N on collection period. To determine the microbial contamination of CP fractions, restricted models were compared with the full model using the model identity test. A value of the corrected fraction A was estimated from the corresponding noncorrected fraction by this equation: Corrected A fraction (ACPC) = 1.99286 + 0.98256 × A fraction without correction (ACPWC). The corrected fraction B was estimated from the corresponding noncorrected fraction and from CP, NDF, neutral detergent insoluble protein (NDIP), and indigestible NDF (iNDF) using the equation corrected B fraction (BCPC) = -17.2181 - 0.0344 × fraction B without correction (BCPWC) + 0.65433 × CP + 1.03787 × NDF + 2.66010 × NDIP - 0.85979 × iNDF. The corrected degradation rate of B fraction (kd)was estimated using the equation corrected degradation rate of B fraction (kdCPC) = 0.04667 + 0.35139 × degradation rate of B fraction without correction (kdCPWC) + 0.0020 × CP - 0.00055839 × NDF - 0.00336 × NDIP + 0.00075089 × iNDF. This equation was obtained to estimate the contamination using CP of the feeds: %C = 79.21 × (1 - e-0.0555t) × e-0.0874CP. It was concluded that A and B fractions and kd of CP could be highly biased by microbial CP contamination, and therefore these corrected values could be obtained mathematically, replacing the use of microbial markers. The percentage of contamination and the corrected apparent degradability of CP could be obtained from values of CP and time of incubation for each feed, which could reduce cost and labor involved when using 15N. © 2013 American Society of Animal Science. All rights reserved.

On exponential stability of functional differential equations with variable impulse perturbations

We consider a class of functional differential equations subject to perturbations, which vary in time, and we study the exponential stability of solutions of these equations using the theory of generalized ordinary differential equations and Lyapunov functionals. We introduce the concept of variational exponential stability for generalized ordinary differential equations and we develop the theory in this direction by establishing conditions for the trivial solutions of generalized ordinary differential equations to be exponentially stable. Then, we apply the results to get corresponding ones for impulsive functional differential equations. We also present an example of a delay differential equation with Perron integrable right-hand side where we apply our result.

On partial derivative-Neumann's problem in generalized functions framework

In this paper, we show that the equation delta u/delta (z) over bar + Gu = f, where the elements involved are in generalized functions context, has a local solution in the generalized functions context.