A two-fluid model for refrigerant flow through adiabatic capillary tubes

This work presents a numerical model to simulate refrigerant flow through capillary tubes, commonly used as expansion devices in refrigeration systems. The flow is divided in a single-phase region, where the refrigerant is in the subcooled liquid state, and a region of two-phase flow. The capillary tube is considered straight and horizontal. The flow is taken as one-dimensional and adiabatic. Steady-state condition is also assumed and the metastable flow phenomena are neglected. The two-fluid model, considering the hydrodynamic and thermal non-equilibrium between the liquid and vapor phases, is applied to the two-phase flow region. Comparisons are made with experimental measurements of the mass flow rate and pressure distribution along two capillary tubes working with refrigerant R-134a in different operating conditions. The results indicate that the present model provides a better estimation than the commonly employed homogeneous model. Some computational results referring to the quality, void fraction, velocities, and temperatures of each phase are presented and discussed.

Postura sentada de crianças com paralisia cerebral: influência da utilização do apoio de membros superiores

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

Análise do equilíbrio estático e de deformidades nos pés de crianças com paralisia cerebral por meio de baropodometria computadorizada

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

The effects of the pressing step on the microstructure and aging of NdFeB bonded magnets

The effects of the compaction step on the (micro)structural features and aging behavior of polymer coated NdFeB-based bonded magnets is reported. Due to the fracture of the material during pressing, it is estimated an increase of at least 14% in the particles' area which is not coated. Such uncoated surfaces, when exposed to the environment, reduce the magnetic performance of the magnets aged/cured in air by 19% in the conditions evaluated in this investigation. Furthermore, XRD results interpreted by Rietveld analyses show a lattice parameter change in the tetragonal structure of the hard magnetic phase after pressing. Such change varies as a function of the height of the compacted part and it is ascribed to macro-elastic stress arising from the pressure distribution in the magnet. An aging/curing step during 24 h is able to relief such macro-elastic stress. (C) 2012 Elsevier B.V. All rights reserved.

Hydroelastic impacts of deformable wedges

This work investigates the slamming phenomenon experienced during the water entry of deformable bodies. Wedges are chosen as reference geometry due to their similarity to a generic hull section. Hull slamming is a phenomenon occurring when a ship re-enters the water after having been partially or completely lifted out the water. While the analysis of rigid structures entering the water has been extensively studied in the past and there are analytical solutions capable of correctly predicting the hydrodynamic pressure distribution and the overall impact dynamics, the effect of the structural deformation on the structural force is still a challenging problem to be solved. In fact, in case of water impact of deformable bodies, the dynamic deflection could interact with the fluid flow, changing the hydrodynamic load. This work investigates the hull-slamming problem by experiments and numerical simulations of the water entry of elastic wedges impacting on an initially calm surface. The effect of asymmetry due to horizontal velocity component or initial tilt angle on the impact dynamics is also studied. The objective of this work is to determine an accurate model to predict the overall dynamics of the wedge and its deformations. More than 1200 experiments were conducted by varying wedge structural stiffness, deadrise angle, impact velocity and mass. On interest are the overall impact dynamics and the local structural deformation of the panels composing the wedge. Alongside with the experimental analysis, numerical simulations based on a coupled Smoothed Particle Hydrodynamics (SPH) and FEM method are developed. The experimental results provide evidence of the mutual interaction between hydrodynamic load and structural deformation. It is found a simple criterion for the onset of fluid structure interaction (FSI), giving reliable information on the cases where FSI should been taken into account.

Finite element based estimation of contact areas and pressures of the human scaphoid in various functional positions of the hand

The scaphoid is the most frequently fractured carpal bone. When investigating fixation stability, which may influence healing, knowledge of forces and moments acting on the scaphoid is essential. The aim of this study was to evaluate cartilage contact forces acting on the intact scaphoid in various functional wrist positions using finite element modeling. A novel methodology was utilized as an attempt to overcome some limitations of earlier studies, namely, relatively coarse imaging resolution to assess geometry, assumption of idealized cartilage thicknesses and neglected cartilage pre-stresses in the unloaded joint. Carpal bone positions and articular cartilage geometry were obtained independently by means of high resolution CT imaging and incorporated into finite element (FE) models of the human wrist in eight functional positions. Displacement driven FE analyses were used to resolve inter-penetration of cartilage layers, and provided contact areas, forces and pressure distribution for the scaphoid bone. The results were in the range reported by previous studies. Novel findings of this study were: (i) cartilage thickness was found to be heterogeneous for each bone and vary considerably between carpal bones; (ii) this heterogeneity largely influenced the FE results and (iii) the forces acting on the scaphoid in the unloaded wrist were found to be significant. As major limitations, accuracy of the method was found to be relatively low, and the results could not be compared to independent experiments. The obtained results will be used in a following study to evaluate existing and recently developed screws used to fix scaphoid fractures.

A longitudinal study of Type-B aortic dissection and endovascular repair scenarios: computational analyses

Conservative medical treatment is commonly first recommended for patients with uncomplicated Type-B aortic dissection (AD). However, if dissection-related complications occur, endovascular repair or open surgery is performed. Here we establish computational models of AD based on radiological three-dimensional images of a patient at initial presentation and after 4-years of best medical treatment (BMT). Computational fluid dynamics analyses are performed to quantitatively investigate the hemodynamic features of AD. Entry and re-entries (functioning as entries and outlets) are identified in the initial and follow-up models, and obvious variations of the inter-luminal flow exchange are revealed. Computational studies indicate that the reduction of blood pressure in BMT patients lowers pressure and wall shear stress in the thoracic aorta in general, and flattens the pressure distribution on the outer wall of the dissection, potentially reducing the progressive enlargement of the false lumen. Finally, scenario studies of endovascular aortic repair are conducted. The results indicate that, for patients with multiple tears, stent-grafts occluding all re-entries would be required to effectively reduce inter-luminal blood communication and thus induce thrombosis in the false lumen. This implicates that computational flow analyses may identify entries and relevant re-entries between true and false lumen and potentially assist in stent-graft planning.

Vessel orientation-dependent sensitivity of optoacoustic imaging using a linear array transducer

For clinical optoacoustic imaging, linear probes are preferably used because they allow versatile imaging of the human body with real-time display and free-hand probe guidance. The two-dimensional (2-D) optoacoustic image obtained with this type of probe is generally interpreted as a 2-D cross-section of the tissue just as is common in echo ultrasound. We demonstrate in three-dimensional simulations, phantom experiments, and in vivo mouse experiments that for vascular imaging this interpretation is often inaccurate. The cylindrical blood vessels emit anisotropic acoustic transients, which can be sensitively detected only if the direction of acoustic radiation coincides with the probe aperture. Our results reveal for this reason that the signal amplitude of different blood vessels may differ even if the vessels have the same diameter and initial pressure distribution but different orientation relative to the imaging plane. This has important implications for the image interpretation, for the probe guidance technique, and especially in cases when a quantitative reconstruction of the optical tissue properties is required.

Two-dimensional and three-dimensional oxygen mapping by 3He-MRI validation in a lung phantom

The aim of this study was to validate oxygen-sensitive 3He-MRI in noninvasive determination of the regional, two- and three-dimensional distribution of oxygen partial pressure. In a gas-filled elastic silicon ventilation bag used as a lung phantom, oxygen sensitive two- and three-dimensional 3He-MRI measurements were performed at different oxygen concentrations which had been equilibrated in a range of normal and pathologic values. The oxygen partial pressure distribution was determined from 3He-MRI using newly developed software allowing for mapping of oxygen partial pressure. The reference bulk oxygen partial pressure inside the phantom was measured by conventional respiratory gas analysis. In two-dimensional measurements, image-based and gas-analysis results correlated with r=0.98; in three-dimensional measurements the between-methods correlation coefficient was r=0.89. The signal-to-noise ratio of three-dimensional measurements was about half of that of two-dimensional measurements and became critical (below 3) in some data sets. Oxygen-sensitive 3He-MRI allows for noninvasive determination of the two- and three-dimensional distribution of oxygen partial pressure in gas-filled airspaces.

Kinetics and kinematics of the tölt: effects of rider interaction and shoeing manipulations

Introduction . Compared to most equine horse breeds which are able to walk, trot and canter /gallop, the gait repertoire of the Icelandic horses additionally includes the lateral gait tölt and frequently also the pace. With respect to the tölt gait, special shoeing, saddling and riding techniques have been developed for Icelandic horses in order to enhance its expressiveness and regularity. Toes are left unnaturally long and heavy shoes and paddings, as well as weighted boots are used to enforce the individual gait predisposition. For the same reason, the rider is placed more caudally to the horse's centre of mass as compared to other riding techniques. The biomechanical impact of these methods on the health of the locomotor system has so far never been subject of systematic research. Objectives . The aims of the presented study are (1) to describe the kinetic and kinematic characteristics of the tölt performed on a treadmill, (2) to understand the mechanical consequences of shoeing manipulation (long hooves, weighted boots) on the loading and protraction movement of the limbs, as well as (3) to study the pressure distribution and effects on the gait pattern of 3 different saddle types used for riding Icelandic horses. Materials and methods . Gait analysis was carried out in 13 Icelandic horses at walk and at slow and medium tölting and trotting speeds on a high-speed treadmill instrumented for measuring vertical ground reaction forces as well as temporal and spatial gait variables. Kinematic data of horse, rider and saddle were measured simultaneously. Gait analysis was first carried out with high, long hooves (SH) without and in combination with weighted boots (ad aim (2)). Afterwards, horses were re-shod according to current horseshoeing standards (SN) and gait analysis was repeated (ad aims (1) and (2)). In a second trial, horses were additionally equipped with a pressure sensitive saddle mat and were ridden with a dressage-like saddle (SDres), an Icelandic saddle (Slcel) and a saddle cushion (SCush) in the standard saddle position (ad aim 3). Results and conclusions . Compared to trot at the same speed, tölting horses had a higher stride rate and lower stride impulses. At the tölt loading of the forelimbs was increased in form of higher peak vertical forces (Fzpeak) due to shorter relative stance durations (StDrel). Conversely, in the hindlimbs, longer StDrel resulted in lower Fzpeak. Despite the higher head-neck position at tölt, there was no measurable shift in weight to the hindlimbs. Footfall rhythm was in most horses laterally coupled at the tölt and frequently had a slight fourbeat and a very short suspension phase at trot; underlining the fact that performance of correct gaits in Icelandic horses needs special training. Gait performance as it is currently judged in competition could be improved using a shoeing with SH, resulting in a 21 ± 5 mm longer dorsal hoof wall, but also a weight gain of 273 ± 50 g at the distal limb due to heavier shoeing material. Compared to SN, SH led to a lower stride rate, a longer stride length and a higher, but not wider, forelimb protraction arc, which were also positively associated with speed. At the tölt, the footfall rhythm showed less tendency to lateral couplets and at the trot, the suspension phase was longer. However, on the long term, SH may have negative implications for the health of the palmar structures of the distal foot by increased limb impulses, higher torques at breakover (up to 20%); as well as peak vertical forces at faster speeds. Compared to the shoeing style, the saddle type had less influence on limb forces or movements. The slight weight shift to the rear with SCush and Slcel may be explained by the more caudal position of the rider relative to the horse's back. With SCush, pressure was highest under the cranial part of the saddle, whereas the saddles with trees had more pressure under the caudal area.

Numerical and Experimental Studies of Sail Aerodynamics

The purpose of this investigation was the determination of the aerodynamic performance of sails and gain knowledge of the phenomena involved in order to improve the aerody¬namic characteristics. In this research, the airflow around different sails in four scenarios was studied. The method to analyze these scenarios was the combination of numerical simulations and experimental tests by taking advantage of the best of each tool. Two different Com¬putational Fluid Dynamic codes were utilized: the ANSYS-CFX and the CD-Adapco’s STAR-CCM+. The experimental tests were conducted in the Atmospheric Boundary Layer Wind Tunnel at the Universidad de Granada (Spain), the Twisted Flow Wind Tunnel at the University of Auckland (New Zealand) and the A9 Wind Tunnel at the Universidad Polit´ecnica de Madrid (Spain). Through this research, it was found the three-dimensional effect of the mast on the aerodynamic performance of an IMS Class boat. The pressure distribution on a Transpac 52 Class mainsail was also determined. Moreover, the aerodynamic perfor¬mance of the 43ft and 60ft Dhow Classes was obtained. Finally, a feasibility study was conducted to use a structural wing in combination with conventional propulsions systems. The main conclusion was that this research clarified gaps on the knowledge of the aerodynamic performance of sails. Moreover, since commercial codes were not specifically designed to study sails, a procedure was developed. On the other hand, innovative experimental techniques were used and applied to model-scale sails. The achievements of this thesis are promising and some of the results are already in use by the industry on a daily basis. El propósito de este estudio era determinar el comportamiento aerodinámico de unas velas y mejorar el conocimiento de los fenómenos que suceden para optimizar las características aerodinámicas de dichas velas. En esta investigación se estudió el flujo de aire alrededor de diferentes velas en cuatro escenarios. El método para analizar estos escenarios fue la combinación de simulaciones numéricas y ensayos experimentales mediante el aprovechamiento de las ventajas de cada herramienta. Se utilizaron dos códigos de dinámica de fluidos computacional: el ANSYS-CFX y el STAR-CCM+ de la empresa CD-Adapco. Los ensayos experimentales se desarrollaron en el túnel de viento de capa límite de la Universidad de Granada (España), el túnel de viento de la Universidad de Auckland (Nueva Zelanda) y en el túnel A9 de la Universidad Politécnica de Madrid (España). Mediante esta investigación, se determinó el efecto tridimensional del mástil en un velero de la clase IMS. También se describió la distribución de presiones sobre una mayor de un Transpac 52. Además, se obtuvo el comportamiento aerodinámico de las clases 43ft y 60ft de los veleros Dhows. Finalmente, se llevó a cabo un estudio de viabilidad de la utilización de un ala estructural en combinación con sistemas de propulsión convencionales. La conclusión principal de esta investigación fue la capacidad de explicar ciertas lagunas en el conocimiento del comportamiento aerodinámico de las velas en diferentes escenarios. Además, dado que los códigos comerciales no están específicamente diseñados para el estudio de velas, se desarrolló un procedimiento a tal efecto. Por otro lado, se han utilizado innovadoras técnicas experimentales y se han aplicado a modelos de velas a escala. Los logros de esta investigación son prometedores y algunos de los resultados obtenidos ya están siendo utilizados por la industria en su día a día.

Consolidation problems

The analysis of deformation in soils is of paramount importance in geotechnical engineering. For a long time the complex behaviour of natural deposits defied the ingenuity of engineers. The time has come that, with the aid of computers, numerical methods will allow the solution of every problem if the material law can be specified with a certain accuracy. Boundary Techniques (B.E.) have recently exploded in a splendid flowering of methods and applications that compare advantegeously with other well-established procedures like the finite element method (F.E.). Its application to soil mechanics problems (Brebbia 1981) has started and will grow in the future. This paper tries to present a simple formulation to a classical problem. In fact, there is already a large amount of application of B.E. to diffusion problems (Rizzo et al, Shaw, Chang et al, Combescure et al, Wrobel et al, Roures et al, Onishi et al) and very recently the first specific application to consolidation problems has been published by Bnishi et al. Here we develop an alternative formulation to that presented in the last reference. Fundamentally the idea is to introduce a finite difference discretization in the time domain in order to use the fundamental solution of a Helmholtz type equation governing the neutral pressure distribution. Although this procedure seems to have been unappreciated in the previous technical literature it is nevertheless effective and straightforward to implement. Indeed for the special problem in study it is perfectly suited, because a step by step interaction between the elastic and flow problems is needed. It allows also the introduction of non-linear elastic properties and time dependent conditions very easily as will be shown and compares well with performances of other approaches.

Mathematical Analysis of the Effect of Rotor Geometry on Cup Anemometer Response

The calibration coefficients of two commercial anemometers equipped with different rotors were studied. The rotor cups had the same conical shape, while the size and distance to the rotation axis varied.The analysis was based on the 2-cup positions analytical model, derived using perturbation methods to include second-order effects such as pressure distribution along the rotating cups and friction.Thecomparison with the experimental data indicates a nonuniformdistribution of aerodynamic forces on the rotating cups, with higher forces closer to the rotating axis. The 2-cup analytical model is proven to be accurate enough to study the effect of complex forces on cup anemometer performance.

A parametric, experimental analysis of conical vortices on curved roofs of low-rise buildings

Different methods to reduce the high suction caused by conical vortices have been reported in the literature: vertical parapets, either solid or porous, placed at the roof edges being the most analysed configuration. Another method for alleviating the high suction peaks due to conical vortices is the use of some non-standard parapet configuration like cantilever parapets. In this paper the influence of roof curvature on the conical vortex pattern appearing on a curved roof (Fig. 1) when subject to oblique winds is experimentally analysed by testing the mean pressure distribution on the curved roofs of low-rise building models in a wind tunnel. Also, the efficiency of cantilever parapets to reduce mean suction loads on curved roofs is experimentally checked. Very high suction loads have been measured on curved roofs, the magnitude of these high suction loads being significantly decreased when cantilever parapets are used. Thus, the suitability of these parapets to reduce wind pressure loads on curved roofs is demonstrated.

Wind tunnel analysis on the influence of cantilever parapets on the wind loads on curved roofs

Different methods to reduce the high suction caused by conical vortices have been reported in the literature: vertical parapets, either solids or porous, placed at the roof edges being the most analysed configuration. Another method for alleviating the high suction peaks due to conical vortices is to round the roof edges. Very recently, the use of some non-standard parapet configurations, like cantilever parapets, has been suggested. In this paper, its efficiency to reduce suction loads on curved roofs is experimentally checked by testing the pressure distribution on the curved roof of a low-rise building model in a wind tunnel. Very high suction loads have been measured on this model, the magnitude of these high suction loads being significantly decreased when cantilever...