Unsteady heat transfer analysis of a film cooling flow

Unsteady heat transfer in a turbine blade film cooling flow is studied using detached eddy simulation (DES). Detailed computation of a single row of 35 degree round holes on a flat plate has been obtained for a blowing ratio of 1.0 and a density ratio of 2.0. The instantaneous flow fields and heat transfer distributions are found to be highly unsteady and oscillatory in nature. The fluctuation of the adiabatic effectiveness and heat transfer coefficient, for example, can be as high as 15 and 50 percent of the time-averaged value, respectively. The correlation between the coherent vortical structures and the unsteady heat transfer is carefully examined. It is shown that the fluctuations in the adiabatic effectiveness and heat transfer coefficient are mainly caused by the spanwise fluctuation of the coolant jet and the thermal turbulent boundary layer accompanying the unsteady flow structures.

Flow due to multiple jets downstream of a barrage: Experiments, 3-D CFD and depth-averaged modelling

The flow through and downstream of a row of seven open draft tubes in a barrage has been investigated through laboratory experiments in a wide flume, a three-dimensional (3D) computational fluid dynamics simulation, and a two-dimensional depth-averaged computation. Agreement between the experiments and the 3D modeling is shown to be good, including the prediction of an asymmetric Coandă effect. One aim is to determine the distance downstream at which depth-averaged modeling provides a reasonable prediction; this is shown to be approximately 20 tube diameters downstream of the barrage. Upstream of this, the depth-averaged modeling inaccurately predicts water level, bed shear, and the 3D flow field. The 3D model shows that bed shear stress can be markedly magnified near the barrage, particularly where the jets become attached.

Difference equation approach to two-thermocouple sensor characterization in constant velocity flow environments

Thermocouples are one of the most popular devices for temperature measurement due to their robustness, ease of manufacture and installation, and low cost. However, when used in certain harsh environments, for example, in combustion systems and engine exhausts, large wire diameters are required, and consequently the measurement bandwidth is reduced. This article discusses a software compensation technique to address the loss of high frequency fluctuations based on measurements from two thermocouples. In particular, a difference equation (DE) approach is proposed and compared with existing methods both in simulation and on experimental test rig data with constant flow velocity. It is found that the DE algorithm, combined with the use of generalized total least squares for parameter identification, provides better performance in terms of time constant estimation without any a priori assumption on the time constant ratios of the thermocouples.

Multiphysics Simulation of Thermal Plasma Produced by Lightning Strike on Aircraft Structures

A novel numerical technique is proposed to model thermal plasma of microseconds/milliseconds time-scale effect. Modelling thermal plasma due to lightning strike will allow the estimation of electric current density, plasma pressure, and heat flux at the surface of the aircraft structure. These input data can then be used for better estimation of the mechanical/thermal induced damage on the aircraft structures for better protection systems design. Thermal plasma generated during laser cutting, electric (laser) welding and other plasma processing techniques have been the focus of many researchers. Thermal plasma is a gaseous state that consists from a mixture of electrons, ions, and natural particles. Thermal plasma can be assumed to be in local thermodynamic equilibrium, which means the electrons and the heavy species have equal temperature. Different numerical techniques have been developed using a coupled Navier Stokes – Heat transfer – Electromagnetic equations based on the assumption that the thermal plasma is a single laminar gas flow. These previous efforts focused on generating thermal plasma of time-scale in the range of seconds. Lighting strike on aircraft structures generates thermal plasma of time-scale of milliseconds/microseconds, which makes the previous physics used not applicable. The difficulty comes from the Navier-Stokes equations as the fluid is simulated under shock load, this introducing significant changes in the density and temperature of the fluid.

Python PSS/E Simulation to Test Efficacy of Proposed PMU Based WAMS and Potential WAMPAC Applications

A PSS/E 32 model of a real section of the Northern Ireland electrical grid was dynamically controlled with Python 2.5. In this manner data from a proposed wide area monitoring system was simulated. The area is of interest as it is a weakly coupled distribution grid with significant distributed generation. The data was used to create an optimization and protection metric that reflected reactive power flow, voltage profile, thermal overload and voltage excursions. Step changes in the metric were introduced upon the operation of special protection systems and voltage excursions. A wide variety of grid conditions were simulated while tap changer positions and switched capacitor banks were iterated through; with the most desirable state returning the lowest optimization and protection metric. The optimized metric was compared against the metric generated from the standard system state returned by PSS/E. Various grid scenarios were explored involving an intact network and compromised networks (line loss) under summer maximum, summer minimum and winter maximum conditions. In each instance the output from the installed distributed generation is varied between 0 MW and 80 MW (120% of installed capacity). It is shown that in grid models the triggering of special protection systems is delayed by between 1 MW and 6 MW (1.5% to 9% of capacity), with 3.5 MW being the average. The optimization and protection metric gives a quantitative value for system health and demonstrates the potential efficacy of wide area monitoring for protection and control.

Simulation using novel equipment designed to explain spirometric abnormalities in respiratory disease enhances learning in higher cognitive domains

Simulation of disorders of respiratory mechanics shown by spirometry provides insight into the pathophysiology of disease but some clinically important disorders have not been simulated and none have been formally evaluated for education. We have designed simple mechanical devices which, along with existing simulators, enable all the main dysfunctions which have diagnostic value in spirometry to be simulated and clearly explained with visual and haptic feedback. We modelled the airways as Starling resistors by a clearly visible mechanical action to simulate intra- and extra-thoracic obstruction. A narrow tube was used to simulate fixed large airway obstruction and inelastic bands to simulate restriction. We hypothesized that using simulators whose action explains disease promotes learning especially in higher domain educational objectives. The main features of obstruction and restriction were correctly simulated. Simulation of variable extra-thoracic obstruction caused blunting and plateauing of inspiratory flow, and simulation of intra-thoracic obstruction caused limitation of expiratory flow with marked dynamic compression. Multiple choice tests were created with questions allocated to lower (remember and understand) or higher cognitive domains (apply, analyse and evaluate). In a cross-over design, overall mean scores increased after 1½ h simulation spirometry (43-68 %, effect size 1.06, P < 0.0001). In higher cognitive domains the mean score was lower before and increased further than lower domains (Δ 30 vs 20 %, higher vs lower effect size 0.22, P < 0.05). In conclusion, the devices successfully simulate various patterns of obstruction and restriction. Using these devices medical students achieved marked enhancement of learning especially in higher cognitive domains.

Measuring and modelling air mass flow rate in the injection stretch blow moulding process

The injection stretch blow moulding process involves the inflation and stretching of a hot preform into a mould to form bottles. A critical process variable and an essential input for process simulations is the rate of pressure increase within the preform during forming, which is regulated by an air flow restrictor valve. The paper describes a set of experiments for measuring the air flow rate within an industrial ISBM machine and the subsequent modelling of it with the FEA package ABAQUS. Two rigid containers were inserted into a Sidel SBO1 blow moulding machine and subjected to different supply pressures and air flow restrictor settings. The pressure and air temperature were recorded for each experiment enabling the mass flow rate of air to be determined along with an important machine characteristic known as the ‘dead volume’. The experimental setup was simulated within the commercial FEA package ABAQUS/Explicit using a combination of structural, fluid and fluid link elements that idealize the air flowing through an orifice behaving as an ideal gas under isothermal conditions. Results between experiment and simulation are compared and show a good correlation.

Shocks in unmagnetized plasma with a shear flow: Stability and magnetic field generation

A pair of curved shocks in a collisionless plasma is examined with a two-dimensional particle-in-cell simulation. The shocks are created by the collision of two electron-ion clouds at a speed that exceeds everywhere the threshold speed for shock formation. A variation of the collision speed along the initially planar collision boundary, which is comparable to the ion acoustic speed, yields a curvature of the shock that increases with time. The spatially varying Mach number of the shocks results in a variation of the downstream density in the direction along the shock boundary. This variation is eventually equilibrated by the thermal diffusion of ions. The pair of shocks is stable for tens of inverse ion plasma frequencies. The angle between the mean flow velocity vector of the inflowing upstream plasma and the shock's electrostatic field increases steadily during this time. The disalignment of both vectors gives rise to a rotational electron flow, which yields the growth of magnetic field patches that are coherent over tens of electron skin depths.

Modelling and visualisation of material flow in friction stir spot welding

Abstract The material flow in friction stir spot welding of aluminium to both aluminium and steel has been investigated, using pinless tools in a lap joint geometry. The flow behaviour was revealed experimentally using dissimilar Al alloys of similar strength. The effect on the material flow of tool surface features, welding conditions (rotation speed, plunge depth, dwell time), and the surface state of the steel sheet (un-coated or galvanized) have been systematically studied. A novel kinematic flow model is presented, which successfully predicts the observed layering of the dissimilar Al alloys under a range of conditions. The model and the experimental observations provide a consistent interpretation of the stick-slip conditions at the tool-workpiece interface, addressing an elusive and long-standing issue in the modelling of heat generation in friction stir processing.

A simulation investigation into the effects of sea level rise on tidal propagation in fjords

Senior thesis written for Oceanography 445

Simulation model for driving dynamics, energy use and power supply

Dissertação para obtenção do grau de Mestre em Engenharia Electrotécnica Ramo de Energia

Simulation and dynamics of freeway traffic

This paper presents the new package entitled Simulator of Intelligent Transportation Systems (SITS) and a computational oriented analysis of traffic dynamics. The SITS adopts a microscopic simulation approach to reproduce real traffic conditions considering different types of vehicles, drivers and roads. A set of experiments with the SITS reveal the dynamic phenomena exhibited by this kind of system. For this purpose a modelling formalism is developed that embeds the statistics and the Laplace transform. The results make possible the adoption of classical system theory tools and point out that it is possible to study traffic systems taking advantage of the knowledge gathered with automatic control algorithms. A complementary perspective for the analysis of the traffic flow is also quantified through the entropy measure.

Project, modelling and simulation of air stripping columns for the treatment of groundwater contaminated with volatile organic compounds

Volatile organic compounds are a common source of groundwater contamination that can be easily removed by air stripping in columns with random packing and using a counter-current flow between the phases. This work proposes a new methodology for the column design for any particular type of packing and contaminant avoiding the necessity of a pre-defined diameter used in the classical approach. It also renders unnecessary the employment of the graphical Eckert generalized correlation for pressure drop estimates. The hydraulic features are previously chosen as a project criterion and only afterwards the mass transfer phenomena are incorporated, in opposition to conventional approach. The design procedure was translated into a convenient algorithm using C++ as programming language. A column was built in order to test the models used either in the design or in the simulation of the column performance. The experiments were fulfilled using a solution of chloroform in distilled water. Another model was built to simulate the operational performance of the column, both in steady state and in transient conditions. It consists in a system of two partial non linear differential equations (distributed parameters). Nevertheless, when flows are steady, the system became linear, although there is not an evident solution in analytical terms. In steady state the resulting system of ODE can be solved, allowing for the calculation of the concentration profile in both phases inside the column. In transient state the system of PDE was numerically solved by finite differences, after a previous linearization.

Microsatellites can be misleading: an empirical and simulation study.

It has been long recognized that highly polymorphic genetic markers can lead to underestimation of divergence between populations when migration is low. Microsatellite loci, which are characterized by extremely high mutation rates, are particularly likely to be affected. Here, we report genetic differentiation estimates in a contact zone between two chromosome races of the common shrew (Sorex araneus), based on 10 autosomal microsatellites, a newly developed Y-chromosome microsatellite, and mitochondrial DNA. These results are compared to previous data on proteins and karyotypes. Estimates of genetic differentiation based on F- and R-statistics are much lower for autosomal microsatellites than for all other genetic markers. We show by simulations that this discrepancy stems mainly from the high mutation rate of microsatellite markers for F-statistics and from deviations from a single-step mutation model for R-statistics. The sex-linked genetic markers show that all gene exchange between races is mediated by females. The absence of male-mediated gene flow most likely results from male hybrid sterility.

Une nouvelle méthode smoothed particle hydrodynamics : simulation des interfaces immergées et de la dynamique Brownienne des molécules avec des interactions hydrodynamiques

Dans cette thèse, nous présentons une nouvelle méthode smoothed particle hydrodynamics (SPH) pour la résolution des équations de Navier-Stokes incompressibles, même en présence des forces singulières. Les termes de sources singulières sont traités d'une manière similaire à celle que l'on retrouve dans la méthode Immersed Boundary (IB) de Peskin (2002) ou de la méthode régularisée de Stokeslets (Cortez, 2001). Dans notre schéma numérique, nous mettons en oeuvre une méthode de projection sans pression de second ordre inspirée de Kim et Moin (1985). Ce schéma évite complètement les difficultés qui peuvent être rencontrées avec la prescription des conditions aux frontières de Neumann sur la pression. Nous présentons deux variantes de cette approche: l'une, Lagrangienne, qui est communément utilisée et l'autre, Eulerienne, car nous considérons simplement que les particules SPH sont des points de quadrature où les propriétés du fluide sont calculées, donc, ces points peuvent être laissés fixes dans le temps. Notre méthode SPH est d'abord testée à la résolution du problème de Poiseuille bidimensionnel entre deux plaques infinies et nous effectuons une analyse détaillée de l'erreur des calculs. Pour ce problème, les résultats sont similaires autant lorsque les particules SPH sont libres de se déplacer que lorsqu'elles sont fixes. Nous traitons, par ailleurs, du problème de la dynamique d'une membrane immergée dans un fluide visqueux et incompressible avec notre méthode SPH. La membrane est représentée par une spline cubique le long de laquelle la tension présente dans la membrane est calculée et transmise au fluide environnant. Les équations de Navier-Stokes, avec une force singulière issue de la membrane sont ensuite résolues pour déterminer la vitesse du fluide dans lequel est immergée la membrane. La vitesse du fluide, ainsi obtenue, est interpolée sur l'interface, afin de déterminer son déplacement. Nous discutons des avantages à maintenir les particules SPH fixes au lieu de les laisser libres de se déplacer. Nous appliquons ensuite notre méthode SPH à la simulation des écoulements confinés des solutions de polymères non dilués avec une interaction hydrodynamique et des forces d'exclusion de volume. Le point de départ de l'algorithme est le système couplé des équations de Langevin pour les polymères et le solvant (CLEPS) (voir par exemple Oono et Freed (1981) et Öttinger et Rabin (1989)) décrivant, dans le cas présent, les dynamiques microscopiques d'une solution de polymère en écoulement avec une représentation bille-ressort des macromolécules. Des tests numériques de certains écoulements dans des canaux bidimensionnels révèlent que l'utilisation de la méthode de projection d'ordre deux couplée à des points de quadrature SPH fixes conduit à un ordre de convergence de la vitesse qui est de deux et à une convergence d'ordre sensiblement égale à deux pour la pression, pourvu que la solution soit suffisamment lisse. Dans le cas des calculs à grandes échelles pour les altères et pour les chaînes de bille-ressort, un choix approprié du nombre de particules SPH en fonction du nombre des billes N permet, en l'absence des forces d'exclusion de volume, de montrer que le coût de notre algorithme est d'ordre O(N). Enfin, nous amorçons des calculs tridimensionnels avec notre modèle SPH. Dans cette optique, nous résolvons le problème de l'écoulement de Poiseuille tridimensionnel entre deux plaques parallèles infinies et le problème de l'écoulement de Poiseuille dans une conduite rectangulaire infiniment longue. De plus, nous simulons en dimension trois des écoulements confinés entre deux plaques infinies des solutions de polymères non diluées avec une interaction hydrodynamique et des forces d'exclusion de volume.