Numerical simulation of the flow and performance evaluation of an aerospike engine

The main focus of this work is to define a numerical methodology to simulate an aerospike engine and then to analyse the performance of DemoP1, which is a small aerospike demonstrator built by Pangea Aerospace. The aerospike is a promising solution to build more efficient engine than the actual one. Its main advantage is the expansion adaptation that allows to reach the optimal expansion in a wide range of ambient pressures delivering more thrust than an equivalent bell-shaped nozzle. The main drawbacks are the cooling system design and the spike manufacturing but nowadays, these issues seem to be overcome with the use of the additive manufacturing method. The simulations are performed with dbnsTurbFoam which is a solver of OpenFOAM. It has been designed to simulate a supersonic compressible turbulent flow. This work is divided in four chapters. The first one is a short introduction. The second one shows a brief summary of the theoretical performance of the aerospike. The third one introduces the numerical methodology to simulate a compressible supersonic flow. In the fourth chapter, the solver has been verified with an experiment found in literature. And in the fifth chapter, the simulations on DemoP1 engine are illustrated.

Uma revisão sobre o efeito Toms: o fenômeno onde macromoléculas atenuam a turbulência em um líquido

It is very well known that the addition of polymers to a liquid increases the shear viscosity of the solution. In other words, the polymer increases the dissipation of the flow energy. Contrarily, in turbulent flow, some particular macromolecules in very low concentration are able to produce large attenuation in the turbulence and thus, decreasing the dissipation of the energy. This article present a brief revision about macroscopic and molecular models used to explain this dynamic effect. Some of the experimental techniques used to quantify the attenuation of the turbulence and the main active substances are also discussed.

LOCAL HEAD LOSS OF NON-COAXIAL EMITTERS INSERTED IN POLYETHYLENE PIPE

The design of a lateral line for drip irrigation requires accurate evaluation of head losses in not only the pipe but in the emitters as well. A procedure was developed to determine localized head losses within the emitters by the formulation of a mathematical model that accounts for the obstruction caused by the insertion point. These localized losses can be significant when compared with tire total head losses within the system due to the large number of emitters typically installed along the lateral line. Air experiment was carried out by altering flow characteristics to create Reynolds numbers (R) from 7,480 to 32,597 to provide turbulent flow and a maximum velocity of 2.0 m s(-1). The geometry of the emitter was determined by an optical projector and sensor An equation was formulated to facilitate the localized head loss calculation using the geometric characteristics of the emitter (emitter length, obstruction ratio, and contraction coefficient). The mathematical model was tested using laboratory measurements on four emitters. The local head loss was accurately estimated for the Uniram (difference of +13.6%) and Drip Net (difference of +7.7%) emitters, while appreciable deviations were found for the Twin Plus (-21.8%) and Tiran (+50%) emitters. The head loss estimated by the model was sensitive to the variations in the obstruction area of the emitter However, the variations in the local head loss did not result in significant variations in the maximum length of the lateral lines. In general, for all the analyzed emitters, a 50% increase in the local head loss for the emitters resulted in less than an 8% reduction in the maximum lateral length.

Friction factor for small diameter polyethylene pipes

On this paper, the results of an experimental study oil the hydraulic friction loss for small-diameter polyethylene pipes are reported. The experiment was carried out using a range of Reynolds number between 6000 to 72000, obtained by varying discharge at 20 degrees C water temperature, with internal pipe diameters of 10.0 mm, 12.9 mm, 16.1 mm, 17.4 mm and 19.7 mm. According to the analysis results and experimental conditions, the friction factor 0 of the Darcy-Weisbach equation call be estimated with c = 0.300 and m = 0.25. The Blasius equation (c = 0.316 and m = 0.25) gives an overestimate of friction loss, although this fact is non-restrictive for micro-irrigation system designs. The analysis shows that both the Blasius and the adjusted equation parameters allow for accurate friction factor estimates, characterized by low mean error (5.1%).

Antioxidant enzymatic defense in salivary glands of streptozotocin-induced diabetic rats: a temporal study

Hyperglycemia induces overproduction of superoxide and it is related to diabetic complications. In this study, we analyzed the antioxidant enzymatic defense and the lipid peroxidation of rat salivary glands in six different periods of diabetic condition. Ninety-six rats were divided into 12 groups: C7/14/21128/45/60 (non-diabetic animals) and D7/14/21/28/45/60 (diabetic animals). Diabetes was induced by streptozotocin and the rats were euthanized after 7, 14, 21, 28, 45, or 60 days. Their parotid (PA) and submandibular (SM) glands were removed soon after the sacrifice and the total protein and malondialdehyde (MDA) concentrations, as well as, the superoxide dismutase (SOD), glutathione peroxidase (GPx), and catalase (CAT) activities were determined. Twenty-one days after the diabetes induction, the SM glands showed an increase in SOD, CAT, and GPx activities, as well as, MDA concentration. Concerning the PA glands, an increase in the CAT activity and MDA content was observed throughout the observation period. The results suggest that diabetes can cause alterations on the salivary glands and that PA and SM glands react differently when exposed to diabetes condition. However, no impairment of antioxidant system was observed in the group whose diabetic condition had been induced 60 days earlier, herein named 60-day group. Copyright (C) 2010 John Wiley & Sons, Ltd.

Efeito da razão de submersão sobre os parâmetros hidrodinâmicos de uma coluna de recirculação de líquido através de um tubo exterior

No trabalho realizado nesta tese, procedeu-se ao estudo hidrodinâmico de uma coluna de borbulhamento de recirculação externa do líquido (CREL), que permitiu ampliar o conhecimento já existente sobre este tipo de colunas. Para realizar este estudo utilizaram-se líquidos viscosos Newtonianos, nomeadamente soluções aquosas de glicerina com viscosidades entre 0,007 e 0,522 Pa.s. Para a gama de caudais de ar injectados, 1,5x10-5 até 1,35x10-4 m3/s, o ar ascendia ao longo da coluna de borbulhamento sob a forma de bolhas tubulares. Após a realização dos ensaios, verificou-se que o regime de escoamento do líquido entre bolhas tubulares variou desde o tipo laminar até transição (1,9 turbulento (12,1 < ReV <3297). Nos ensaios experimentais foram utilizadas duas colunas de borbulhamento com diferentes diâmetros internos de 22 e 32 mm. A razão de submersão, , seleccionada foi de 0,85. Este valor foi escolhido por ser significativamente diferente dos usados em estudos anteriores (0,935 por Santos (2004), Santos (2005) e 0,91 por Rangel (2008)) de forma a poder analisar a influência da razão de submersão sobre os parâmetros hidrodinâmicos. Nos ensaios realizados foram determinadas a retenção do gás na coluna e a velocidade superficial do líquido em escoamento no interior da coluna de borbulhamento. Os valores da retenção do gás foram comparados com valores teóricos previstos por equações empíricas de Nicklin et al (1962). Estes valores também foram comparados com valores obtidos experimentalmente em trabalhos anteriores, onde se utilizaram razões de submersão de 0,935 e 0,91. Após análise dos valores da retenção verificou-se que estes eram mais elevados para as soluções de maior viscosidade e menor diâmetro da coluna. Para as experiências realizadas com as soluções com viscosidades entre 0,52 Pa.s e 0,053 Pa.s, os valores experimentais eram inferiores aos previstos por Nicklin et al (1962) para regime laminar tanto na coluna de 22 mm como na coluna de 32 mm. Para a solução com 0,007 Pa.s de viscosidade os valores experimentais situavam-se na zona de transição entre o regime laminar e o regime turbulento em ambas as colunas. Verificou-se que a diminuição da razão de submersão de 0,935 para 0,85 aumentou significativamente a retenção em todas as experiências realizadas. A velocidade superficial média de circulação do líquido variava na razão inversa do diâmetro e da viscosidade. A diminuição da retenção provocou uma diminuição significativa na velocidade do líquido.

Multiple Mechanical Complications in ST-Segment Elevation Myocardial Infarction with Angiographically Normal Coronary Arteries

This case report discusses an unusual presentation of ST-segment elevation myocardial infarction (STEMI) with normal coronary arteries and severe mechanical complications successfully treated with surgery. An 82-year-old man presented STEMI with angiographically normal coronary arteries and no major echocardiographic alterations at discharge. At the first month follow-up, he complained of fatigue and dyspnea, and contrast echocardiography complemented by cardiac magnetic resonance imaging revealed a large left ventricular apical aneurysm with a thrombus communicating by two jets of a turbulent flow to an aneurysmatic formation of the right ventricular apex. The patient underwent a Dor procedure, which was successful. Ventricular septal defects and ventricular aneurysms are rare but devastating complications of STEMI, with almost all patients presenting multivessel coronary artery disease. Interestingly in this case, the angiographic pattern was normal.

Runoff and sediment transport in a degraded area

Gully erosion occurs by the combined action of splash, sheetwash and rill-wash (interrill and rill erosion). These erosion processes have a great capacity for both sediment production and sediment transport. The objectives of this experiment were to evaluate hydrological and sediment transport in a degraded area, severely dissected by gullies; to assess the hydraulic flow characteristics and their aggregate transport capacity; and to measure the initial splash erosion rate. In the study area in Guarapuava, State of Paraná, Brazil (lat 25º 24' S; long 51º24' W; 1034 m asl), the soil was classified as Cambissolo Húmico alumínico, with the following particle-size composition: sand 0.116 kg kg-1; silt 0.180 kg kg-1; and clay 0.704 kg kg-1. The approach of this research was based on microcatchments formed in the ground, to study the hydrological response and sediment transport. A total of eight rill systems were simulated with dry and wet soil. An average rainfall of 33.7 ± 4.0 mm was produced for 35 to 54 min by a rainfall simulator. The equipment was installed, and a trough was placed at the end of the rill to collect sediments and water. During the simulation, the following variables were measured: time to runoff, time to ponding, time of recession, flow velocity, depth, ratio of the initial splash and grain size. The rainsplash of dry topsoil was more than twice as high as under moist conditions (5 g m-2 min-1 and 2 g m-2 min-1, respectively). The characteristics of the flow hydraulics indicate transition from laminar to turbulent flow [Re (Reynolds number) 1000-2000]. In addition, it was observed that a flow velocity of 0.12 m s-1 was the threshold for turbulent flow (Re > 2000), especially at the end of the rainfall simulation. The rill flow tended to be subcritical [Fr (Froude Number) < 1.0]. The variation in hydrological attributes (infiltration and runoff) was lower, while the sediment yield was variable. The erosion in the rill systems was characterized as limited transport, although the degraded area generated an average of 394 g m-2 of sediment in each simulation.

Inertial effects on reactive particles advected by turbulence

We study the problem of the advection of passive particles with inertia in a two-dimensional, synthetic, and stationary turbulent flow. The asymptotic analytical result and numerical simulations show the importance of inertial bias in collecting the particles preferentially in certain regions of the flow, depending on their density relative to that of the flow. We also study how these aggregates are affected when a simple chemical reaction mechanism is introduced through a Eulerian scheme. We find that inertia can be responsible for maintaining a stationary concentration pattern even under nonfavorable reactive conditions or destroying it under favorable ones.

Reaction-diffusion fronts under stochastic advection

We study front propagation in stirred media using a simplified modelization of the turbulent flow. Computer simulations reveal the existence of the two limiting propagation modes observed in recent experiments with liquid phase isothermal reactions. These two modes respectively correspond to a wrinkled although sharp propagating interface and to a broadened one. Specific laws relative to the enhancement of the front velocity in each regime are confirmed by our simulations.

Inertial effects on reactive particles advected by turbulence

We study the problem of the advection of passive particles with inertia in a two-dimensional, synthetic, and stationary turbulent flow. The asymptotic analytical result and numerical simulations show the importance of inertial bias in collecting the particles preferentially in certain regions of the flow, depending on their density relative to that of the flow. We also study how these aggregates are affected when a simple chemical reaction mechanism is introduced through a Eulerian scheme. We find that inertia can be responsible for maintaining a stationary concentration pattern even under nonfavorable reactive conditions or destroying it under favorable ones.

Reaction-diffusion fronts under stochastic advection

We study front propagation in stirred media using a simplified modelization of the turbulent flow. Computer simulations reveal the existence of the two limiting propagation modes observed in recent experiments with liquid phase isothermal reactions. These two modes respectively correspond to a wrinkled although sharp propagating interface and to a broadened one. Specific laws relative to the enhancement of the front velocity in each regime are confirmed by our simulations.

Experimental insights into deformation dynamics and intermittency in rapid granular shear flows

This work concerns the experimental study of rapid granular shear flows in annular Couette geometry. The flow is induced by continuous driving of the horizontal plate at the top of the granular bed in an annulus. The compressive pressure, driving torque, instantaneous bed height and rotational speed of the shearing plate are measured. Moreover, local stress fluctuations are measured in a medium made of steel spheres 2 and 3 mm in diameter. Both monodisperse packing and bidisperse packing are investigated to reveal the influence of size diversity in intermittent features of granular materials. Experiments are conducted in an annulus that can contain up to 15 kg of spherical steel balls. The shearing granular medium takes place via the rotation of the upper plate which compresses the material loaded inside the annulus. Fluctuations of compressive force are locally measured at the bottom of the annulus using a piezoelectric sensor. Rapid shear flow experiments are pursued at different compressive forces and shear rates and the sensitivity of fluctuations are then investigated by different means through monodisperse and bidisperse packings. Another important feature of rapid granular shear flows is the formation of ordered structures upon shearing. It requires a certain range for the amount of granular material (uniform size distribution) loaded in the system in order to obtain stable flows. This is studied more deeply in this thesis. The results of the current work bring some new insights into deformation dynamics and intermittency in rapid granular shear flows. The experimental apparatus is modified in comparison to earlier investigations. The measurements produce data for various quantities continuously sampled from the start of shearing to the end. Static failure and dynamic shearing ofa granular medium is investigated. The results of this work revealed some important features of failure dynamics and structure formation in the system. Furthermore, some computer simulations are performed in a 2D annulus to examine the nature of kinetic energy dissipation. It is found that turbulent flow models can statistically represent rapid granular flows with high accuracy. In addition to academic outcomes and scientific publications our results have a number of technological applications associated with grinding, mining and massive grain storages.

Quadrant/octant sequencing and the role of coherent structures in bed load sediment entrainment

To permit the tracking of turbulent flow structures in an Eulerian frame from single-point measurements, we make use of a generalization of conventional two-dimensional quadrant analysis to three-dimensional octants. We characterize flow structures using the sequences of these octants and show how significance may be attached to particular sequences using statistical mull models. We analyze an example experiment and show how a particular dominant flow structure can be identified from the conditional probability of octant sequences. The frequency of this structure corresponds to the dominant peak in the velocity spectra and exerts a high proportion of the total shear stress. We link this structure explicitly to the propensity for sediment entrainment and show that greater insight into sediment entrainment can be obtained by disaggregating those octants that occur within the identified macroturbulence structure from those that do not. Hence, this work goes beyond critiques of Reynolds stress approaches to bed load entrainment that highlight the importance of outward interactions, to identifying and prioritizing the quadrants/octants that define particular flow structures. Key Points <list list-type=''bulleted'' id=''jgrf20196-list-0001''> <list-item id=''jgrf20196-li-0001''>A new method for analysing single point velocity data is presented <list-item id=''jgrf20196-li-0002''>Flow structures are identified by a sequence of flow states (termed octants) <list-item id=''jgrf20196-li-0003''>The identified structure exerts high stresses and causes bed-load entrainment

Lämmönsiirron numeerinen laskenta profiloidulle kanavalle.

Työn tavoitteena oli tutkia aaltomaisen profiloinnin vaikutusta suorien jäähdytyskanavien lämmönsiirtoon ja painehäviöön. Erilaisia profiileja oli kymmenen kappaletta ja ne olivat 5mm leveitä ja 30cm pitkiä kukin. Ne laskettiin kolmeulotteisina tapauksina FINFLO-virtausratkaisijalla kolmella eri Reynoldsin luvulla, jotka vastasivat laminaarista, osittain turbulenttista ja lähes kokonaan turbulenttista virtausta. Lämmönsiirtoaine oli kuiva +30°C ilma ja profiloinnin toteutustapa oli toisiaan sivuavat ympyräkaaret kolmella erilaisella säteen arvolla ja kolmella erilaisella aallonpituuden arvolla. Lisäksi laskettiin saman levyisen tasokanavan arvot jokaisella Reynoldsin luvulla kaksiulotteisina tapauksina. Näitä profiloimattomia kanavia pidettiin referenssitapauksina. Tuloksena havaittiin että profiloimalla saadaan yksiselitteisesti suurempi lämpöteho ulos samasta tilavuudesta. Lämmönsiirtokerroin kasvaa profiloinnin avulla parhaimmillaan n. 20% käytetystä turbulenssimallista tai lämmönsiirtokertoimen määritelmästä riippumatta. Painehäviö kasvaa myös aina, mutta kitkakerroin voi hieman pienentyä. Profiilin varsinaisena hyvyyskriteerinä pidettiin lämmönsiirtokertoimen ja kitkakertoimen suhdetta h/f. Se osoittautui riippuvaksi Reynoldsin luvusta ja turbulenssimallista; ASM ja Chien k-έ -mallit ennustavat transitioetäisyyden eri tavalla. Laminaarisilla virtauksilla h/f :n vaihtelu oli vähäistä; suhde vaihteli vain ±5% eri profiilien kesken. ASM-mallilla havaittiin sekundääripyörteilyä, ehkä siksi että se mallintaa anisotrooppisen turbulenssin. Chien k-έ malli ennusti suuremman ja aikaisemmin alkavan turbulenttisuuden kuin ASM. Lisäksi havaittiin mm. että tietyillä profiileilla muodostuu kanavan kapeimpaan kohtaan selvä virtausnopeuden paikallinen minimi seinämän läheisyyden takia.