Framework for Establishing Limits of Tabular Aerodynamic Models for Flight Dynamics Analysis

This paper describes the use of the Euler equations for the generation and testing of tabular aerodynamic models for flight dynamics analysis. Maneuvers for the AGARD Standard Dynamics Model sharp leading-edge wind-tunnel geometry are considered as a test case. Wind-tunnel data is first used to validate the prediction of static and dynamic coefficients at both low and high angles, featuring complex vortical flow, with good agreement obtained at low to moderate angles of attack. Then the generation of aerodynamic tables is described based on a data fusion approach. Time-optimal maneuvers are generated based on these tables, including level flight trim, pull-ups at constant and varying incidence, and level and 90 degrees turns. The maneuver definition includes the aircraft states and also the control deflections to achieve the motion. The main point of the paper is then to assess the validity of the aerodynamic tables which were used to define the maneuvers. This is done by replaying them, including the control surface motions, through the time accurate computational fluid dynamics code. The resulting forces and moments are compared with the tabular values to assess the presence of inadequately modeled dynamic or unsteady effects. The agreement between the tables and the replay is demonstrated for slow maneuvers. Increasing rate maneuvers show discrepancies which are ascribed to vortical flow hysteresis at the higher rate motions. The framework is suitable for application to more complex viscous flow models, and is powerful for the assessment of the validity of aerodynamics models of the type currently used for studies of flight dynamics.

Comparison of Experimental PIV Data and CFD Simulations for Flow in a Diesel Particulate Filter Inlet Diffuser

Flow maldistribution of the exhaust gas entering a Diesel Particulate Filter (DPF) can cause uneven soot distribution during loading and excessive temperature gradients during the regeneration phase. Minimising the magnitude of this maldistribution is therefore an important consideration in the design of the inlet pipe and diffuser, particularly in situations where packaging constraints dictate bends in the inlet pipe close to the filter, or a sharp diffuser angle. This paper describes the use of Particle Image Velocimetry (PIV) to validate a Computational Fluid Dynamic (CFD) model of the flow within the inlet diffuser of a DPF so that CFD can be used with confidence as a tool to minimise this flow maldistribution. PIV is used to study the flow of gas into a DPF over a range of steady state flow conditions. The distribution of flow approaching the front face of the substrate was of particular interest to this study. Optically clear diffusing cones were designed and placed between pipe and substrate to allow PIV analysis to take place. Stereoscopic PIV was used to eliminate any error produced by the optical aberrations caused by looking through the curved wall of the inlet cone. In parallel to the experiments, numerical analysis was carried out using a CFD program with an incorporated DPF model. Boundary conditions for the CFD simulations were taken from the experimental data, allowing an experimental validation of the numerical results. The CFD model incorporated a DPF model, the cement layers seen in segmented filters and the intumescent matting that is commonly used to pack the filter into a metal casing. The mesh contained approximately 580,000 cells and used the realizable ?-e turbulence model. The CFD simulation predicted both pressure drop across the DPF and the velocity field within the cone and at the DPF face with reasonable accuracy, providing confidence in the use the CFD in future work to design new, more efficient cones.

Numerically Investigating the Effects of Cross-Links in Scaled Microchannel Heat Sinks

Thermal management as a method of heightening performance in miniaturized electronic devices using microchannel heat sinks has recently become of interest to researchers and the industry. One of the current challenges is to design heat sinks with uniform flow distribution. A number of experimental studies have been conducted to seek appropriate designs for microchannel heat sinks. However, pursuing this goal experimentally can be an expensive endeavor. The present work investigates the effect of cross-links on adiabatic two-phase flow in an array of parallel channels. It is carried out using the three dimensional mixture model from the computational fluid dynamics software, FLUENT 6.3. A straight channel and two cross-linked channel models were simulated. The cross-links were located at 1/3 and 2/3 of the channel length, and their widths were one and two times larger than the channel width. All test models had 45 parallel rectangular channels, with a hydraulic diameter of 1.59 mm. The results showed that the trend of flow distribution agrees with experimental results. A new design, with cross-links incorporated, was proposed and the results showed a significant improvement of up to 55% on flow distribution compared with the standard straight channel configuration without a penalty in the pressure drop. Further discussion about the effect of cross-links on flow distribution, flow structure, and pressure drop was also documented.

Effects of Inlet Temperature Uniformity and Nonuniformity on the Tip Leakage Flow and Rotor Blade Tip and Casing Heat Transfer Characteristics

High thermal load appears at the blade tip and casing of a gas turbine engine. It becomes a significant design challenge to protect the turbine materials from this severe situation. As a result of geometric complexity and experimental limitations, computational fluid dynamics tools have been used to predict blade tip leakage flow aerodynamics and heat transfer at typical engine operating conditions. In this paper, the effect of turbine inlet temperature on the tip leakage flow structure and heat transfer has been studied numerically. Uniform low (444 K) and high (800 K) inlet temperatures and nonuniform (parabolic) temperature profiles have been considered at a fixed rotor rotation speed (9500 rpm). The results showed that the change of flow properties at a higher inlet temperature yields significant variations in the leakage flow aerodynamics and heat transfer relative to the lower inlet temperature condition. Aerodynamic behavior of the tip leakage flow varies significantly with the distortion of turbine inlet temperature. For more realistic inlet condition, the velocity range is insignificant at all the time instants. At a high inlet temperature, reverse secondary flow is strongly opposed by the tip leakage flow and the heat transfer fluctuations are reduced greatly.

Transonic Nonlinear Aeroelastic Simulations Using An Harmonic Balance Method

This paper presents an approach to compute transonic Limit Cycle O
scillations using a coupled Harmonic Balance formulation based on the Euler equations for fluid dynamics and finite element models. The paper will investigate the role of aerodynamic (shocks) and structural nonlinearities in driving the limit cycle behaviour. Part icular attention will be given to nonlinear interactions for subcritical LCOs. The Aero elastic Harmonic Balance formulation, allows for solutions of the coupled structural dynamics and CFD system at a reduced cost.

Evaluation of an in situ spatial resolution instrument for fixed beds through the assessment of the invasiveness of probes and a comparison with a micro-kinetic model

This paper reports the detailed description and validation of a fully automated, computer controlled analytical method to spatially probe the gas composition and thermal characteristics in packed bed systems. This method has been designed to limit the invasiveness of the probe, a characteristic assessed using CFD. The thermocouple is aligned with the sampling holes to enable simultaneous recording of the gas composition and temperature profiles. This analysis technique has been validated by studying CO oxidation over a 1% Pt/Al2O3 catalyst. The resultant profiles have been compared with a micro-kinetic model, to further assess the strength of the technique. 

Development and validation of Fluid Structure Interaction Methods for an Oscillating Wave Surge Converter

The term fatigue loads on the Oyster Oscillating Wave Surge Converter (OWSC) is used to describe hydrostatic loads due to water surface elevation with quasi-static changes of state. Therefore a procedure to implement hydrostatic pressure distributions into finite element analysis of the structure is desired. Currently available experimental methods enable one to measure time variant water surface elevation at discrete locations either on or around the body of the scale model during tank tests. This paper discusses the development of a finite element analysis procedure to implement time variant, spatially distributed hydrostatic pressure derived from discretely measured water surface elevation. The developed method can process differently resolved (temporal and spatial) input data and approximate the elevation over the flap faces with user defined properties. The structural loads, namely the forces and moments on the body can then be investigated by post processing the numerical results. This method offers the possibility to process surface elevation or hydrostatic pressure data from computational fluid dynamics simulations and can thus be seen as a first step to a fluid-structure interaction model.

Determination of Non-linear Damping Coefficients of bottom-hinged Oscillating Wave Surge Converters Using Numerical Free Decay Tests

Linear wave theory models are commonly applied to predict the performance of bottom-hinged oscillating wave surge converters (OWSC) in operational sea states. To account for non-linear effects, the additional input of coefficients not included in the model itself becomes necessary. In ocean engineering it is
common practice to obtain damping coefficients of floating structures from free decay tests. This paper presents results obtained from experimental tank tests and numerical computational fluid dynamics simulations of OWSC’s. Agreement between numerical and experimental methods is found to be very good, with CFD providing more data points at small amplitude rotations.
Analysis of the obtained data reveals that linear quadratic-damping, as commonly used in time domain models, is not able to accurately model the occurring damping over the whole regime of rotation amplitudes. The authors
conclude that a hyperbolic function is most suitable to express the instantaneous damping ratio over the rotation amplitude and would be the best choice to be used in coefficient based time domain models.

Automatic extraction of Reduced order models from CFD simulations for building energy modeling

Thermal comfort is defined as “that condition of mind which expresses satisfaction with the thermal environment’ [1] [2]. Field studies have been completed in order to establish the governing conditions for thermal comfort [3]. These studies showed that the internal climate of a room was the strongest factor in establishing thermal comfort. Direct manipulation of the internal climate is necessary to retain an acceptable level of thermal comfort. In order for Building Energy Management Systems (BEMS) strategies to be efficiently utilised it is necessary to have the ability to predict the effect that activating a heating/cooling source (radiators, windows and doors) will have on the room. The numerical modelling of the domain can be challenging due to necessity to capture temperature stratification and/or different heat sources (radiators, computers and human beings). Computational Fluid Dynamic (CFD) models are usually utilised for this function because they provide the level of details required. Although they provide the necessary level of accuracy these models tend to be highly computationally expensive especially when transient behaviour needs to be analysed. Consequently they cannot be integrated in BEMS. This paper presents and describes validation of a CFD-ROM method for real-time simulations of building thermal performance. The CFD-ROM method involves the automatic extraction and solution of reduced order models (ROMs) from validated CFD simulations. The test case used in this work is a room of the Environmental Research Institute (ERI) Building at the University College Cork (UCC). ROMs have shown that they are sufficiently accurate with a total error of less than 1% and successfully retain a satisfactory representation of the phenomena modelled. The number of zones in a ROM defines the size and complexity of that ROM. It has been observed that ROMs with a higher number of zones produce more accurate results. As each ROM has a time to solution of less than 20 seconds they can be integrated into the BEMS of a building which opens the potential to real time physics based building energy modelling.

Prediction of Transonic Limit-Cycle Oscillations Using an Aeroelastic Harmonic Balance Method

This work proposes a novel approach to compute transonic limit-cycle oscillations using high-fidelity analysis. Computational-Fluid-Dynamics based harmonic balance methods have proven to be efficient tools to predict periodic phenomena. This paper’s contribution is to present a new methodology to determine the unknown frequency of oscillations, enabling harmonic balance methods to accurately capture limit-cycle oscillations; this is achieved by defining a frequency-updating procedure based on a coupled computational-fluid-dynamics/computational-structural-dynamics harmonic balance formulation to find the limit-cycle oscillation condition. A pitch/plunge airfoil and delta wing aerodynamic and respective linear structural models are used to validate the new method against conventional time-domain simulations. Results show consistent agreement between the proposed and time-marching methods for both limit-cycle oscillation amplitude and frequency while producing at least a one-order-of-magnitude reduction in computational time.

A Numerical and Experimental Assessment of the Use of a Turbine Utilizing Splitter Blades for an Automotive Variable Geometry Turbocharger

Automotive manufacturers require improved part load engine performance to further improve fuel economy. For a swing vane VGS (Variable Geometry Stator) turbine this means a more closed stator vane, to deal with the low MFRs (Mass Flow Rates), high PRs (Pressure Ratios) and low rotor rotational speeds. During these conditions the turbine is operating at low velocity ratios. As more energy is available at high pressure ratios and during lower turbocharger rotational speeds, a turbine which is efficient at these conditions is desirable. Another key aspect for automotive manufacturers is engine responsiveness. High inertia designs result in “turbo lag” which means an increased time before the target boost pressure is reached. Therefore, designs with improved performance at low velocity ratios, reduced inertia or an increased swallowing capacity are the current targets for turbocharger manufacturers.

To try to meet these design targets a CFD (Computational Fluid Dynamics) study was performed on a turbine wheel using splitter blades. A number of parameters were investigated. These included splitter blade merdional length, blade number and blade angle distribution.

The numerical study was performed on a scaled automotive VGS. Three different stator vane positions have been analysed. A single passage CFD model was developed and used to provide information on the flow features affecting performance in both the stator vanes and turbine.

Following the CFD investigation the design with the best compromise in terms of performance, inertia and increased MFP (Mass Flow Parameter) was selected for manufacture and testing. Tests were performed on a scaled, low temperature turbine test rig. The aerodynamic flow path of the gas stand was the same as that investigated during the CFD. The test results revealed a design which had similar performance at the closed stator vane positions when compared to the baseline wheel. At the maximum MFR stator vane condition a drop of −0.6% pts in efficiency was seen. However, 5.5% increase in MFP was obtained with the additional benefit of a drop in rotor inertia of 3.7%, compared to the baseline wheel.

Coletores de distribuição e sua otimização hidráulica, aplicados a uma instalação de produção de água arrefecida em sistemas de AVAC ( Fórum Picoas)

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Mecânica

Simultaneous Drying and Separation of Composite Agricultural Material for Hay Processing

In composite agricultural materials such as grass, tee, medicinal plants; leaves and stems have a different drying time. By this behavior, after leaving the dryer, the stems may have greater moisture content than desired, while the leaves one minor, which can cause either the appearance of fungi or the collapse of the over-dried material. Taking into account that a lot of grass is dehydrated in forced air dryers, especially rotary drum dryers, this research was developed in order to establish conditions enabling to make a separation of the components during the drying process in order to provide a homogeneous product at the end. For this, a rotary dryer consisting of three concentric cylinders and a circular sieve aligned with the more internal cylinder was proposed; so that, once material enters into the dryer in the area of the inner cylinder, stems pass through sieve to the middle and then continue towards the external cylinder, while the leaves continue by the inner cylinder. For this project, a mixture of Ryegrass and White Clover was used. The characteristics of the components of a mixture were: Drying Rate in thin layer and in rotation, Bulk density, Projected Area, Terminal velocity, weight/Area Ratio, Flux through Rotary sieve. Three drying temperatures; 40°C, 60° C and 80° C, and three rotation speeds; 10 rpm, 20 rpm and 40 rpm were evaluated. It was found that the differences in drying time are the less at 80 °C when the dryer rotates at 40 rpm. Above this speed, the material adheres to the walls of the dryer or sieve and does not flow. According to the measurements of terminal velocity of stems and leaves of the components of the mixture, the speed of the air should be less than 1.5 m s-1 in the inner drum for the leaves and less than 4.5 m s-1 in middle and outer drums for stems, in such way that only the rotational movement of the dryer moves the material and achieves a greater residence time. In other hand, the best rotary sieve separation efficiencies were achieved when the material is dry, but the results are good in all the moisture contents. The best rotary speed of sieve is within the critical rotational speed, i.e. 20 rpm. However, the rotational speed of the dryer, including the sieve in line with the inner cylinder should be 10 rpm or less in order to achieve the greatest residence times of the material inside the dryer and the best agitation through the use of lifting flights. With a finite element analysis of a dryer prototype, using an air flow allowing speeds of air already stated, I was found that the best performance occurs when, through a cover, air enters the dryer front of the Middle cylinder and when the inner cylinder is formed in its entirety through a sieve. This way, air flows in almost equal amounts by both the middle and external cylinders, while part of the air in the Middle cylinder passes through the sieve towards the inner cylinder. With this, leaves do not adhere to the sieve and flow along drier, thanks to the rotating movement of the drums and the showering caused by the lifting flights. In these conditions, the differences in drying time are reduced to 60 minutes, but the residence time is higher for the stems than for leaves, therefore the components of the mixture of grass run out of the dryer with the same desired moisture content.