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.

Euclidean 4d exact solitons in a Skyrme type model

We introduce a Skyrme type, four-dimensional Euclidean field theory made of a triplet of scalar fields n→, taking values on the sphere S2, and an additional real scalar field φ, which is dynamical only on a three-dimensional surface embedded in R4. Using a special ansatz we reduce the 4d non-linear equations of motion into linear ordinary differential equations, which lead to the construction of an infinite number of exact soliton solutions with vanishing Euclidean action. The theory possesses a mass scale which fixes the size of the solitons in way which differs from Derrick's scaling arguments. The model may be relevant to the study of the low energy limit of pure SU(2) Yang-Mills theory. © 2004 Elsevier B.V. All rights reserved.

Plane waves in noncommutative fluids

We study the dynamics of the noncommutative fluid in the Snyder space perturbatively at the first order in powers of the noncommutative parameter. The linearized noncommutative fluid dynamics is described by a system of coupled linear partial differential equations in which the variables are the fluid density and the fluid potentials. We show that these equations admit a set of solutions that are monochromatic plane waves for the fluid density and two of the potentials and a linear function for the third potential. The energy-momentum tensor of the plane waves is calculated. © 2013 Elsevier B.V.

Restauração de imagens digitais com texturas utilizando técnicas de decomposição e equações diferenciais parciais

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

Estabilidade de equações de diferenças quase lineares

Pós-graduação em Matemática Universitária - IGCE

Amphiphilic planar membranes in ionic equilibrium: a study of pH position-dependent values

The pH values near a planar dissociating membrane are studied under a mean field approximation using the Poisson-Boltzmann equation and its linear form. The equations are solved in planar symmetry with the consideration that the charge density on the dissociating membrane surface results from an equilibrium process with the neighboring electrolyte. Results for the membrane dissociation degree are presented as a function of the electrolyte ionic strength and membrane surface charge density. Our calculations indicate that pH values have an appreciable variation within 2 nm from the membrane. It is shown that the dissociation process is enhanced due to the presence of bivalent ions and that pH values acquire better stability than in an electrolyte containing univalent ions.

Estimativa da área foliar de Panicum maximum usando dimensões lineares do limbo foliar

The objective of this study was to obtain a mathematical equation to estimate the leaf area of Panicum maximum using linear measures of leaf blade. Correlation studies were conducted involving the real leaf area (Sf), the main vein leaf length (C), and the maximum leaf width (L). The linear and geometric equations related to C provided good leaf area estimates. For practical reasons, the use of an equation involving only the C*L product is suggested. Thus, an estimate of P. maximum leaf area can be obtained by the equation Sf = 0.6058 (C*L), with the coefficient of determination R = 0.8586.

Análise da estabilidade do movimento rotacional de satélites artificiais com quatérnions e sob a influência de torques externos

The aim of this work is to analyze the stability of the rotational motion’s artificial satellite using the Routh Hurwitz Algorithm (CRH) and the quaternions to describe the satellite’s attitude. This algorithm allows the investigation of the stability of the motion using the coefficients of the characteristic equation associated with the equation of the rotational motion in the linear form. The equations of the rotational motion are given by the four cinematic equations for the quaternion and the three equations of Euler for the spin velocity’s components. In the Euler equations are included the components of the gravity gradient torque (TGG) and the solar radiation torque (TRS). The TGG is generated by the difference of the Earth gravity force direction and intensity actuating on each satellite mass element and it depends on the mass distribution and the form of the satellite. The TRS is created by changing of the linear momentum, which happens due to the interactions of solar photons with the satellite surface. The equilibrium points are gotten by the equation of rotational motion and the CRH is applied in the linear form of these equations. Simulations are developed for small and medium satellites, but the gotten equilibrium points are not stable by CRH. However, when some of the eigenvalues of the characteristic equation are analyzed, it is found some equilibrium points which can be pointed out as stables for an interval of the time, due to small magnitude of the real part of these eigenvalue

A hybrid method for the solution of a sparse power system matrices on vector computers

This paper describes a methodology for solving a linear system of equations on vector computer. The methodology combines direct and inverse factors. The decomposition and implementation of the direct solution in a CRAY Y-MPZE/232, and the performance results are discussed.

Modelagem matemática: uma abordagem do método gráfico e do método simplex na resolução de problemas de otimização

Pós-graduação em Matemática em Rede Nacional - IBILCE

Miscibility in a degenerate fermionic mixture induced by linear coupling

We consider a one-dimensional mean-field-hydrodynamic model of a two-component degenerate Fermi gas in an external trap, each component representing a spin state of the same atom. We demonstrate that the interconversion between them (linear coupling), imposed by a resonant electromagnetic wave, transforms the immiscible binary gas into a miscible state, if the coupling constant, kappa, exceeds a critical value, kappa(cr). The effect is predicted in a variational approximation, and confirmed by numerical solutions. Unlike the recently studied model of a binary Bose-Einsten condensate with the linear coupling, the components in the immiscible phase of the binary fermion mixture never fill two separated domains with a wall between them, but rather form antilocked (pi-phase-shifted) density waves. Another difference from the bosonic mixture is spontaneous breaking of symmetry between the two components in terms of the numbers of atoms in them, N(1) and N(2). The latter effect is characterized by the parameter nu equivalent to(N(1)-N(2))/(N(1)+N(2)) (only N(1)+N(2) is a conserved quantity), the onset of miscibility at kappa >=kappa(cr) meaning a transition to nu equivalent to 0. At kappa linear coupling, and the developing spatial structure resembles the known Larkin-Ovchinnikov-Fulde-Ferrell states in the Fermi mixture. Dynamical states, when kappa is suddenly switched from zero to a value exceeding kappa(cr), are considered too. In the latter case, the system features oscillatory relaxation to the mixed state.

Classes of exact Klein-Gordon equations with spatially dependent masses: Regularizing the one-dimensional inversely linear potential

In this work we consider the effect of a spatially dependent mass over the solution of the Klein-Gordon equation in 1 + 1 dimensions, particularly the case of inversely linear scalar potential, which usually presents problems of divergence of the ground-state wave function at the origin, and possible nonexistence of the even-parity wave functions. Here we study this problem, showing that for a certain dependence of the mass with respect to the coordinate, this problem disappears. (c) 2006 Elsevier B.V. All rights reserved.

A set of linear equations to calculate the steady-state operation of an electrical distribution system

In this paper, the calculation of the steady-state operation of a radial/meshed electrical distribution system (EDS) through solving a system of linear equations (non-iterative load flow) is presented. The constant power type demand of the EDS is modeled through linear approximations in terms of real and imaginary parts of the voltage taking into account the typical operating conditions of the EDS's. To illustrate the use of the proposed set of linear equations, a linear model for the optimal power flow with distributed generator is presented. Results using some test and real systems show the excellent performance of the proposed methodology when is compared with conventional methods. © 2011 IEEE.