Treatment of singularities in 2-D domains using BIEM

The singularities which arise when there is a sudden change of boundary conditions are modelled using spectral shape interpolation functions. The procedure can be used for elasticity as well as potential theory and to any degree of accuracy with respect to the smooth part of the curve.

SERBA: a B.I.E. program with linear elements for 2-D elastostatics analysis

The great developments that have occurred during the last few years in the finite element method and its applications has kept hidden other options for computation. The boundary integral element method now appears as a valid alternative and, in certain cases, has significant advantages. This method deals only with the boundary of the domain, while the F.E.M. analyses the whole domain. This has the following advantages: the dimensions of the problem to be studied are reduced by one, consequently simplifying the system of equations and preparation of input data. It is also possible to analyse infinite domains without discretization errors. These simplifications have the drawbacks of having to solve a full and non-symmetric matrix and some difficulties are incurred in the imposition of boundary conditions when complicated variations of the function over the boundary are assumed. In this paper a practical treatment of these problems, in particular boundary conditions imposition, has been carried out using the computer program shown below. Program SERBA solves general elastostatics problems in 2-dimensional continua using the boundary integral equation method. The boundary of the domain is discretized by line or elements over which the functions are assumed to vary linearly. Data (stresses and/or displacements) are introduced in the local co-ordinate system (element co-ordinates). Resulting stresses are obtained in local co-ordinates and displacements in a general system. The program has been written in Fortran ASCII and implemented on a 1108 Univac Computer. For 100 elements the core requirements are about 40 Kwords. Also available is a Fortran IV version (3 segments)implemented on a 21 MX Hewlett-Packard computer,using 15 Kwords.

Kinematic analysis of a novel 2-d.o.f. orientation device

This paper presents the development of a new parallel robot designed for helping with bone milling surgeries. The robot is a small modular wrist with 2 active degrees of freedom, and it is proposed to be used as an orientation device located at the end of a robotic arm designed for bone milling processes. A generic kinematic geometry is proposed for this device. This first article shows the developments on the workspace optimization and the analysis of the force field required to complete a reconstruction of the inferior jawbone. The singularities of the mechanism are analyzed, and the actuator selection is justified with the torque requirements and the study of the force space. The results obtained by the simulations allow building a first prototype using linear motors. Bone milling experiment video is shown as additional material.

An efficient nonlinear transformation for the numerical computation of the singular integrals appearing in the 2-D Boundary Element Method

We discuss several methods, based on coordinate transformations, for the evaluation of singular and quasisingular integrals in the direct Boundary Element Method. An intrinsec error of some of these methods is detected. Two new transformations are suggested which improve on those currently available.

Cálculo de estructuras continuas: V. 1. Métodos generales : estructuras 1-D y 2-D

Estas notas, que se publican en dos volúmenes, surgen como resultado de las clases de doctorado que, durante el curso académico 1979-80, impartió Avelino Samartín en la Escuela Técnica Superior de Ingenieros de Caminos, Canales y Puertos de la Universidad de Santander. El objetivo del curso era ambicioso ya que se intentaban presentar, de un modo panorámico, los distintos procedimientos analíticos y semianalíticos de Cálculo de Estructuras así como su aplicación al análisis efectivo de diferentes casos prácticos: vigas, placas y láminas. La publicación que se presenta es el volumen uno, que corresponde a la primera parte del curso, referente a los métodos generales de cálculo y su aplicación a casos sencillos : estructuras 1-D y 2-D. También se incluye, en el capítulo 5, una combinación natural de elementos 2-D - laja y placa-, correspondiente a la lámina plegada. Se deja para un segundo volumen el estudio de otras estructuras laminares específicas, así como la introducción a la teoría general de láminas. A lo largo del curso -y en estas notas se refleja- se ha intentado mostrar la identidad de la metodología existente en todo el cálculo estructural , desde la estructura más simple -columna- a la más compleja -lámina-, tanto en su planteamiento como en su resolución. Con esta publicación se desea -aparte de cubrir un evidente vacío bibliográfico, en nuestra lengua, sobre el tema- impulsar el uso de los diferentes métodos analíticos que actualmente existen y que, por causas bien conocidas - facilidad de computación - han sido, en opinión de los autores, inmerecidamente relegados al olvido, cuando podrían, muchas veces, sustituir , o al menos complementar, otras técnicas numéricas más potentes.

Analysis, reliability and optimization of reinforced concrete 2-D structures

In this paper a consistent analysis of reinforced concrete (RC) two-dimensional (2-D) structures,namely slab structures subjected to in-plane and out-plane forces, is presented. By using this method of analysis the well established methodology for dimensioning and verifying RC sections of beam structures is extended to 2-D structures. The validity of the proposed analysis results is checked by comparing them with some published experimental test results. Several examples show some of these proposed analysis features, such as the influence of the reinforcement layout on the service and ultimate behavior of a slab structure and the non straightforward problem of the optimal dimension at a slab point subjected to several loading cases. Also, in these examples, the method applications to design situations as multiple steel families and non orthogonal reinforcement layout are commented.

Simulation of plasma flows in divergent magnetic nozzles

Images from the simulation code DIMAGNO illustrate the roles of pressure, electric, and magnetic forces in the 2-D plasma expansion in a magnetic nozzle and the generation of thrust

Revisiting the Classics to recover the Physical Sense in electrical noise

This paper shows a physically cogent model for electrical noise in resistors that has been obtained from Thermodynamical reasons. This new model derived from the works of Johnson and Nyquist also agrees with the Quantum model for noisy systems handled by Callen and Welton in 1951, thus unifying these two Physical viewpoints. This new model is a Complex or 2-D noise model based on an Admittance that considers both Fluctuation and Dissipation of electrical energy to excel the Real or 1-D model in use that only considers Dissipation. By the two orthogonal currents linked with a common voltage noise by an Admittance function, the new model is shown in frequency domain. Its use in time domain allows to see the pitfall behind a paradox of Statistical Mechanics about systems considered as energy-conserving and deterministic on the microscale that are dissipative and unpredictable on the macroscale and also shows how to use properly the Fluctuation-Dissipation Theorem.

Toward Multi-Scale Object-Based Data Fusion

This paper proposes a new methodology for object based 2-D data fu- sion, with a multiscale character. This methodology is intended to be use in agriculture, specifically in the characterization of the water status of different crops, so as to have an appropriate water management at a farm-holding scale. As a first approach to its evaluation, vegetation cover vigor data has been integrated with texture data. For this purpose, NDVI maps have been calculated using a multispectral image and Lacunarity maps from the panchromatic image. Preliminary results show this methodology is viable in the integration and management of large volumes of data, which characterize the behavior of agricultural covers at farm-holding scale.

Distributed bees algorithm for task allocation in swarm of robots

In this paper, we propose the distributed bees algorithm (DBA) for task allocation in a swarm of robots. In the proposed scenario, task allocation consists in assigning the robots to the found targets in a 2-D arena. The expected distribution is obtained from the targets' qualities that are represented as scalar values. Decision-making mechanism is distributed and robots autonomously choose their assignments taking into account targets' qualities and distances. We tested the scalability of the proposed DBA algorithm in terms of number of robots and number of targets. For that, the experiments were performed in the simulator for various sets of parameters, including number of robots, number of targets, and targets' utilities. Control parameters inherent to DBA were tuned to test how they affect the final robot distribution. The simulation results show that by increasing the robot swarm size, the distribution error decreased.

Improved Integration Methods for P-Adaptive Boundary Element Techniques

Non linear transformations are a good alternative for the numerical evaluation of singular and quasisingular integrals appearing in Boundary Element Method specially in the p-adaptive version. Some aspects of its numerical implementation in 2-D Potential codes is discussed and some examples are shown.

CO2 soil flux baseline at the technological development plant for CO2 injection at Hontomin (Burgos, Spain)

From the end of 2013 and during the following two years, 20 kt of CO2sc are planned to be injected in a saline reservoir (1500 m depth) at the Hontomín site (NE Spain). The target aquifers are Lower Jurassic limestone formations which are sealed by Lower Cretaceous clay units at the Hontomín site (NE Spain). The injection of CO2 is part of the activities committed in the Technology Development phase of the EC-funded OXYCFB300 project (European Energy Program for Recovery – EEPR, http://www.compostillaproject.eu), which include CO2 injection strategies, risk assessment, and testing and validating monitoring methodologies and techniques. Among the monitoring works, the project is intended to prove that present-day technology is able to monitor the evolution of injected CO2 in the reservoir and to detect potential leakage. One of the techniques is the measurement of CO2 flux at the soil–atmosphere interface, which includes campaigns before, during and after the injection operations. In this work soil CO2 flux measurements in the vicinity of oil borehole, drilled in the eighties and named H-1 to H-4, and injection and monitoring wells were performed using an accumulation chamber equipped with an IR sensor. Seven surveys were carried out from November 2009 to summer 2011. More than 4000 measurements were used to determine the baseline flux of CO2 and its seasonal variations. The measured values were low (from 5 to 13 g m−2 day−1) and few outliers were identified, mainly located close to the H-2 oil well. Nevertheless, these values cannot be associated to a deep source of CO2, being more likely related to biological processes, i.e. soil respiration. No anomalies were recognized close to the deep fault system (Ubierna Fault) detected by geophysical investigations. There, the CO2 flux is indeed as low as other measurement stations. CO2 fluxes appear to be controlled by the biological activity since the lowest values were recorded during autumn-winter seasons and they tend to increase in warm periods. Two reference CO2 flux values (UCL50 of 5 g m−2 d−1 for non-ploughed areas in autumn–winter seasons and 3.5 and 12 g m−2 d−1 for in ploughed and non-ploughed areas, respectively, in spring–summer time, and UCL99 of 26 g m−2 d−1 for autumn–winter in not-ploughed areas and 34 and 42 g m−2 d−1 for spring–summer in ploughed and not-ploughed areas, respectively) were calculated. Fluxes higher than these reference values could be indicative of possible leakage during the operational and post-closure stages of the storage project.

B.I.E.M. in solid mechanics

In solid mechanics the weak formulation produces an integral equation ready for a discretization and with less restrictive requiremets than the standard field equations. Fundamentally the weak formulation is a expresion of a green formula. An alternative is to choose another green formula materializing a reciprocity relationship between the basis unknowns and an auxiliary family of functions. The degree of smoothness requiered to practice the discretization is then translated to the auxiliar functions. The subsequent discretization (constant, linear etc.)produces a set of equations on the boundary of the domain. For linear 3-D problems the BIEM appears then as a powerful alternative to FEM, because of the reduction to 2-D thanks to the features previously described.

Impact of N2 plasma power discharge on AlGaN/GaN HEMT performance

The effects of power and time conditions of in situ N2 plasma treatment, prior to silicon nitride (SiN) passivation, were investigated on an AlGaN/GaN high-electron mobility transistor (HEMT). These studies reveal that N2 plasma power is a critical parameter to control the SiN/AlGaN interface quality, which directly affects the 2-D electron gas density. Significant enhancement in the HEMT characteristics was observed by using a low power N2 plasma pretreatment. In contrast, a marked gradual reduction in the maximum drain-source current density (IDS max) and maximum transconductance (gm max), as well as in fT and fmax, was observed as the N2 plasma power increases (up to 40% decrease for 210 W). Different mechanisms were proposed to be dominant as a function of the discharge power range. A good correlation was observed between the device electrical characteristics and the surface assessment by atomic force microscopy and Kelvin force microscopy techniques.

Vertical-cavity semiconductor optical amplifiers (VCSOAs) as optical sensing elements

Semiconductor Optical Amplifiers (SOAs) have mainly found application in optical telecommunication networks for optical signal regeneration, wavelength switching or wavelength conversion. The objective of this paper is to report the use of semiconductor optical amplifiers for optical sensing taking into account their optical bistable properties. As it was previously reported, some semiconductor optical amplifiers, including Fabry-Perot and Distributed-Feedback Semiconductor Optical Amplifiers (FPSOAs and DFBSOAs), may exhibit optical bistability. The characteristics of the attained optical bistability in this kind of devices are strongly dependent on different parameters including wavelength, temperature or applied bias current and small variations lead to a change on their bistable properties. As in previous analyses for Fabry-Perot and DFB SOAs, the variations of these parameters and their possible application for optical sensing are reported in this paper for the case of the Vertical-Cavity Semiconductor Optical Amplifier (VCSOA). When using a VCSOA, the input power needed for the appearance of optical bistability is one order of magnitude lower than that needed in edge-emitting devices. This feature, added to the low manufacturing costs of VCSOAs and the ease to integrate them in 2-D arrays, makes the VCSOA a very promising device for its potential use in optical sensing applications.