Patterns in parabolic problems with nonlinear boundary conditions

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

An application of the Element Free Galerkin Method to the analysis of quantum well structures

In order to obtain the quantum-mechanical properties of layered semicondutor structures (quantum well and superlattice structures, for instance), solutions of the Schrodinger equation should be obtained for arbitrary potential profiles. In this paper, it is shown that such problems may be also studied by the Element Free Galerkin Method.

A MODEL OF THE DEUTERON BASED ON A BREIT TYPE EQUATION

A relativistic treatment of the deuteron and its observables based on a two-body Dirac (Breit) equation, with phenomenological interactions, associated to one-boson exchanges with cutoff masses, is presented. The 16-component wave function for the deuteron (J(pi) = 1+) solution contains four independent radial functions which obey a system of four coupled differential equations of first order. This radial system is numerically integrated, from infinity to the origin, by fixing the value of the deuteron binding energy and using appropriate boundary conditions at infinity. Specific examples of mixtures containing scalar, pseudoscalar and vector like terms are discussed in some detail and several observables of the deuteron are calculated. Our treatment differs from more conventional ones in that nonrelativistic reductions of the order c-2 are not used.

A note on the controllability for the wave equation in nonsmooth plane domains

We study exact boundary controllability for a two-dimensional wave equation in a region which is an angular sector of a circle or an angular sector of an annular region. The control, of Neumann type, acts on the curved part of the boundary, while in the straight part we impose homogeneous Dirichlet boundary condition. The initial state has finite energy and the control is square integrable. (c) 2005 Elsevier B.V. All rights reserved.

An implicit BEM formulation to model strong discontinuities in solids

A boundary element method (BEM) formulation to predict the behavior of solids exhibiting displacement (strong) discontinuity is presented. In this formulation, the effects of the displacement jump of a discontinuity interface embedded in an internal cell are reproduced by an equivalent strain field over the cell. To compute the stresses, this equivalent strain field is assumed as the inelastic part of the total strain. As a consequence, the non-linear BEM integral equations that result from the proposed approach are similar to those of the implicit BEM based on initial strains. Since discontinuity interfaces can be introduced inside the cell independently on the cell boundaries, the proposed BEM formulation, combined with a tracking scheme to trace the discontinuity path during the analysis, allows for arbitrary discontinuity propagation using a fixed mesh. A simple technique to track the crack path is outlined. This technique is based on the construction of a polygonal line formed by segments inside the cells, in which the assumed failure criterion is reached. Two experimental concrete fracture tests were analyzed to assess the performance of the proposed formulation.

Linear analysis of building floor structures by a BEM formulation based on Reissner's theory

In this work, the plate bending formulation of the boundary element method (BEM) based on the Reissner's hypothesis is extended to the analysis of zoned plates in order to model a building floor structure. In the proposed formulation each sub-region defines a beam or a slab and depending on the way the sub-regions are represented, one can have two different types of analysis. In the simple bending problem all sub-regions are defined by their middle surface. on the other hand, for the coupled stretching-bending problem all sub-regions are referred to a chosen reference surface, therefore eccentricity effects are taken into account. Equilibrium and compatibility conditions are automatically imposed by the integral equations, which treat this composed structure as a single body. The bending and stretching values defined on the interfaces are approximated along the beam width, reducing therefore the number of degrees of freedom. Then, in the proposed model the set of equations is written in terms of the problem values on the beam axis and on the external boundary without beams. Finally some numerical examples are presented to show the accuracy of the proposed model.

Ferroelectric state analysis in grain boundary of Na0.85Li0.15NbO3 ceramic

The electric and dielectric properties of the grain boundary of Na0.85Li0.15NbO3 lead-free ferroelectric-semiconductor perovskite were investigated. The impedance spectroscopy was carried out as a function of a thermal cycle. The sodium lithium niobate was synthesized by a chemical route based on the evaporation method. Dense ceramic, relative density of 97%, was prepared at 1423 K for 2 h in air atmosphere. ac measurements were carried out in the frequency range of 5 Hz-13 MHz and from 673 to 1023 K. Theoretical adjust of the impedance data was performed to deriving the electric parameters of the grain boundary. The electric conductivity follows the Arrhenius law, with activation energy values equal to 1.55 and 1.54 eV for heating and cooling cycle, respectively. The nonferroelectric state of the grain boundary and its correlation with symmetry are discussed in the temperature domain. (C) 2003 American Institute of Physics.

Evaluating the residual stress in PbTiO3 thin films prepared by a polymeric chemical method

We report a study of residual stress in PbTiO3 (PT) thin films prepared on Si substrates by a polymeric chemical method. The E(1TO) frequency was used to evaluate the residual stress through an empirical equation available for bulk PT. We find that the residual stress in PT films increases as the film thickness decreases and conclude that it originates essentially from the contributions of extrinsic and intrinsic factors. Polarized Raman experiments showed that the PT films prepared by a polymeric chemical method are somewhat a-domain (polar axis c parallel to the substrate) oriented.

An application of the Adomian's decomposition method to the analysis of MHD duct flows

A new approach is proposed in this work for the treatment of boundary value problems through the Adomian's decomposition method. Although frequently claimed as accurate and having fast convergence rates, the original formulation of Adomian's method does not allow the treatment of homogeneous boundary conditions along closed boundaries. The technique here presented overcomes this difficulty, and is applied to the analysis of magnetohydrodynamic duct flows. Results are in good agreement with finite element method calculations and analytical solutions for square ducts. Therefore, new possibilities appear for the application of Adomian's method in electromagnetics.

A well-balanced flow equation for noise removal and edge detection

In this paper, an anisotropic nonlinear diffusion equation for image restoration is presented. The model has two terms: the diffusion and the forcing term. The balance between these terms is made in a selective way, in which boundary points and interior points of the objects that make up the image are treated differently. The optimal smoothing time concept, which allows for finding the ideal stop time for the evolution of the partial differential equation is also proposed. Numerical results show the proposed model's high performance.

Dirac equation in magnetic-solenoid field

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

Multiple-time higher-order perturbation analysis of the regularized long-wavelength equation

By considering the long-wavelength limit of the regularized long wave (RLW) equation, we study its multiple-time higher-order evolution equations. As a first result, the equations of the Korteweg-de Vries hierarchy are shown to play a crucial role in providing a secularity-free perturbation theory in the specific case of a solitary-wave solution. Then, as a consequence, we show that the related perturbative series can be summed and gives exactly the solitary-wave solution of the RLW equation. Finally, some comments and considerations are made on the N-soliton solution, as well as on the limitations of applicability of the multiple-scale method in obtaining uniform perturbative series.

A study about the solutions of Schrödinger equation with an asymptotically linear potential

We determine the solutions of the Schrödinger equation for an asymptotically linear potential. Analytical solutions are obtained by superalgebra in quantum mechanics and we establish when these solutions are possible. Numerical solutions for the spectra are obtained by the shifted 1/N expansion method.

The Role of the Korteweg-de Vries Hierarchy in the N-Soliton Dynamics of the Shallow Water Wave Equation

We apply a multiple-time version of the reductive perturbation method to study long waves as governed by the shallow water wave model equation. As a consequence of the requirement of a secularity-free perturbation theory, we show that the well known N-soliton dynamics of the shallow water wave equation, in the particular case of α = 2β, can be reduced to the N-soliton solution that satisfies simultaneously all equations of the Korteweg-de Vries hierarchy.

Fourth-order method for solving the Navier-Stokes equations in a constricting channel

A fourth-order numerical method for solving the Navier-Stokes equations in streamfunction/vorticity formulation on a two-dimensional non-uniform orthogonal grid has been tested on the fluid flow in a constricted symmetric channel. The family of grids is generated algebraically using a conformal transformation followed by a non-uniform stretching of the mesh cells in which the shape of the channel boundary can vary from a smooth constriction to one which one possesses a very sharp but smooth corner. The generality of the grids allows the use of long channels upstream and downstream as well as having a refined grid near the sharp corner. Derivatives in the governing equations are replaced by fourth-order central differences and the vorticity is eliminated, either before or after the discretization, to form a wide difference molecule for the streamfunction. Extra boundary conditions, necessary for wide-molecule methods, are supplied by a procedure proposed by Henshaw et al. The ensuing set of non-linear equations is solved using Newton iteration. Results have been obtained for Reynolds numbers up to 250 for three constrictions, the first being smooth, the second having a moderately sharp corner and the third with a very sharp corner. Estimates of the error incurred show that the results are very accurate and substantially better than those of the corresponding second-order method. The observed order of the method has been shown to be close to four, demonstrating that the method is genuinely fourth-order. © 1977 John Wiley & Sons, Ltd.