970 resultados para preconditioning convection-diffusion equation matrix equation
Resumo:
This paper addresses the statistical mechanics of ideal polymer chains next to a hard wall. The principal quantity of interest, from which all monomer densities can be calculated, is the partition function, G N(z) , for a chain of N discrete monomers with one end fixed a distance z from the wall. It is well accepted that in the limit of infinite N , G N(z) satisfies the diffusion equation with the Dirichlet boundary condition, G N(0) = 0 , unless the wall possesses a sufficient attraction, in which case the Robin boundary condition, G N(0) = - x G N ′(0) , applies with a positive coefficient, x . Here we investigate the leading N -1/2 correction, D G N(z) . Prior to the adsorption threshold, D G N(z) is found to involve two distinct parts: a Gaussian correction (for z <~Unknown control sequence '\lesssim' aN 1/2 with a model-dependent amplitude, A , and a proximal-layer correction (for z <~Unknown control sequence '\lesssim' a described by a model-dependent function, B(z).
Resumo:
In this paper we conclude the analysis started in [J.M. Arrieta, AN Carvalho, G. Lozada-Cruz, Dynamics in dumbbell domains I. Continuity of the set of equilibria, J. Differential Equations 231 (2006) 551-597] and continued in [J.M. Arrieta, AN Carvalho, G. Lozada-Cruz, Dynamics in dumbbell domains II. The limiting problem, J. Differential Equations 247 (1) (2009) 174-202 (this issue)] concerning the behavior of the asymptotic dynamics of a dissipative reaction-diffusion equation in a dumbbell domain as the channel shrinks to a line segment. In [J.M. Arrieta, AN Carvalho. G. Lozada-Cruz, Dynamics in dumbbell domains I. Continuity of the set of equilibria, J. Differential Equations 231 (2006) 551-597], we have established an appropriate functional analytic framework to address this problem and we have shown the continuity of the set of equilibria. In [J.M. Arrieta, AN Carvalho, G. Lozada-Cruz. Dynamics in dumbbell domains II. The limiting problem, J. Differential Equations 247 (1) (2009) 174-202 (this issue)], we have analyzed the behavior of the limiting problem. In this paper we show that the attractors are Upper semicontinuous and, moreover, if all equilibria of the limiting problem are hyperbolic, then they are lower semicontinuous and therefore, continuous. The continuity is obtained in L(p) and H(1) norms. (C) 2008 Elsevier Inc. All rights reserved.
Resumo:
We consider a 1-dimensional reaction-diffusion equation with nonlinear boundary conditions of logistic type with delay. We deal with non-negative solutions and analyze the stability behavior of its unique positive equilibrium solution, which is given by the constant function u equivalent to 1. We show that if the delay is small, this equilibrium solution is asymptotically stable, similar as in the case without delay. We also show that, as the delay goes to infinity, this equilibrium becomes unstable and undergoes a cascade of Hopf bifurcations. The structure of this cascade will depend on the parameters appearing in the equation. This equation shows some dynamical behavior that differs from the case where the nonlinearity with delay is in the interior of the domain. (C) 2009 Elsevier Inc. All rights reserved.
Resumo:
The goal of this paper is to analyze the character of the first Hopf bifurcation (subcritical versus supercritical) that appears in a one-dimensional reaction-diffusion equation with nonlinear boundary conditions of logistic type with delay. We showed in the previous work [Arrieta et al., 2010] that if the delay is small, the unique non-negative equilibrium solution is asymptotically stable. We also showed that, as the delay increases and crosses certain critical value, this equilibrium becomes unstable and undergoes a Hopf bifurcation. This bifurcation is the first one of a cascade occurring as the delay goes to infinity. The structure of this cascade will depend on the parameters appearing in the equation. In this paper, we show that the first bifurcation that occurs is supercritical, that is, when the parameter is bigger than the delay bifurcation value, stable periodic orbits branch off from the constant equilibrium.
Resumo:
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
Resumo:
The dengue virus is transmitted in regions previously infested with the mosquito Aedes aegypti. To assess the spreading and establishment of the dengue disease vector, a mathematical model is developed that takes into account the diffusion and advection phenomena. A discrete model based on the cellular automata approach, which is a good framework to deal with small populations, is also developed to be compared with the continuous modeling.
Resumo:
A procedure to model optical diffused-channel waveguides is presented in this work. The dielectric waveguides present anisotropic refractive indexes which are calculated from the proton concentration. The proton concentration inside the channel is calculated by the anisotropic 2D-linear diffusion equation and converted to the refractive indexes using mathematical relations obtained from experimental data, the arbitrary refractive index profile is modeled by a. nodal expansion in the base functions. The TE and TM-like propagation properties (effective index) and the electromagnetic fields for well-annealed proton-exchanged (APE) LiNbO3 waveguides are computed by the finite element method.
Resumo:
Fresh persimmon has a high moisture content (about 85% wet basis) making it highly perishable and requiring adequate drying conditions to obtain an acceptable dehydrated product. Drying kinetics of persimmon cv. Rama Forte was studied in a fixed bed dryer at temperatures ranging from 50 to 80 degreesC and air velocity of 0.8 m/s. Shrinkage during drying was described by a linear correlation with respect to water content. Evaluation of effective diffusivity as a function of moisture content, with undergoing shrinkage during drying was based on Fourier series solution of Fick's diffusion equation. Effective diffusivity values at moisture contents between 0.09 - 4.23 kg water/kg dry matter were found to be in the range of 2.6 x 10(-10) m(2)/s to 5.4 x 10(-10) m(2)/s, and its dependence on air drying temperature was represented by an Arrhenius type equation. Activation energy increased with decreasing water content in persimmons.
Resumo:
The pressure field in a high-power klystron amplifier is investigated to scale the ionic vacuum pump used to maintain the ultra high-vacuum in the device in order to increase its life-time. The investigation is conducted using an 1.3 GHz, 100 A - 240 keV high-power klystron with five reentrant coaxial cavities, assembled in a cylindrical drift tube 1.2 m long. The diffusion equation is solved to the regime molecular flow to obtain the pressure profile along the axis of the klystron drift tube. The model, solved by both analytical and numerical procedures, is able to determine the pressure values in steady-state case. This work considers the specific conductance and all important gas sources, as in the degassing of the drift tube and cavities walls, cathode, and collector. For the drift tube degassing rate equals to q(deg) = 2x10(-12) (-)mbar.L.s(-1) cm(-2) (degassing rate per unit area), to cavities q(cavity) = 3x10(-13) mbar.L.s(-1)cm(-2), to the cathode q(cathode) = 6x10(-9)_mbar.L.s(-1) and to the collector q(collector) = 6x10(-9) mbar.L.s(-1), it was found that a 10 L.s(-1) ionic vacuum pump connected in the output waveguide wall is suitable. In this case, the pressure obtained in the cathode is p(cathode) = 6.3x10(-9) mbar, in the collector p(collector) = 2.7x10(-9) mbar, and in the output waveguide p = 2.1x10(-9) mbar. Although only the steady-state case is analyzed, some aspects that may be relevant in a transient situation, for instance, when the beam hits the drift tube walls, producing a gas burst, is also commented.
Resumo:
We analyze the dynamics of a reaction-diffusion equation with homogeneous Neumann boundary conditions in a dumbbell domain. We provide an appropriate functional setting to treat this problem and, as a first step, we show in this paper the continuity of the set of equilibria and of its linear unstable manifolds. (c) 2006 Elsevier B.V. All rights reserved.
Resumo:
Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)
Resumo:
The pressure field of a high-power klystron amplifier in the cathode and anode region was investigated. The investigation was performed using a 1.3 GHz, 100 A and 240 kV high-power klystron with five reentrant coaxial cavities, assembled in cylindrical drift tube 1.2 m long. The diffusion equation in mathematical model was also solved by using a 3-D finite element method code, in order to obtain pressure profile in region of interest. The results show that density profile of molecules between cathode-anode region was determined, where cathode pressure is approximately 10% higher than anode pressure.
Resumo:
Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)
Resumo:
Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)
Resumo:
Apresentamos aqui uma metodologia alternativa para modelagem de ferramentas de indução diretamente no domínio do tempo. Este trabalho consiste na solução da equação de difusão do campo eletromagnético através do método de diferenças finitas. O nosso modelo consiste de um meio estratificado horizontalmente, através do qual simulamos um deslocamento da ferramenta na direção perpendicular às interfaces. A fonte consiste de uma bobina excitada por uma função degrau de corrente e o registro do campo induzido no meio é feito através de uma bobina receptora localizada acima da bobina transmissora. Na solução da equação de difusão determinamos o campo primário e o campo secundário separadamente. O campo primário é obtido analiticamente e o campo secundário é determinado utilizando-se o método de Direção Alternada Implícita, resultando num sistema tri-diagonal que é resolvido através do método recursivo proposto por Claerbout. Finalmente, determina-se o valor máximo do campo elétrico secundário em cada posição da ferramenta ao longo da formação, obtendo-se assim uma perfilagem no domínio do tempo. Os resultados obtidos mostram que este método é bastante eficiente na determinação do contato entre camadas, inclusive para camadas de pequena espessura.