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.

Análisis de la influencia de los diferentes parámetros de diseño en la aparición de fatiga superficial en contactos mecánicos

La fatiga superficial es uno de los principales problemas en las transmisiones mecánicas y es uno de los focos de atención de las investigaciones de los últimos anos en Tribología. La disminución de viscosidad de los lubricantes para la mejora de la eficiencia, el aumento de las potencias a transmitir, el aumento de la vida de los componentes o la mejora de su fiabilidad han supuesto que los fenómenos de fatiga superficial hayan cobrado especial relevancia, especialmente los fenómenos de pitting y micropitting en cajas multiplicadoras/reductoras de grandes potencias de aplicación, por ejemplo, en el sector eólico. Como todo fenómeno de fatiga, el pitting y micropitting son debidos a la aplicación de cargas ciclicas. Su aparición depende de las presiones y tensiones cortantes en el contacto entre dos superficies que al encontrarse en rodadura y deslizamiento varian con el tiempo. La principal consecuencia de la fatiga superficial es la aparición de hoyuelos de diferente magnitud segun la escala del fenómeno (pitting o micropitting) en la superficie del material. La aparición de estos hoyuelos provoca la perdida de material, induce vibraciones y sobrecargas en el elemento que finalmente acaba fallando. Debido a la influencia de la presión y tensión cortante en el contacto, la aparición de fatiga depende fuertemente del lubricante que se encuentre entre las dos superficies y de las condiciones de funcionamiento en las cuales este trabajando. Cuando el contacto trabaja en condiciones de lubricacion mixta-elastohidrodinamica tiende a aparecer micropitting debido a las altas tensiones localizadas en las proximidades de las asperezas, mientras que si el régimen es de lubricación completa el tipo de fatiga superficial suele ser pitting debido a las tensiones mas suavizadas y menos concentradas. En esta Tesis Doctoral se han analizado todos estos factores de influencia que controlan el pitting y el micropitting prestando especial atención al efecto del lubricante. Para ello, se ha dado un enfoque conjunto a ambos fenómenos resolviendo las ecuaciones involucradas en el contacto elastohidrodinamico no-Newtoniano (la ecuación de Reynolds, la deformación elástica de los sólidos y la reologia del lubricante) para conocer la presión y la tensión cortante en el contacto. Conocidas estas, se resuelve el campo de tensiones en el interior del material y, finalmente, se aplican criterios de fatiga multiaxial (Crossland, Dang Van y Liu-Mahadevan) para conocer si el material falla o no falla. Con la metodología desarrollada se ha analizado el efecto sobre las tensiones y la aparición de la fatiga superficial del coeficiente viscosidad-presion, de la compresibilidad, del espesor especifico de película y de la fricción así como de la influencia de las propiedades a fatiga del material y de las condiciones de funcionamiento (radios de contacto, velocidad, deslizamiento, carga y temperatura). Para la validación de los resultados se han utilizado resultados teóricos y experimentales de otros autores junto con normas internacionales de amplia utilización en el mundo industrial, entre otras, para el diseño y calculo de engranajes. A parte del trabajo realizado por simulación y cálculo de los diferentes modelos desarrollados, se ha realizado un importante trabajo experimental que ha servido no solo para validar la herramienta desarrollada sino que además ha permitido incorporar al estudio factores no considerados en los modelos, como los aditivos del lubricante. Se han realizado ensayos de medida del coeficiente de fricción en una maquina de ensayo puntual con la que se ha validado el cálculo del coeficiente de fricción y se ha desarrollado un proceso de mejora del coeficiente de fricción mediante texturizado superficial en contactos puntuales elastohidrodinamicos mediante fotolitografia y ataque quimico. Junto con los ensayos de medida de fricción en contacto puntual se han realizado ensayos de fricción y fatiga superficial en contacto lineal mediante una maquina de discos que ha permitido evaluar la influencia de diferentes aditivos (modificadores de fricción, antidesgaste y extrema-presion) en la aparición de fatiga superficial (pitting y micropitting) y la fricción en el contacto. Abstract Surface fatigue is one of the most important problems of mechanical transmissions and therefore has been one of the main research topics on Tribology during the last years. On the one hand, industrial demand on fuel economy has led to reduce lubricant viscosity in order to improve efficiency. On the other hand, the requirements of power and life of machine elements are continuously increasing, together with the improvements in reliability. As a consequence, surface fatigue phenomena have become critical in machinery, in particular pitting and micropitting in high power gearboxes of every kind of machines, e.g., wind turbines or cranes. In line with every fatigue phenomena, pitting and micropitting are caused by cyclic loads. Their appearance depends on the evolution of pressures and shear stresses with time, throughout the contact between surfaces under rolling and sliding conditions. The main consequence of surface fatigue is the appearance of pits on the surface. The size of the pits is related to the scale of the fatigue: pitting or micropitting. These pits cause material loss, vibrations and overloads until the final failure is reached. Due to the great influence of the pressures and shear stresses in surface fatigue, the appearance of pits depends directly on the lubricant and the operating conditions. When the contact works under mixed regime (or under elastohydrodynamic but close to mixed regime) the main fatigue failure is micropitting because of the high pressures located near the asperities. In contrast, when the contact works under elastohydrodynamic fully flooded conditions the typical fatigue failure is pitting. In this Ph.D. Thesis, the main factors with influence on pitting and micropitting phenomena are analyzed, with special attention to the effect of the lubricant. For this purpose, pitting and micropitting are studied together by solving the equations involved in the non-Newtonian elastohydrodynamic contact. Thus, pressure and shear stress distributions are found by taking into account Reynolds equation, elastic deflection of the solids and lubricant rheology. Subsequently, the stress field inside the material can be calculated and different multiaxial fatigue criteria (Crossland, Dang Van and Liu- Mahadevan) can be applied to predict whether fatigue failure is reached. The influences of the main parameters on pressure and surface fatigue have been studied, taking into account the lubricant compressibility and its viscosity-pressure coefficient, the specific film thickness, the friction coefficient and the fatigue properties of the contacting materials, together with the operating conditions (contact radius, mean velocity, sliding velocity, load and temperature). Several theoretical and experimental studies of different authors have been used to validate all the results obtained, together with international standards used worldwide in gear design industry. Moreover, an experimental stage has been carried out in order to validate the calculation methods and introduce additional influences not included previously, e.g., lubricant additives. The experimentation includes different friction tests in point contacts performed with a tribological equipment in order to validate the results given by the calculations. Furthermore, the reduction and optimization of the friction coefficient is analyzed by means of textured surfaces, obtained combining photolithography and chemical etching techniques. Besides the friction tests with point contact, friction and surface fatigue tests have also been performed with line contact in a tribological test rig. This equipment is also used to study the influence of different types of additives (friction modifiers, anti-wear and extreme-pressure additives) on surface fatigue (pitting and micropitting).

Physics and applications of electrodynamic space tethers

Basic effects and dynamical and electrical contact issues in the physics of (electrodynamic space) bare tethers are discussed. Scientific experiments and powerpropulsion applications, including a paradoxical use of bare tethers in outer-planet exploration,are considered.

On super free fall.

Villermaux & Pomeau (J. Fluid Mech., vol. 642, 2010, p. 147) analysed the motion of the interface of an inviscid liquid column released from rest in a vertical tube whose area expands gradually downwards, with application to an inverted conical container for which experimental measurements were carried out. An error in the analysis is found and corrected in the present investigation, which provides the new governing equation for the super-accelerated interface motion down gradually varying tubes in general, and integrated results for interface trajectories, velocities and accelerations down a conical tube in particular. Interestingly, the error does not affect any of the conclusions given in the 2010 paper. Further new results are reported here such as the equation governing the centre of mass and proof that the end point acceleration is exactly that of gravity

Laminar-turbulent transition induced by a discrete roughness element in a supersonic boundary layer

The linear instability and breakdown to turbulence induced by an isolated roughness element in a boundary layer at Mach 2:5, over an isothermal flat plate with laminar adiabatic wall temperature, have been analysed by means of direct numerical simulations, aided by spatial BiGlobal and three-dimensional parabolized (PSE-3D) stability analyses. It is important to understand transition in this flow regime since the process can be slower than in incompressible flow and is crucial to prediction of local heat loads on next-generation flight vehicles. The results show that the roughness element, with a height of the order of the boundary layer displacement thickness, generates a highly unstable wake, which is composed of a low-velocity streak surrounded by a three-dimensional high-shear layer and is able to sustain the rapid growth of a number of instability modes. The most unstable of these modes are associated with varicose or sinuous deformations of the low-velocity streak; they are a consequence of the instability developing in the three-dimensional shear layer as a whole (the varicose mode) or in the lateral shear layers (the sinuous mode). The most unstable wake mode is of the varicose type and grows on average 17% faster tan the most unstable sinuous mode and 30 times faster than the most unstable boundary layer mode occurring in the absence of a roughness element. Due to the high growthrates registered in the presence of the roughness element, an amplification factor of N D 9 is reached within 50 roughness heights from the roughness trailing edge. The independently performed Navier–Stokes, spatial BiGlobal and PSE-3D stability results are in excellent agreement with each other, validating the use of simplified theories for roughness-induced transition involving wake instabilities. Following the linear stages of the laminar–turbulent transition process, the roll-up of the three-dimensional shear layer leads to the formation of a wedge of turbulence, which spreads laterally at a rate similar to that observed in the case of compressible turbulent spots for the same Mach number.

On linear instability mechanisms in incompressible open cavity flow

A theoretical study of linear global instability of incompressible flow over a rectangular spanwise-periodic open cavity in an unconfined domain is presented. Comparisons with the limited number of results available in the literature are shown. Subsequently, the parameter space is scanned in a systematic manner, varying Reynolds number, incoming boundary-layer thickness and length-to-depth aspect ratio. This permits documenting the neutral curves and leading eigenmode characteristics of this flow. Correlations constructed using the results obtained collapse all available theoretical data on the three-dimensional instabilities.

The breaking of axisymmetric slender liquid bridges

Liquids held by surface tension forces can bridge the gap between two solid bodies placed not too far apart from each other. The equilibrium conditions and stability criteria for static, cylindrical liquid bridges are well known. However, the behaviour of an unstable liquid bridge, regarding both its transition toward breaking and the resulting configuration, is a matter for discussion. The dynamical problem of axisymmetric rupture of a long liquid bridge anchored at two equal coaxial disks is treated in this paper through the adoption of one-dimensional theories which are widely used in capillary jet problems

One-dimensional self-similar solution of the dynamics of axisymmetric slender liquid bridges

In this paper the dynamics of axisymmetric, slender, viscous liquid bridges having volume close to the cylindrical one, and subjected to a small gravitational field parallel to the axis of the liquid bridge, is considered within the context of one-dimensional theories. Although the dynamics of liquid bridges has been treated through a numerical analysis in the inviscid case, numerical methods become inappropriate to study configurations close to the static stability limit because the evolution time, and thence the computing time, increases excessively. To avoid this difficulty, the problem of the evolution of these liquid bridges has been attacked through a nonlinear analysis based on the singular perturbation method and, whenever possible, the results obtained are compared with the numerical ones.

Numerical and experimental study of the dynamics of axisymmetric slender liquid bridges

A one-dimensional inviscid slice model has been used to study numerically the influence of axial microgravity on the breaking of liquid bridges having a volume close to that of gravitationless minimum volume stability limit. Equilibrium shapes and stability limits have been obtained as well as the dependence of the volume of the two drops formed after breaking on both the length and the volume of the liquid bridge. The breaking process has also been studied experimentally. Good agreement has been found between theory and experiment for neutrally buoyant systems

One-dimensional linear analysis of the compound jet

The stability of an infinitely long compound liquid column is analysed by using a one-dimensional inviscid slice model. Results obtained from this one-dimensional linear analysis are applicable to the study of compound capillary jets, which are used in the ink-jet printing technique. Stability limits and the breaking regimes of such fluid configurations are established, and, whenever possible, theoretical results are compared with experimental ones.

Quasi-cylindrical approximation to the swirling flow in an atomizer chamber

A quasi-cylindrical approximation is used to analyse the axisymmetric swirling flow of a liquid with a hollow air core in the chamber of a pressure swirl atomizer. The liquid is injected into the chamber with an azimuthal velocity component through a number of slots at the periphery of one end of the chamber, and flows out as an anular sheet through a central orifice at the other end, following a conical convergence of the chamber wall. An effective inlet condition is used to model the effects of the slots and the boundary layer that develops at the nearby endwall of the chamber. An analysis is presented of the structure of the liquid sheet at the end of the exit orifice, where the flow becomes critical in the sense that upstream propagation of long-wave perturbations ceases to be possible. This nalysis leads to a boundary condition at the end of the orifice that is an extension of the condition of maximum flux used with irrotational models of the flow. As is well known, the radial pressure gradient induced by the swirling flow in the bulk of the chamber causes the overpressure that drives the liquid towards the exit orifice, and also leads to Ekman pumping in the boundary layers of reduced azimuthal velocity at the convergent wall of the chamber and at the wall opposite to the exit orifice. The numerical results confirm the important role played by the boundary layers. They make the thickness of the liquid sheet at the end of the orifice larger than predicted by rrotational models, and at the same time tend to decrease the overpressure required to pass a given flow rate through the chamber, because the large axial velocity in the boundary layers takes care of part of the flow rate. The thickness of the boundary layers increases when the atomizer constant (the inverse of a swirl number, proportional to the flow rate scaled with the radius of the exit orifice and the circulation around the air core) decreases. A minimum value of this parameter is found below which the layer of reduced azimuthal velocity around the air core prevents the pressure from increasing and steadily driving the flow through the exit orifice. The effects of other parameters not accounted for by irrotational models are also analysed in terms of their influence on the boundary layers.

The influence of the outer bath in the dynamics of axisymmetric liquid bridges

The main effects on the dynamics of a liquid bridge due to the presence of an outer liquid, as occur in experiments using the Plateau-tank technique, are considered. The one-dimensional nonlinear model developed here allows us to perform the computation of both breaking processes and oscillatory motions of slender liquid bridges, although in this paper only the results concerning breaking processes are reported. Additionally,the oscillatory motions are studied both experimentally and by using a new linear model. Results from both sources show good agreement

On the controllability of distributed systems

To “control” a system is to make it behave (hopefully) according to our “wishes,” in a way compatible with safety and ethics, at the least possible cost. The systems considered here are distributed—i.e., governed (modeled) by partial differential equations (PDEs) of evolution. Our “wish” is to drive the system in a given time, by an adequate choice of the controls, from a given initial state to a final given state, which is the target. If this can be achieved (respectively, if we can reach any “neighborhood” of the target) the system, with the controls at our disposal, is exactly (respectively, approximately) controllable. A very general (and fuzzy) idea is that the more a system is “unstable” (chaotic, turbulent) the “simplest,” or the “cheapest,” it is to achieve exact or approximate controllability. When the PDEs are the Navier–Stokes equations, it leads to conjectures, which are presented and explained. Recent results, reported in this expository paper, essentially prove the conjectures in two space dimensions. In three space dimensions, a large number of new questions arise, some new results support (without proving) the conjectures, such as generic controllability and cases of decrease of cost of control when the instability increases. Short comments are made on models arising in climatology, thermoelasticity, non-Newtonian fluids, and molecular chemistry. The Introduction of the paper and the first part of all sections are not technical. Many open questions are mentioned in the text.

Obtenção de gel mucoadesivo de nistatina para o tratamento da candidíase oral. Desenvolvimento e caracterização de dispersões sólidas de nistatina

A nistatina (NYS) é o fármaco de primeira escolha no tratamento da candidíase oral, que frequentemente acomete mais os indivíduos imunocomprometidos e pacientes com outras desordens (diabetes não tratada, neoplasias, imunodeficiências). No mercado brasileiro, a NYS é encontrada na forma de suspensão oral aquosa, onde o procedimento para sua administração consiste em bochechar o medicamento. Apesar de haver a indicação de que se mantenha o contato direto entre fármaco e a mucosa oral, na qual se encontra a Candida spp., o que aumentaria expressivamente o sucesso terapêutico, a suspensão não apresenta tal propriedade. Assim, a NYS que é fármaco com ação efetiva contra a candidíase oral, é considerada pertencente à Classe IV do Sistema de Classificação Biofarmacêutica, ou seja, apresenta baixa solubilidade e baixa permeabilidade. A baixa solubilidade pode comprometer sua disponibilidade na cavidade oral, e consequentemente, sua ação farmacológica. Diante desse quadro, o objetivo do presente trabalho foi o desenvolvimento de dispersões sólidas de NYS para o tratamento da candidíase oral, e sua posterior incorporação em gel mucoadesivo oral, favorecendo a formulação no local de ação. As dispersões sólidas são sistemas farmacêuticos, onde um fármaco pouco solúvel em água encontra-se dispersado em um carreador, no estado sólido. Os carreadores normalmente são hidrofílicos, o que permite que esses sistemas sejam empregados para aumentar a solubilidade aquosa do fármaco. Assim, foram desenvolvidas as dispersões sólidas de NYS, pelo método de eliminação do solvente, empregando como carreadores, lactose, HPMC, poloxamer 407 e poloxamer 188. Essas foram submetidas à caracterização por análise térmica, usando os ensaios de calorimetria exploratória diferencial (DSC) e termogravimetria/termogravimetria derivada (TG/DTG). Dentre essas dispersões sólidas, aquelas que se mostraram com comportamento térmico sugerindo a formação de um novo \"sistema\", foram analisadas por meio de ensaio de solubilidade. Dessa forma, a formulação NYS DS G2 (49) se destacou, pois apresentou maior solubilidade em água (4,484 mg/mL); em pH 5,5 (4,249 mg/mL) e em pH 7,0 (4,293 mg/mL), ou seja, houve um aumento de 1,426 vezes em água; 4,227 vezes em pH 5,5; e 2,743 vezes em pH 7,0. Essa formulação foi, por fim avaliada por difração de raio-X e espectroscopia de infravermelho com transformada de Fourier, técnicas que corroboraram com a análise térmica quanto à indicação de formação da dispersão sólida. Por sua vez, essa dispersão sólida foi incorporada em 4 bases de géis mucoadesivos de carbopol ® 934 PNF, alterando apenas a concentração do polímero (0,5; 1,0; 1,5; 2,0 %p/p). Foi observado que a liberação de NYS DS G2 (49) foi superior, quando comparada à liberação de NYS MP a partir do gel, e através do ensaio de mucoadesão, percebeu-se que os géis desenvolvidos apresentaram propriedades mucoadesivas compatíveis com relatos na literatura, independentemente da quantidade de carbopol ® empregada. As características reológicas foram distintas, e foi observado que as formulações Gel A e Gel B, que possuem menor quantidade de polímero, tiverem um indicativo de comportamento de fluido newtoniano, diferente dos demais, o que pode não ser desejado para esse tipo de forma farmacêutica tópica e semi-sólida. Ao final desse trabalho, pode-se concluir que foi possível desenvolver um sistema farmacêutico na forma de dispersão sólida com maior solubilidade que a NYS pura, e sua incorporação em uma forma farmacêutica mucoadesiva, e que a liberação da NYS na forma DS foi muito superior que o fármaco na forma \"convencional\", o que permite que a NYS esteja mais disponível na cavidade oral, e também junto à mucosa bucal, o que levaria a efeito farmacológico mais efetivo do antifúngico.

Transporte hidráulico de misturas bentonita-água em condutos forçados

A argila bentonítica é amplamente utilizada em transporte de sólidos produzidos durante perfuração de poços. Teve por objetivo estudar o escoamento de misturas bentonita-água e determinar suas propriedades reológicas e parâmetros hidráulicos úteis nos projetos de instalações de recalque de misturas sólido-líquido. Foi montado um circuito fechado de tubulações para estudar dados de perda de carga e perfis de velocidade. Realizaram-se ensaios com misturas bentonita-água sob varias concentrações, algumas transportando areia. Observaram-se que a reologia da mistura bentonita-água é melhor descrita pela formulação de Herschell-Bulkley para fluidos não-Newtonianos. O coeficiente de atrito para descrever a perda de carga da mistura bentonita-água observada em tubulações no laboratório coloca-se entre as previsões de Tomita (1959) e Szilas et al (1981). A variação da velocidade da mistura na seção transversal do tubo é melhor aproximada pela equação de Bogue-Metzner (1963).