Two-dimensional mesh optimization in the finite element method

The solution to the problem of finding the optimum mesh design in the finite element method with the restriction of a given number of degrees of freedom, is an interesting problem, particularly in the applications method. At present, the usual procedures introduce new degrees of freedom (remeshing) in a given mesh in order to obtain a more adequate one, from the point of view of the calculation results (errors uniformity). However, from the solution of the optimum mesh problem with a specific number of degrees of freedom some useful recommendations and criteria for the mesh construction may be drawn. For 1-D problems, namely for the simple truss and beam elements, analytical solutions have been found and they are given in this paper. For the more complex 2-D problems (plane stress and plane strain) numerical methods to obtain the optimum mesh, based on optimization procedures have to be used. The objective function, used in the minimization process, has been the total potential energy. Some examples are presented. Finally some conclusions and hints about the possible new developments of these techniques are also given.

A novel method for measuring the dynamic fracture-initiation toughness under impulsive loadings

The design and development of a new method for performing fracture toughness tests under impulsive loadings using explosives is presented. The experimental set-up was complemented with pressure transducers and strain gauges in order to measure, respectively, the blast wave that reached the specimen and the loading history. Fracture toughness tests on a 7017-T73 aluminium alloy were carried out by using this device under impulsive loadings. Previous studies reported that such aluminium alloy had very little strain rate sensitivity, which made it an ideal candidate for comparison at different loading rates. The fracture-initiation toughness values of the 7017-T73 aluminium alloy obtained at impulsive loadings did not exhibit a significant variation from the cases studied at lower loading rates. Therefore, the method and device developed for measuring the dynamic fracture-initiation toughness under impulsive loadings was considered suitable for such a purpose

Conditions for perfect focusing multiple point sources with the SMS design method

In this work, we demonstrate how it is possible to sharply image multiple object points. The Simultaneous Multiple Surface (SMS) design method has usually been presented as a method to couple N wave-front pairs with N surfaces, but recent findings show that when using N surfaces, we can obtain M image points when Nmethod, from its basics, to imaging two object points through one surface, the transition from two to three objet points obtained by increasing the parallelism, and getting to the designs of six surfaces imaging up to eight object points. These designs are limited with the condition that the surfaces cannot be placed at the aperture stop. In the process of maximizing the object points to sharp image, we try to exhaust the degrees of freedom of aspherics and free-forms. We conjecture that maximal SMS designs are very close to a good solution, hence using them as a starting point for the optimization will lead us faster to a final optical system. We suggest here different optimization strategies which combined with the SMS method are proven to give the best solution. Through the example of imaging with the high aspect ratio, we compare the results obtained optimizing the rotational lens and using a combination of SMS method and optimization, showing that the second approach is giving significantly smaller value of overall RMS spot diameter.

Diseño de un sistema regional de alerta de tsunamis como parte integral de la operación portuaria

En este trabajo se presenta el desarrollo de una metodología para obtener un universo de funciones de Green y el algoritmo correspondiente, para estimar la altura de tsunamis a lo largo de la costa occidental de México en función del momento sísmico y de la extensión del área de ruptura de sismos interplaca localizados entre la costa y la Trinchera Mesoamericana. Tomando como caso de estudio el sismo ocurrido el 9 de octubre de 1995 en la costa de Jalisco-Colima, se estudiaron los efectos del tsunami originados en la hidrodinámica del Puerto de Manzanillo, México, con una propuesta metodológica que contempló lo siguiente: El primer paso de la metodología contempló la aplicación del método inverso de tsunamis para acotar los parámetros de la fuente sísmica mediante la confección de un universo de funciones de Green para la costa occidental de México. Tanto el momento sísmico como la localización y extensión del área de ruptura de sismos se prescribe en segmentos de planos de falla de 30 X 30 km. A cada uno de estos segmentos del plano de falla corresponde un conjunto de funciones de Green ubicadas en la isobata de 100 m, para 172 localidades a lo largo de la costa, separadas en promedio 12 km entre una y otra. El segundo paso de la metodología contempló el estudio de la hidrodinámica (velocidades de las corrientes y niveles del mar en el interior del puerto y el estudio del runup en la playa) originada por el tsunami, la cual se estudió en un modelo hidráulico de fondo fijo y en un modelo numérico, representando un tsunami sintético en la profundidad de 34 m como condición inicial, el cual se propagó a la costa con una señal de onda solitaria. Como resultado de la hidrodinámica del puerto de Manzanillo, se realizó un análisis de riesgo para la definición de las condiciones operativas del puerto en términos de las velocidades en el interior del mismo, y partiendo de las condiciones iniciales del terremoto de 1995, se definieron las condiciones límites de operación de los barcos en el interior y exterior del puerto. In this work is presented the development of a methodology in order to obtain a universe of Green's functions and the corresponding algorithm in order to estimate the tsunami wave height along the west coast of Mexico, in terms of seismic moment and the extent of the area of the rupture, in the interplate earthquakes located between the coast and the Middle America Trench. Taking as a case of study the earthquake occurred on October 9, 1995 on the coast of Jalisco-Colima, were studied the hydrodynamics effects of the tsunami caused in the Port of Manzanillo, Mexico, with a methodology that contemplated the following The first step of the methodology contemplated the implementation of the tsunami inverse method to narrow the parameters of the seismic source through the creation of a universe of Green's functions for the west coast of Mexico. Both the seismic moment as the location and extent of earthquake rupture area prescribed in segments fault planes of 30 X 30 km. Each of these segments of the fault plane corresponds a set of Green's functions located in the 100 m isobath, to 172 locations along the coast, separated on average 12 km from each other. The second step of the methodology contemplated the study of the hydrodynamics (speed and directions of currents and sea levels within the port and the study of the runup on the beach Las Brisas) caused by the tsunami, which was studied in a hydraulic model of fix bed and in a numerical model, representing a synthetic tsunami in the depth of 34 m as an initial condition which spread to the coast with a solitary wave signal. As a result of the hydrodynamics of the port of Manzanillo, a risk analysis to define the operating conditions of the port in terms of the velocities in the inner and outside of the port was made, taken in account the initial conditions of the earthquake and tsunami ocurred in Manzanillo port in 1995, were defined the limits conditions of operation of the ships inside and outside the port.

Real-time detection of DNA hybridization and melting on oligonucleotide arrays by using optical wave guides.

The challenge of the Human Genome Project is to increase the rate of DNA sequence acquisition by two orders of magnitude to complete sequencing of the human genome by the year 2000. The present work describes a rapid detection method using a two-dimensional optical wave guide that allows measurement of real-time binding or melting of a light-scattering label on a DNA array. A particulate label on the target DNA acts as a light-scattering source when illuminated by the evanescent wave of the wave guide and only the label bound to the surface generates a signal. Imaging/visual examination of the scattered light permits interrogation of the entire array simultaneously. Hybridization specificity is equivalent to that obtained with a conventional system using autoradiography. Wave guide melting curves are consistent with those obtained in the liquid phase and single-base discrimination is facile. Dilution experiments showed an apparent lower limit of detection at 0.4 nM oligonucleotide. This performance is comparable to the best currently known fluorescence-based systems. In addition, wave guide detection allows manipulation of hybridization stringency during detection and thereby reduces DNA chip complexity. It is anticipated that this methodology will provide a powerful tool for diagnostic applications that require rapid cost-effective detection of variations from known sequences.

Electrical investigation of p-type 4H-SiC heavily doped by ion implantation

In the last decades, an increasing interest in the research field of wide bandgap semiconductors was observed, mostly due to the progressive approaching of silicon-based devices to their theoretical limits. 4H-SiC is an example among these, and is a mature compound for applications. The main advantages offered 4H-SiC in comparison with silicon are an higher breakdown field, an higher thermal conductivity, a higher operating temperature, very high hardness and melting point, biocompatibility, but also low switching losses in high frequencies applications and lower on-resistances in unipolar devices. Then, 4H-SiC power devices offer great performance improvement; moreover, they can work in hostile environments where silicon power devices cannot function. Ion implantation technology is a key process in the fabrication of almost all kinds of SiC devices, owing to the advantage of a spatially selective doping. This work is dedicated to the electrical investigation of several differently-processed 4H-SiC ion- implanted samples, mainly through Hall effect and space charge spectroscopy experiments. It was also developed the automatic control (Labview) of several experiments. In the work, the effectiveness of high temperature post-implant thermal treatments (up to 2000°C) were studied and compared considering: (i) different methods, (ii) different temperatures and (iii) different duration of the annealing process. Preliminary p + /n and Schottky junctions were also investigated as simple test devices. 1) Heavy doping by ion implantation of single off-axis 4H-SiC layers The electrical investigation is one of the most important characterization of ion-implanted samples, which must be submitted to mandatory post-implant thermal treatment in order to both (i) recover the lattice after ion bombardment, and (ii) address the implanted impurities into lattice sites so that they can effectively act as dopants. Electrical investigation can give fundamental information on the efficiency of the electrical impurity activation. To understand the results of the research it should be noted that: (a) To realize good ohmic contacts it is necessary to obtain spatially defined highly doped regions, which must have conductivity as low as possible. (b) It has been shown that the electrical activation efficiency and the electrical conductivity increase with the annealing temperature increasing. (c) To maximize the layer conductivity, temperatures around 1700°C are generally used and implantation density high till to 10 21 cm -3 . In this work, an original approach, different from (c), is explored by the using very high annealing temperature, around 2000°C, on samples of Al + -implant concentration of the order of 10 20 cm -3 . Several Al + -implanted 4H-SiC samples, resulting of p-type conductivity, were investigated, with a nominal density varying in the range of about 1-5∙10 20 cm -3 and subjected to two different high temperature thermal treatments. One annealing method uses a radiofrequency heated furnace till to 1950°C (Conventional Annealing, CA), the other exploits a microwave field, providing a fast heating rate up to 2000°C (Micro-Wave Annealing, MWA). In this contest, mainly ion implanted p-type samples were investigated, both off-axis and on-axis <0001> semi-insulating 4H-SiC. Concerning p-type off-axis samples, a high electrical activation of implanted Al (50-70%) and a compensation ratio below 10% were estimated. In the work, the main sample processing parameters have been varied, as the implant temperature, CA annealing duration, and heating/cooling rates, and the best values assessed. MWA method leads to higher hole density and lower mobility than CA in equivalent ion implanted layers, resulting in lower resistivity, probably related to the 50°C higher annealing temperature. An optimal duration of the CA treatment was estimated in about 12-13 minutes. A RT resistivity on the lowest reported in literature for this kind of samples, has been obtained. 2) Low resistivity data: variable range hopping Notwithstanding the heavy p-type doping levels, the carrier density remained less than the critical one required for a semiconductor to metal transition. However, the high carrier densities obtained was enough to trigger a low temperature impurity band (IB) conduction. In the heaviest doped samples, such a conduction mechanism persists till to RT, without significantly prejudice the mobility values. This feature can have an interesting technological fall, because it guarantee a nearly temperature- independent carrier density, it being not affected by freeze-out effects. The usual transport mechanism occurring in the IB conduction is the nearest neighbor hopping: such a regime is effectively consistent with the resistivity temperature behavior of the lowest doped samples. In the heavier doped samples, however, a trend of the resistivity data compatible with a variable range hopping (VRH) conduction has been pointed out, here highlighted for the first time in p-type 4H-SiC. Even more: in the heaviest doped samples, and in particular, in those annealed by MWA, the temperature dependence of the resistivity data is consistent with a reduced dimensionality (2D) of the VRH conduction. In these samples, TEM investigation pointed out faulted dislocation loops in the basal plane, whose average spacing along the c-axis is comparable with the optimal length of the hops in the VRH transport. This result suggested the assignment of such a peculiar behavior to a kind of spatial confinement into a plane of the carrier hops. 3) Test device the p + -n junction In the last part of the work, the electrical properties of 4H-SiC diodes were also studied. In this case, a heavy Al + ion implantation was realized on n-type epilayers, according to the technological process applied for final devices. Good rectification properties was shown from these preliminary devices in their current-voltage characteristics. Admittance spectroscopy and deep level transient spectroscopy measurements showed the presence of electrically active defects other than the dopants ones, induced in the active region of the diodes by ion implantation. A critical comparison with the literature of these defects was performed. Preliminary to such an investigation, it was assessed the experimental set up for the admittance spectroscopy and current-voltage investigation and the automatic control of these measurements.

Excitations in Bose-Einstein condensates: collective modes, quantum turbulence and matter wave statistics

In this thesis, we present the generation and studies of a 87Rb Bose-Einstein condensate (BEC) perturbed by an oscillatory excitation. The atoms are trapped in a harmonic magnetic trap where, after an evaporative cooling process, we produce the BEC. In order to study the effect caused by oscillatory excitations, a quadrupole magnetic field time oscillatory is superimposed to the trapping potential. Through this perturbation, collective modes were observed. The dipole mode is excited even for low excitation amplitudes. However, a minimum excitation energy is needed to excite the condensate quadrupole mode. Observing the excited cloud in TOF expansion, we note that for excitation amplitude in which the quadrupole mode is excited, the cloud expands without invert its aspect ratio. By looking these clouds, after long time-of-flight, it was possible to see vortices and, sometimes, a turbulent state in the condensed cloud. We calculated the momentum distribution of the perturbed BECs and a power law behavior, like the law to Kolmogorov turbulence, was observed. Furthermore, we show that using the method that we have developed to calculate the momentum distribution, the distribution curve (including the power law exponent) exhibits a dependence on the quadrupole mode oscillation of the cloud. The randomness distribution of peaks and depletions in density distribution image of an expanded turbulent BEC, remind us to the intensity profile of a speckle light beam. The analogy between matter-wave speckle and light speckle is justified by showing the similarities in the spatial propagation (or time expansion) of the waves. In addition, the second order correlation function is evaluated and the same dependence with distance was observed for the both waves. This creates the possibility to understand the properties of quantum matter in a disordered state. The propagation of a three-dimensional speckle field (as the matter-wave speckle described here) creates an opportunity to investigate the speckle phenomenon existing in dimensions higher than 2D (the case of light speckle).

3D plane-based egomotion for SLAM on semi-structured environment

Several works deal with 3D data in SLAM problem. Data come from a 3D laser sweeping unit or a stereo camera, both providing a huge amount of data. In this paper, we detail an efficient method to extract planar patches from 3D raw data. Then, we use these patches in an ICP-like method in order to address the SLAM problem. Using ICP with planes is not a trivial task. It needs some adaptation from the original ICP. Some promising results are shown for outdoor environment.

Pitfalls on computing liquid-liquid phase equilibria using the k-value method

Poster presented in the 11th Mediterranean Congress of Chemical Engineering, Barcelona, October 21-24, 2008.

Efficient sampling of band-limited signals from sine wave crossings

This correspondence presents an efficient method for reconstructing a band-limited signal in the discrete domain from its crossings with a sine wave. The method makes it possible to design A/D converters that only deliver the crossing timings, which are then used to interpolate the input signal at arbitrary instants. Potentially, it may allow for reductions in power consumption and complexity in these converters. The reconstruction in the discrete domain is based on a recently-proposed modification of the Lagrange interpolator, which is readily implementable with linear complexity and efficiently, given that it re-uses known schemes for variable fractional-delay (VFD) filters. As a spin-off, the method allows one to perform spectral analysis from sine wave crossings with the complexity of the FFT. Finally, the results in the correspondence are validated in several numerical examples.

Shear wave velocity estimation in the metropolitan area of Málaga (S Spain)

We carry out a seismic noise study based on array measurements at three sites in the Málaga basin, South Spain, for the further estimation of shear wave velocity profiles. For this purpose, we use both the H/V method and the f–k technique in order to characterize the different materials present in the zone, i.e., Quaternary sediments and Pliocene sedimentary rocks above the bedrock. The H/V analysis shows frequency peaks going from 1 Hz, in areas close to the border of the basin, to 0.3 Hz in places located toward the center of the formation. The f–k analysis allows obtaining the dispersion curves associated with each site and subsequently, estimating the Vs profiles by inversion of the respective group velocities. In this way, the basin basement can be characterized by S-wave velocities greater than 2000 m/s. Regarding the basin fill, it is divided into three layers defined by different wave velocity intervals. The shallowest one is featured by velocities ranging from 150 to 400 m/s and comprises the Quaternary sediments, while velocities going from 550–700 to1200–1600 m/s characterize the two underlying layers composed by Pliocene sediments. Finally, the information provided by the three Vs profiles is integrated in a 2D cross-section of the basin to have a spatial view of its sedimentary structure. The results obtained here, in addition to providing useful information about the infill of the basin near the metropolitan area of Málaga, will be very helpful for future seismic zonation studies in the region.

Light capsules shaped by curvilinear meta-surfaces

We propose a simple yet efficient method for generating in-plane hollow beams with a nearly full circular light shell without the contribution of backward propagating waves. The method relies on modulating the phase in the near field of a centrosymmetric optical wave front, such as that from a high-numerical-aperture focused wave field. We illustrate how beam acceleration may be carried out by using an ultranarrow non-flat meta-surface formed by engineered plasmonic nanoslits. A mirror-symmetric, with respect to the optical axis, circular caustic surface is numerically demonstrated that can be used as an optical bottle.