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.

Effects of oscillatory flow on fertilization in the green sea urchin Strongylocentrotus droebachiensis

Broadcast spawning invertebrates that live in shallow, high-energy coastal habitats are subjected to oscillatory water motion that creates unsteady flow fields above the surface of animals. The frequency of the oscillatory fluctuations is driven by the wave period, which will influence the stability of local flow structures and may affect fertilization processes. Using an oscillatory water tunnel, we quantified the percentage of eggs fertilized on or near spawning green sea urchins, Strongylocentrotus droebachiensis. Eggs were sampled in the water column, wake eddy, substratum and aboral surface under a range of different periods (T = 4.5 – 12.7 s) and velocities of oscillatory flow. The root-mean-square wave velocity (rms(uw)) was a good predictor of fertilization in oscillatory flow, although the root-mean-square of total velocity (rms(u)), which incorporates all the components of flow (current, wave and turbulence), also provided significant predictions. The percentage of eggs fertilized varied between 50 – 85% at low flows (rms(uw) < 0.02 m s-1), depending on the location sampled, but declined to below 10% for most locations at higher rms(uw). The water column was an important location for fertilization with a relative contribution greater than that of the aboral surface, especially at medium and high rms(uw) categories. We conclude that gametes can be successfully fertilized on or near the parent under a range of oscillatory flow conditions.

Numerical Methods for the Unsteady Compressible Navier-Stokes Equations

We consider numerical methods for the compressible time dependent Navier-Stokes equations, discussing the spatial discretization by Finite Volume and Discontinuous Galerkin methods, the time integration by time adaptive implicit Runge-Kutta and Rosenbrock methods and the solution of the appearing nonlinear and linear equations systems by preconditioned Jacobian-Free Newton-Krylov, as well as Multigrid methods. As applications, thermal Fluid structure interaction and other unsteady flow problems are considered. The text is aimed at both mathematicians and engineers.

Transient flow of the oil-refrigerant mixture through the radial clearance in rolling piston compressors

Unsteady flow of oil and refrigerant gas through radial clearance in rolling piston compressors has been modeled as a heterogeneous mixture, where the properties are determined from the species conservation transport equation coupled with momentum and energy equations. Time variations of pressure, tangential velocity of the rolling piston and radial clearance due to pump setting have been included in the mixture flow model. Those variables have been obtained by modeling the compression process, rolling piston dynamics and by using geometric characteristics of the pump, respectively. An important conclusion concerning this work is the large variation of refrigerant concentration in the oil-filled radial clearance during the compression cycle. That is particularly true for large values of mass flow rates, and for those cases the flow mixture cannot be considered as having uniform concentration. In presence of low mass flow rates homogeneous flow prevail and the mixture tend to have a uniform concentration. In general, it was observed that for calculating the refrigerant mass flow rate using the difference in refrigerant concentration between compression and suction chambers, a time average value for the gas concentration should be used at the clearance inlet.

Study on transient ice formation of laminar flow inside externally supercooled rectangular duct

The transient process of solidification of laminar liquid flow (water) submitted to super-cooling was investigated both theoretically and experimentally. In this study an alternative analytical formulation and numerical approach were adopted resulting in the unsteady model with temperature dependent thermophysical properties in the solid region. The proposed model is based upon the fundamental equations of energy balance in the solid and liquid regions as well as across the solidification front. The basic equations and the associated boundary and initial conditions were made dimensionless by using the Landau transformation to immobilize the moving front and render the problem to a fixed plane type problem. A laminar velocity profile is admitted in the liquid domain and the resulting equations were discretized using the finite difference approach. The numerical predictions obtained were compared with the available results based on other models and concepts such as Neumann analytical model, the apparent thermal capacity model due to Bonacina and the conventional fixed grid energy model due to Goodrich. To obtain further comparisons and more validation of the model and the numerical solution, an experimental rig was constructed and instrumented permitting very well controlled experimental measurements. The numerical predictions were compared with the experimental results and the agreement was found satisfactory.

Aeroelastic behavior of a wing with a magnetorheological hysteretic model to the damping

This paper shows the application of a hysteretic model for the Magnetorheological Damper (MRD) placed in the plunge degree-of-freedom of aeroelastic model of a wing. This hysteretic MRD model was developed by the researchers of the French Aerospace Lab. (ONERA) and describe, with a very good precision, the hysteretic behavior of the MRD. The aeroelastic model used in this paper do not have structural nonlinearities, the only nonlinearities showed in the model, are in the unsteady flow equations and are the same proposed by Theodorsen and Wagner in their unsteady aerodynamics theory; and the nonlinearity introduced by the hysteretic model used. The main objective of this paper is show the mathematical modeling of the problem and the equations that describes the aeroelastic response of our problem; and the gain obtained with the introduction of this hysteretic model in the equations with respect to other models that do not show the this behavior, through of pictures that represents the time response and Phase diagrams. These pictures are obtained using flow velocities before and after the flutter velocity. Finally, an open-loop control was made to show the effect of the MRD in the aeroelastic behavior.

ANALYSIS AND VISUALIZATION OF FLOW FIELDS USING INFORMATION-THEORETIC TECHNIQUES AND GRAPH-BASED REPRESENTATIONS

Three-dimensional flow visualization plays an essential role in many areas of science and engineering, such as aero- and hydro-dynamical systems which dominate various physical and natural phenomena. For popular methods such as the streamline visualization to be effective, they should capture the underlying flow features while facilitating user observation and understanding of the flow field in a clear manner. My research mainly focuses on the analysis and visualization of flow fields using various techniques, e.g. information-theoretic techniques and graph-based representations. Since the streamline visualization is a popular technique in flow field visualization, how to select good streamlines to capture flow patterns and how to pick good viewpoints to observe flow fields become critical. We treat streamline selection and viewpoint selection as symmetric problems and solve them simultaneously using the dual information channel [81]. To the best of my knowledge, this is the first attempt in flow visualization to combine these two selection problems in a unified approach. This work selects streamline in a view-independent manner and the selected streamlines will not change for all viewpoints. My another work [56] uses an information-theoretic approach to evaluate the importance of each streamline under various sample viewpoints and presents a solution for view-dependent streamline selection that guarantees coherent streamline update when the view changes gradually. When projecting 3D streamlines to 2D images for viewing, occlusion and clutter become inevitable. To address this challenge, we design FlowGraph [57, 58], a novel compound graph representation that organizes field line clusters and spatiotemporal regions hierarchically for occlusion-free and controllable visual exploration. We enable observation and exploration of the relationships among field line clusters, spatiotemporal regions and their interconnection in the transformed space. Most viewpoint selection methods only consider the external viewpoints outside of the flow field. This will not convey a clear observation when the flow field is clutter on the boundary side. Therefore, we propose a new way to explore flow fields by selecting several internal viewpoints around the flow features inside of the flow field and then generating a B-Spline curve path traversing these viewpoints to provide users with closeup views of the flow field for detailed observation of hidden or occluded internal flow features [54]. This work is also extended to deal with unsteady flow fields. Besides flow field visualization, some other topics relevant to visualization also attract my attention. In iGraph [31], we leverage a distributed system along with a tiled display wall to provide users with high-resolution visual analytics of big image and text collections in real time. Developing pedagogical visualization tools forms my other research focus. Since most cryptography algorithms use sophisticated mathematics, it is difficult for beginners to understand both what the algorithm does and how the algorithm does that. Therefore, we develop a set of visualization tools to provide users with an intuitive way to learn and understand these algorithms.

The method of characteristics in compressible flow.

Reproduced from typewritten copy.

Hypersonic dynamic stability.

pt. 1. Summary [by] L. A. Marshall.--pt. 2. Conical body experimental program [by] R. B. Hobbs, Jr.--pt. 3. Unsteady flow field program [by] H. Rie, E. A. Linkiewicz [and] F. D. Bosworth.

Scaling of natural convection of an inclined flat plate : ramp cooling condition

A scaling analysis is performed for the transient boundary layer established adjacent to an inclined flat plate following a ramp cooling boundary condition. The imposed wall temperature decreases linearly up to a specific value over a specific time. It is revealed that if the ramp time is sufficiently large then the boundary layer reaches quasi-steady mode before the growth of the temperature is finished. However, if the ramp time is shorter then the steady state of the boundary layer may be reached after the growth of the temperature is completed. In this case, the ultimate steady state is the same as if the start up had been instantaneous. Note that the cold boundary layer adjacent to the plate is potentially unstable to Rayleigh-Bénard instability if the Rayleigh number exceeds a certain critical value for this cooling case. The onset of instability may set in at different stages of the boundary layer development. A proper identification of the time when the instability may set in is discussed. A numerical verification of the time for the onset of instability is presented in this study. Different flow regimes based on the stability of the boundary layer have also been discussed with numerical results.

Natural convection boundary-layer adjacent to an inclined flat plate subject to sudden and ramp heating

The natural convection thermal boundary layer adjacent to an inclined flat plate subject to sudden heating and a temperature boundary condition which follows a ramp function up until a specified time and then remains constant is investigated. The development of the flow from start-up to a steady-state has been described based on scaling analyses and verified by numerical simulations. Different flow regimes based on the Rayleigh number are discussed with numerical results for both boundary conditions. For ramp heating, the boundary layer flow depends on the comparison of the time at which the ramp heating is completed and the time at which the boundary layer completes its growth. If the ramp time is long compared with the steady state time, the layer reaches a quasi steady mode in which the growth of the layer is governed solely by the thermal balance between convection and conduction. On the other hand, if the ramp is completed before the layer becomes steady; the subsequent growth is governed by the balance between buoyancy and inertia, as for the case of instantaneous heating.

Scaling of natural convection of an inclined flat plate : sudden cooling condition

The natural convection boundary layer adjacent to an inclined plate subject to sudden cooling boundary condition has been studied. It is found that the cold boundary layer adjacent to the plate is potentially unstable to Rayleigh-Bénard instability if the Rayleigh number exceeds a certain critical value. A scaling relation for the onset of instability of the boundary layer is achieved. The scaling relations have been developed by equating important terms of the governing equations based on the development of the boundary layer with time. The flow adjacent to the plate can be classified broadly into a conductive, a stable convective or an unstable convective regime determined by the Rayleigh number. Proper scales have been established to quantify the flow properties in each of these flow regimes. An appropriate identification of the time when the instability may set in is discussed. A numerical verification of the time for the onset of instability is also presented in this study. Different flow regimes based on the stability of the boundary layer have been discussed with numerical results.

Natural convection in attic-shaped spaces subject to sudden and ramp heating boundary conditions

In this study, a discussion of the fluid dynamics in the attic space is reported, focusing on its transient response to sudden and linear changes of temperature along the two inclined walls. The transient behaviour of an attic space is relevant to our daily life. The instantaneous and non-instantaneous (ramp) heating boundary condition is applied on the sloping walls of the attic space. A theoretical understanding of the transient behaviour of the flow in the enclosure is performed through scaling analysis. A proper identification of the timescales, the velocity and the thickness relevant to the flow that develops inside the cavity makes it possible to predict theoretically the basic flow features that will survive once the thermal flow in the enclosure reaches a steady state. A time scale for the heating-up of the whole cavity together with the heat transfer scales through the inclined walls has also been obtained through scaling analysis. All scales are verified by the numerical simulations.

Natural convection in attics subject to instantaneous and ramp cooling boundary conditions

A fundamental study of the fluid dynamics inside an attic shaped triangular enclosure with cold upper walls and adiabatic horizontal bottom wall is reported in this study. The transient behaviour of the attic fluid which is relevant to our daily life is examined based on a scaling analysis. The transient phenomenon begins with the instantaneous cooling and the cooling with linear decreases of temperature up to some specific time (ramp time) and then maintain constant of the upper sloped walls. It is shown that both inclined walls develop a thermal boundary layer whose thicknesses increase towards steady-state or quasi-steady values. A proper identification of the timescales, the velocity and the thickness relevant to the flow that develops inside the cavity makes it possible to predict theoretically the basic flow features that will survive once the thermal flow in the enclosure reaches a steady state. A time scale for the cooling-down of the whole cavity together with the heat transfer scales through the inclined walls has also been obtained through scaling analysis. All scales are verified by the numerical simulations.

Scaling of free convection heat transfer in a triangular cavity for Pr>1

Unsteady natural convection inside a triangular cavity subject to a non-instantaneous heating on the inclined walls in the form of an imposed temperature which increases linearly up to a prescribed steady value over a prescribed time is reported. The development of the flow from start-up to a steady-state has been described based on scaling analyses and direct numerical simulations. The ramp temperature has been chosen in such a way that the boundary layer is reached a quasi-steady mode before the growth of the temperature is completed. In this mode the thermal boundary layer at first grows in thickness, then contracts with increasing time. However, if the imposed wall temperature growth period is sufficiently short, the boundary layer develops differently. It is seen that the shape of many houses are isosceles triangular cross-section. The heat transfer process through the roof of the attic-shaped space should be well understood. Because, in the building energy, one of the most important objectives for design and construction of houses is to provide thermal comfort for occupants. Moreover, in the present energy-conscious society it is also a requirement for houses to be energy efficient, i.e. the energy consumption for heating or air-conditioning houses must be minimized.