A numerical model of a liquid-feed solid polymer electrolyte DMFC and its experimental validation

A one-dimensional, biphasic, multicomponent steady-state model based on phenomenological transport equations for the catalyst layer, diffusion layer, and polymeric electrolyte membrane has been developed for a liquid-feed solid polymer electrolyte direct methanol fuel cell (SPE- DMFC). The model employs three important requisites: (i) implementation of analytical treatment of nonlinear terms to obtain a faster numerical solution as also to render the iterative scheme easier to converge, (ii) an appropriate description of two-phase transport phenomena in the diffusive region of the cell to account for flooding and water condensation/evaporation effects, and (iii) treatment of polarization effects due to methanol crossover. An improved numerical solution has been achieved by coupling analytical integration of kinetics and transport equations in the reaction layer, which explicitly include the effect of concentration and pressure gradient on cell polarization within the bulk catalyst layer. In particular, the integrated kinetic treatment explicitly accounts for the nonhomogeneous porous structure of the catalyst layer and the diffusion of reactants within and between the pores in the cathode. At the anode, the analytical integration of electrode kinetics has been obtained within the assumption of macrohomogeneous electrode porous structure, because methanol transport in a liquid-feed SPE- DMFC is essentially a single-phase process because of the high miscibility of methanol with water and its higher concentration in relation to gaseous reactants. A simple empirical model accounts for the effect of capillary forces on liquid-phase saturation in the diffusion layer. Consequently, diffusive and convective flow equations, comprising Nernst-Plank relation for solutes, Darcy law for liquid water, and Stefan-Maxwell equation for gaseous species, have been modified to include the capillary flow contribution to transport. To understand fully the role of model parameters in simulating the performance of the DMCF, we have carried out its parametric study. An experimental validation of model has also been carried out. (C) 2003 The Electrochemical Society.

Localization induced base isolation in fractionally damped non linear system

In the present article we take up the study of nonlinear localization induced base isolation of a 3 degree of freedom system having cubic nonlinearities under sinusoidal base excitation. The damping forces in the system are described by functions of fractional derivative of the instantaneous displacements, typically linear and quadratic damping are considered here separately. Under the assumption of smallness of certain system parameters and nonlinear terms an approximate estimate of the response at each degree of freedom of the system is obtained by the Method of Multiple Scales approach. We then consider a similar system where the nonlinear terms and certain other parameters are no longer small. Direct numerical simulation is made use of to obtain the amplitude plot in the frequency domain for this case, which helps us to establish the efficacy of this method of base isolation for a broad class of systems. Base isolation obtained this way has no counterpart in the linear theory.

Magnetostatic relations including fluctuation fields - the local expansion

Using the approach of local expansion, we analyze the magnetostatic relations in the case of conventional turbulence. The turbulent relations are obtained consisten tly for themomentum equation and induction equation of both the average and fluctuation relations.In comparison with the magnetostatic relations as discussed usually, turbulent fluctuationfields produce forces, one of which 1/(4π)(α1×B0)×B0 may have parallel and perpendicular components in the direction of magnetic field, the other of which 1/(4π)K×B0 is introduced by the boundary value of turbulence and is perpendicular to the magnetic field. In the case of 2-dimensional configuration of magnetic field, the basic equation will be reduced into a second-order elliptic equation, which includes some linear and nonlinear terms introduced by turbulent fluctuation fields. Turbulent fields may change the configuration of magnetic field and even shear it non-uniformly. The study on the influence of turbulent fields is significant since they are observed in many astrophysical environments.

Obtaining bifurcation diagrams with a thermoacoustic network model

Linear techniques can predict whether the non-oscillating (steady) state of a thermoacoustic system is stable or unstable. With a sufficiently large impulse, however, a thermoacoustic system can reach a stable oscillating state even when the steady state is also stable. A nonlinear analysis is required to predict the existence of this oscillating state. Continuation methods are often used for this but they are computationally expensive. In this paper, an acoustic network code called LOTAN is used to obtain the steady and the oscillating solutions for a horizontal Rijke tube. The heat release is modelled as a nonlinear function of the mass flow rate. Several test cases from the literature are analysed in order to investigate the effect of various nonlinear terms in the flame model. The results agree well with the literature, showing that LOTAN can be used to map the steady and oscillating solutions as a function of the control parameters. Furthermore, the nature of the bifurcation between steady and oscillating states can be predicted directly from the nonlinear terms inside the flame model. Copyright © 2012 by ASME.

THEORETICAL INVESTIGATION OF THE DYNAMIC PROCESS OF THE ILLUMINATION OF GAAS

The dynamic process of light illumination of GaAs is studied numerically in this paper to understand the photoquenching characteristics of the material. This peculiar behavior of GaAs is usally ascribed to the existence of EL2 states and their photodriven metastable states. To understand the conductivity quenching, we have introduced nonlinear terms describing the recombination of the nonequilibrium free electrons and holes into the calculation. Though some photoquenching such as photocapacitance, infrared absorption, and electron-paramagnetic-resonance quenching can be explained qualitatively by only considering the internal transfer between the EL2 state and its metastability, it is essential to take the recombination into consideration for a clear understanding of the photoquenching process. The numerical results and approximate analytical approach are presented in this paper for the first time to our knowledge. The calculation gives quite a reasonable explanation for n-type semiconducting GaAs to have infrared absorption quenching while lacking photoconductance quenching. Also, the calculation results have allowed us to interpret the enhanced photoconductance phenomenon following the conductance quenching in typical semi-insulating GaAs and have shown the expected thermal recovery temperature of about 120 K. The numerical results are in agreement with the reported experiments and have diminished some ambiguities in previous works.

Analytical method of predicating the instabilities of a micro arch-shaped beam under electrostatic loading

An arch-shaped beam with different configurations under electrostatic loading experiences either the direct pull-in instability or the snap-through first and then the pull-in instability. When the pull-in instability occurs, the system collides with the electrode and adheres to it, which usually causes the system failure. When the snap-through instability occurs, the system experiences a discontinuous displacement to flip over without colliding with the electrode. The snap-through instability is an ideal actuation mechanism because of the following reasons: (1) after snap-through the system regains the stability and capability of withstanding further loading; (2) the system flips back when the loading is reduced, i.e. the system can be used repetitively; and (3) when approaching snap-through instability the system effective stiffness reduces toward zero, which leads to a fast flipping-over response. To differentiate these two types of instability responses for an arch-shaped beam is vital for the actuator design. For an arch-shaped beam under electrostatic loading, the nonlinear terms of the mid-plane stretching and the electrostatic loading make the analytical solution extremely difficult if not impossible and the related numerical solution is rather complex. Using the one mode expansion approximation and the truncation of the higher-order terms of the Taylor series, we present an analytical solution here. However, the one mode approximation and the truncation error of the Taylor series can cause serious error in the solution. Therefore, an error-compensating mechanism is also proposed. The analytical results are compared with both the experimental data and the numerical multi-mode analysis. The analytical method presented here offers a simple yet efficient solution approach by retaining good accuracy to analyze the instability of an arch-shaped beam under electrostatic loading.

Dynamic recurrent neural network for system identification and control

A dynamic recurrent neural network (DRNN) that can be viewed as a generalisation of the Hopfield neural network is proposed to identify and control a class of control affine systems. In this approach, the identified network is used in the context of the differential geometric control to synthesise a state feedback that cancels the nonlinear terms of the plant yielding a linear plant which can then be controlled using a standard PID controller.

Dynamic input/output linearization using recurrent neural networks

A dynamic recurrent neural network (DRNN) is used to input/output linearize a control affine system in the globally linearizing control (GLC) structure. The network is trained as a part of a closed loop that involves a PI controller, the goal is to use the network, as a dynamic feedback, to cancel the nonlinear terms of the plant. The stability of the configuration is guarantee if the network and the plant are asymptotically stable and the linearizing input is bounded.

Processos difusivos generalizados

We investigate several diffusion equations which extend the usual one by considering the presence of nonlinear terms or a memory effect on the diffusive term. We also considered a spatial time dependent diffusion coefficient. For these equations we have obtained a new classes of solutions and studied the connection of them with the anomalous diffusion process. We start by considering a nonlinear diffusion equation with a spatial time dependent diffusion coefficient. The solutions obtained for this case generalize the usual one and can be expressed in terms of the q-exponential and q-logarithm functions present in the generalized thermostatistics context (Tsallis formalism). After, a nonlinear external force is considered. For this case the solutions can be also expressed in terms of the q-exponential and q-logarithm functions. However, by a suitable choice of the nonlinear external force, we may have an exponential behavior, suggesting a connection with standard thermostatistics. This fact reveals that these solutions may present an anomalous relaxation process and then, reach an equilibrium state of the kind Boltzmann- Gibbs. Next, we investigate a nonmarkovian linear diffusion equation that presents a kernel leading to the anomalous diffusive process. Particularly, our first choice leads to both a the usual behavior and anomalous behavior obtained through a fractionalderivative equation. The results obtained, within this context, correspond to a change in the waiting-time distribution for jumps in the formalism of random walks. These modifications had direct influence in the solutions, that turned out to be expressed in terms of the Mittag-Leffler or H of Fox functions. In this way, the second moment associated to these distributions led to an anomalous spread of the distribution, in contrast to the usual situation where one finds a linear increase with time

Propriedades físicas de bicamadas magnéticas

The magnetic order of bylayers composed by a ferromagnetic film (F) coupled with an antiferromagnetic film (AF) is studied. Piles of coupled monolayers describe the films and the interfilm coupling is described by an exchange interaction between the magnetic moments at the interface. The F has a cubic anisotropy while the AF has a uniaxial anisotropy. We analyze the effects of an external do magnetic field applied parallel to the interface. We consider the intralayer coupling is strong enough to keep parallel all moments of the monolayer an then they are described by one vector proportional to the magnetization of the layer. The interlayer coupling is represented by an exchange interaction between these vectors. The magnetic energy of the system is the sum of the exchange. Anisotropy and Zeeman energies and the equilibrium configuration is one that gives the absolute minimum of the total energy. The magnetization of the system is calculated and the influence of the external do field combined with the interfilm coupling and the unidirectional anisotropy is studied. Special attention is given to the region near of the transition fields. The torque equation is used to study dynamical behavior of these systems. We consider small oscillations around the equilibrium position and we negleet nonlinear terms to obtain the natural frequencies of the system. The dependence of the frequencies with the external do field and their behavior in the phase transition region is analized

Modelagem Matemática e Computacional de fenômenos eletrocinéticos em meios porosos carregados eletricamente

In this work we present a mathematical and computational modeling of electrokinetic phenomena in electrically charged porous medium. We consider the porous medium composed of three different scales (nanoscopic, microscopic and macroscopic). On the microscopic scale the domain is composed by a porous matrix and a solid phase. The pores are filled with an aqueous phase consisting of ionic solutes fully diluted, and the solid matrix consists of electrically charged particles. Initially we present the mathematical model that governs the electrical double layer in order to quantify the electric potential, electric charge density, ion adsorption and chemical adsorption in nanoscopic scale. Then, we derive the microscopic model, where the adsorption of ions due to the electric double layer and the reactions of protonation/ deprotanaç~ao and zeta potential obtained in modeling nanoscopic arise in microscopic scale through interface conditions in the problem of Stokes and Nerst-Planck equations respectively governing the movement of the aqueous solution and transport of ions. We developed the process of upscaling the problem nano/microscopic using the homogenization technique of periodic structures by deducing the macroscopic model with their respectives cell problems for effective parameters of the macroscopic equations. Considering a clayey porous medium consisting of kaolinite clay plates distributed parallel, we rewrite the macroscopic model in a one-dimensional version. Finally, using a sequential algorithm, we discretize the macroscopic model via the finite element method, along with the interactive method of Picard for the nonlinear terms. Numerical simulations on transient regime with variable pH in one-dimensional case are obtained, aiming computational modeling of the electroremediation process of clay soils contaminated

Stability and collapse of a coupled Bose-Einstein condensate

The dynamics of a coupled Bose-Einstein condensate involving trapped atoms in two quantum states is studied using the time-dependent Gross-Pitaevskii equation including an interaction which can transform atoms from one state to the other. We find interesting oscillation of the number of atoms in each of the states. For all repulsive interactions, stable condensates are formed. When some of the atomic interactions are attractive, the possibility of collapse is studied by including an absorptive contact interaction and a quartic three-body recombination term. One or both components of the condensate may undergo collapse when one or more of the nonlinear terms are attractive in nature. (C) 2001 Elsevier B.V. B.V. All rights reserved.

Stability and collapse of a hybrid Bose-Einstein condensate of atoms and molecules

The dynamics of stability and collapse of a trapped atomic Bose-Einstein condensate (BEC) coupled to a molecular one is studied using the time-dependent Gross-Pitaevskii (GP) equation including a nonlinear interaction term which can transform two atoms into a molecule and vice versa. We find an interesting oscillation of the number of atoms and molecules for a BEC of fixed mass. This oscillation is a consequence of continuous transformation in the condensate of two atoms into a molecule and vice versa. For the study of collapse an absorptive contact interaction and an imaginary quartic three-body recombination term are included in the GP equation. It is possible to have a collapse of one or both components when one or more of the nonlinear terms in the GP equation are attractive in nature, respectively.

Spontaneous symmetry breaking of Bose-Fermi mixtures in double-well potentials

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

A bouncing ball model with two nonlinearities: a prototype for Fermi acceleration

Some dynamical properties of a bouncing ball model under the presence of an external force modelled by two nonlinear terms are studied. The description of the model is made by the use of a two-dimensional nonlinear measure-preserving map on the variable's velocity of the particle and time. We show that raising the straight of a control parameter which controls one of the nonlinearities, the positive Lyapunov exponent decreases in the average and suffers abrupt changes. We also show that for a specific range of control parameters, the model exhibits the phenomenon of Fermi acceleration. The explanation of both behaviours is given in terms of the shape of the external force and due to a discontinuity of the moving wall's velocity.