A collocation method for solving singular integro-differential equations

This paper describes a collocation method for numerically solving Cauchy-type linear singular integro-differential equations. The numerical method is based on the transformation of the integro-differential equation into an integral equation, and then applying a collocation method to solve the latter. The collocation points are chosen as the Chebyshev nodes. Uniform convergence of the resulting method is then discussed. Numerical examples are presented and solved by the numerical techniques.

A singular integro-differential equation model for dryout in LMFBR boiler tubes

A 2D steady model for the annular two-phase flow of water and steam in the steam-generating boiler pipes of a liquid metal fast breeder reactor is proposed The model is based on thin-layer lubrication theory and thin aerofoil theory. The exchange of mass between the vapour core and the liquid film due to evaporation of the liquid film is accounted for using some simple thermodynamics models, and the resultant change of phase is modelled by proposing a suitable Stefan problem Appropriate boundary conditions for the now are discussed The resulting non-lineal singular integro-differential equation for the shape of the liquid film free surface is solved both asymptotically and numerically (using some regularization techniques) Predictions for the length to the dryout point from the entry of the annular regime are made The influence of both the traction tau provided by the fast-flowing vapour core on the liquid layer and the mass transfer parameter eta on the dryout length is investigated

Non-linear boundary element formulation applied to contact analysis using tangent operator

This work presents a non-linear boundary element formulation applied to analysis of contact problems. The boundary element method (BEM) is known as a robust and accurate numerical technique to handle this type of problem, because the contact among the solids occurs along their boundaries. The proposed non-linear formulation is based on the use of singular or hyper-singular integral equations by BEM, for multi-region contact. When the contact occurs between crack surfaces, the formulation adopted is the dual version of BEM, in which singular and hyper-singular integral equations are defined along the opposite sides of the contact boundaries. The structural non-linear behaviour on the contact is considered using Coulomb`s friction law. The non-linear formulation is based on the tangent operator in which one uses the derivate of the set of algebraic equations to construct the corrections for the non-linear process. This implicit formulation has shown accurate as the classical approach, however, it is faster to compute the solution. Examples of simple and multi-region contact problems are shown to illustrate the applicability of the proposed scheme. (C) 2011 Elsevier Ltd. All rights reserved.

Dual boundary element formulation applied to analysis of multi-fractured domains

This paper deals with analysis of multiple random crack propagation in two-dimensional domains using the boundary element method (BEM). BEM is known to be a robust and accurate numerical technique for analysing this type of problem. The formulation adopted in this work is based on the dual BEM, for which singular and hyper-singular integral equations are used. We propose an iterative scheme to predict the crack growth path and the crack length increment at each time step. The proposed scheme able us to simulate localisation and coalescence phenomena, which is the main contribution of this paper. Considering the fracture mechanics analysis, the displacement correlation technique is applied to evaluate the stress intensity factors. The propagation angle and the equivalent stress intensity factor are calculated using the theory of maximum circumferential stress. Examples of simple and multi-fractured domains, loaded up to the rupture, are considered to illustrate the applicability of the proposed scheme. (C) 2010 Elsevier Ltd. All rights reserved.

A BEM model applied to failure analysis of multi-fractured structures

Due to manufacturing or damage process, brittle materials present a large number of micro-cracks which are randomly distributed. The lifetime of these materials is governed by crack propagation under the applied mechanical and thermal loadings. In order to deal with these kinds of materials, the present work develops a boundary element method (BEM) model allowing for the analysis of multiple random crack propagation in plane structures. The adopted formulation is based on the dual BEM, for which singular and hyper-singular integral equations are used. An iterative scheme to predict the crack growth path and crack length increment is proposed. This scheme enables us to simulate the localization and coalescence phenomena, which are the main contribution of this paper. Considering the fracture mechanics approach, the displacement correlation technique is applied to evaluate the stress intensity factors. The propagation angle and the equivalent stress intensity factor are calculated using the theory of maximum circumferential stress. Examples of multi-fractured domains, loaded up to rupture, are considered to illustrate the applicability of the proposed method. (C) 2011 Elsevier Ltd. All rights reserved.

Multiple random crack propagation using a boundary element formulation

This paper proposes a boundary element method (BEM) model that is used for the analysis of multiple random crack growth by considering linear elastic fracture mechanics problems and structures subjected to fatigue. The formulation presented in this paper is based on the dual boundary element method, in which singular and hyper-singular integral equations are used. This technique avoids singularities of the resulting algebraic system of equations, despite the fact that the collocation points coincide for the two opposite crack faces. In fracture mechanics analyses, the displacement correlation technique is applied to evaluate stress intensity factors. The maximum circumferential stress theory is used to evaluate the propagation angle and the effective stress intensity factor. The fatigue model uses Paris` law to predict structural life. Examples of simple and multi-fractured structures loaded until rupture are considered. These analyses demonstrate the robustness of the proposed model. In addition, the results indicate that this formulation is accurate and can model localisation and coalescence phenomena. (C) 2010 Elsevier Ltd. All rights reserved.

Equifocality of a singular Riemannian foliation

A singular foliation on a complete Riemannian manifold M is said to be Riemannian if each geodesic that is perpendicular to a leaf at one point remains perpendicular to every leaf it meets. We prove that the regular leaves are equifocal, i.e., the end point map of a normal foliated vector field has constant rank. This implies that we can reconstruct the singular foliation by taking all parallel submanifolds of a regular leaf with trivial holonomy. In addition, the end point map of a normal foliated vector field on a leaf with trivial holonomy is a covering map. These results generalize previous results of the authors on singular Riemannian foliations with sections.

BEM formulation for von Karman plates

This work deals with nonlinear geometric plates in the context of von Karman`s theory. The formulation is written such that only the boundary in-plane displacement and deflection integral equations for boundary collocations are required. At internal points, only out-of-plane rotation, curvature and in-plane internal force representations are used. Thus, only integral representations of these values are derived. The nonlinear system of equations is derived by approximating all densities in the domain integrals as single values, which therefore reduces the computational effort needed to evaluate the domain value influences. Hyper-singular equations are avoided by approximating the domain values using only internal nodes. The solution is obtained using a Newton scheme for which a consistent tangent operator was derived. (C) 2009 Elsevier Ltd. All rights reserved.

Integral Piezoactuator System with Optimum Placement of Functionally Graded Material - A Topology Optimization Paradigm

Piezoactuators consist of compliant mechanisms actuated by two or more piezoceramic devices. During the assembling process, such flexible structures are usually bonded to the piezoceramics. The thin bonding layer(s) between the compliant mechanism and the piezoceramic may induce undesirable behavior, including unusual interfacial nonlinearities. This constitutes a drawback of piezoelectric actuators and, in some applications, such as those associated to vibration control and structural health monitoring (e. g., aircraft industry), their use may become either unfeasible or at least limited. A possible solution to this standing problem can be achieved through the functionally graded material concept and consists of developing `integral piezoactuators`, that is those with no bonding layer(s) and whose performance can be improved by tailoring their structural topology and material gradation. Thus, a topology optimization formulation is developed, which allows simultaneous distribution of void and functionally graded piezoelectric materials (including both piezo and non-piezoelectric materials) in the design domain in order to achieve certain specified actuation movements. Two concurrent design problems are considered, that is the optimum design of the piezoceramic property gradation, and the design of the functionally graded structural topology. Two-dimensional piezoactuator designs are investigated because the applications of interest consist of planar devices. Moreover, material gradation is considered in only one direction in order to account for manufacturability issues. To broaden the range of such devices in the field of smart structures, the design of integral Moonie-type functionally graded piezoactuators is provided according to specified performance requirements.

Singular Perturbation for the Discounted Continuous Control of Piecewise Deterministic Markov Processes

This paper deals with the expected discounted continuous control of piecewise deterministic Markov processes (PDMP`s) using a singular perturbation approach for dealing with rapidly oscillating parameters. The state space of the PDMP is written as the product of a finite set and a subset of the Euclidean space a""e (n) . The discrete part of the state, called the regime, characterizes the mode of operation of the physical system under consideration, and is supposed to have a fast (associated to a small parameter epsilon > 0) and a slow behavior. By using a similar approach as developed in Yin and Zhang (Continuous-Time Markov Chains and Applications: A Singular Perturbation Approach, Applications of Mathematics, vol. 37, Springer, New York, 1998, Chaps. 1 and 3) the idea in this paper is to reduce the number of regimes by considering an averaged model in which the regimes within the same class are aggregated through the quasi-stationary distribution so that the different states in this class are replaced by a single one. The main goal is to show that the value function of the control problem for the system driven by the perturbed Markov chain converges to the value function of this limit control problem as epsilon goes to zero. This convergence is obtained by, roughly speaking, showing that the infimum and supremum limits of the value functions satisfy two optimality inequalities as epsilon goes to zero. This enables us to show the result by invoking a uniqueness argument, without needing any kind of Lipschitz continuity condition.

Near-infrared integral field spectroscopy of the Homunculus nebula around eta Carinae using Gemini/CIRPASS

This work presents the first integral field spectroscopy of the Homunculus nebula around eta Carinae in the near-infrared spectral region (J band). We confirmed the presence of a hole on the polar region of each lobe, as indicated by previous near-IR long-slit spectra and mid-IR images. The holes can be described as a cylinder of height (i.e. the thickness of the lobe) and diameter of 6.5 and 6.0 x 10(16) cm, respectively. We also mapped the blue-shifted component of He I lambda 10830 seen towards the NW lobe. Contrary to previous works, we suggested that this blue-shifted component is not related to the Paddle but it is indeed in the equatorial disc. We confirmed the claim of N. Smith and showed that the spatial extent of the Little Homunculus matches remarkably well the radio continuum emission at 3 cm, indicating that the Little Homunculus can be regarded as a small H II region. Therefore, we used the optically thin 1.3 mm radio flux to derive a lower limit for the number of Lyman-continuum photons of the central source in eta Car. In the context of a binary system, and assuming that the ionizing flux comes entirely from the hot companion star, the lower limit for its spectral type and luminosity class ranges from O5.5 III to O7 I. Moreover, we showed that the radio peak at 1.7 arcsec NW from the central star is in the same line-of-sight of the `Sr-filament` but they are obviously spatially separated, while the blue-shifted component of He I lambda 10830 may be related to the radio peak and can be explained by the ultraviolet radiation from the companion star.

Phase portraits of quadratic polynomial vector fields having a rational first integral of degree 3

In this paper, we classify all the global phase portraits of the quadratic polynomial vector fields having a rational first integral of degree 3. (C) 2008 Elsevier Ltd. All rights reserved.

Measurement of thermal neutron cross section and resonance integral for the (165)Ho(n, gamma)(166g)Ho reaction

The ground state thermal neutron cross section and the resonance integral for the (165)Ho(n, gamma)(166)Ho reaction in thermal and 1/E regions, respectively, of a thermal reactor neutron spectrum have been measured experimentally by activation technique. The reaction product, (166)Ho in the ground state, is gaining considerable importance as a therapeutic radionuclide and precisely measured data of the reaction are of significance from the fundamental point of view as well as for application. In this work, the spectrographically pure holmium oxide (Ho(2)O(3)) powder samples were irradiated with and without cadmium covers at the IEA-RI reactor (IPEN, Sao Paulo), Brazil. The deviation of the neutron spectrum shape from 1/E law was measured by co-irradiating Co, Zn, Zr and Au activation detectors with thermal and epithermal neutrons followed by regression and iterative procedures. The magnitudes of the discrepancies that can occur in measurements made with the ideal 1/E law considerations in the epithermal range were studied. The measured thermal neutron cross section at the Maxwellian averaged thermal energy of 0.0253 eV is 59.0 +/- 2.1 b and for the resonance integral 657 +/- 36b. The results are measured with good precision and indicated a consistency trend to resolve the discrepant status of the literature data. The results are compared with the values in main libraries such as ENDF/B-VII, JEF-2.2 and JENDL-3.2, and with other measurements in the literature.

Pseudoclassical description of scalar particle in non-abelian background and path-integral representations

Path-integral representations for a scalar particle propagator in non-Abelian external backgrounds are derived. To this aim, we generalize the procedure proposed by Gitman and Schvartsman of path-integral construction to any representation of SU(N) given in terms of antisymmetric generators. And for arbitrary representations of SU(N), we present an alternative construction by means of fermionic coherent states. From the path-integral representations we derive pseudoclassical actions for a scalar particle placed in non-Abelian backgrounds. These actions are classically analyzed and then quantized to prove their consistency.