Special Issue “ Fluid Flow, Energy Transfer and Design”

Fluids are important because of their preponderance in our lives. Fluid mechanics touches almost every aspect of our daily lives, and it plays a central role in many branches of science and technology. Therefore, it is a challenging and exciting field of scientific activity due to the complexity of the subject studied and the breadth of the applications. The quest for advances in fluid mechanics, as in other scientific fields, emerge from analytical, computational (CFD) and experimental studies. The improvement in our ability to describe, predict and control the phenomena played (and plays) key roles in the technological breakthroughs. The present theme issue of “Fluid and Heat Flow: Simulation and Optimization” collects a selection of papers. selection of papers presented at Special Session “Fluid Flow, Energy Transfer and Design”

Prediction of fluidization behaviour and a quasi-stationary approach to drying kinetics of irregular particulate food materials

Changes in fluidization behaviour behaviour was characterised for parallelepiped particles with three aspect ratios, 1:1, 2:1 and 3:1 and spherical particles. All drying experiments were conducted at 500C and 15 % RH using a heat pump dehumidifier system. Fluidization experiments were undertaken for the bed heights of 100, 80, 60 and 40 mm and at 10 moisture content levels. Due to irregularities in shape minimum fluidisation velocity of parallelepiped particulates (potato) could not fitted to any empirical model. Also a generalized equation was used to predict minimum fluidization velocity. The modified quasi-stationary method (MQSM) has been proposed to describe drying kinetics of parallelepiped particulates at 30o C, 40o C and 50o C that dry mostly in the falling rate period in a batch type fluid bed dryer.

Velocity and Concentration field of discharge from a boat propeller

Two-stroke outboard boat engines using total loss lubrication deposit a significant proportion of their lubricant and fuel directly into the water. The purpose of this work is to document the velocity and concentration field characteristics of a submerged swirling water jet emanating from a propeller in order to provide information on its fundamental characteristics. Measurements of the velocity and concentration field were performed in a turbulent jet generated by a model boat propeller (0.02 m diameter) operating at 1500 rpm and 3000 rpm. The measurements were carried out in the Zone of Established Flow up to 50 propeller diameters downstream of the propeller. Both the mean axial velocity profile and the mean concentration profile showed self-similarity. Further, the stand deviation growth curve was linear. The effects of propeller speed and dye release location were also investigated.

Computational fluid dynamics for improved bioreactor design and 3D culture

The complex relationship between the hydrodynamic environment and surrounding tissues directly impacts on the design and production of clinically useful grafts and implants. Tissue engineers have generally seen bioreactors as 'black boxes' within which tissue engineering constructs (TECs) are cultured. It is accepted that a more detailed description of fluid mechanics and nutrient transport within process equipment can be achieved by using computational fluid dynamics (CFD) technology. This review discusses applications of CFD for tissue engineering-related bioreactors -- fluid flow processes have direct implications on cellular responses such as attachment, migration and proliferation. We conclude that CFD should be seen as an invaluable tool for analyzing and visualizing the impact of fluidic forces and stresses on cells and TECs.

Bioreactor Systems for Tissue Engineering

The hydrodynamic environment “created” by bioreactors for the culture of a tissue engineered construct (TEC) is known to influence cell migration, proliferation and extra cellular matrix production. However, tissue engineers have looked at bioreactors as black boxes within which TECs are cultured mainly by trial and error, as the complex relationship between the hydrodynamic environment and tissue properties remains elusive, yet is critical to the production of clinically useful tissues. It is well known in the chemical and biotechnology field that a more detailed description of fluid mechanics and nutrient transport within process equipment can be achieved via the use of computational fluid dynamics (CFD) technology. Hence, the coupling of experimental methods and computational simulations forms a synergistic relationship that can potentially yield greater and yet, more cohesive data sets for bioreactor studies. This review aims at discussing the rationale of using CFD in bioreactor studies related to tissue engineering, as fluid flow processes and phenomena have direct implications on cellular response such as migration and/or proliferation. We conclude that CFD should be seen by tissue engineers as an invaluable tool allowing us to analyze and visualize the impact of fluidic forces and stresses on cells and TECs.

Investigating the hydrodynamic performance of a gross pollutant trap

Field studies show that the internal screens in a gross pollutant trap (GPT) are often clogged with organic matter, due to infrequent cleaning. The hydrodynamic performance of a GPT with fully blocked screens was comprehensively investigated under a typical range of onsite operating conditions. Using an acoustic Doppler velocimeter (ADV), velocity profiles across three critical sections of the GPT were measured and integrated to examine the net fluid flow at each section. The data revealed that when the screens are fully blocked, the flow structure within the GPT radically changes. Consequently, the capture/retention performance of the device rapidly deteriorates. Good agreement was achieved between the experimental and the previous 2D computational fluid dynamics (CFD) velocity profiles for the lower GPT inlet flow conditions.

Exit humidity of wet scrubbers for underground coal mines

A wet scrubber is a device used in underground coal mines for the exhaust treatment system of various internal combustion engines (generally diesel) primarily as a spark arrestor with a secondary function to remove pollutants from the exhaust gas. A pool of scrubbing liquid (generally water based) is used in conjunction with a Diesel Particulate Filter (DPF). Scrubbers are widely used in underground applications of diesel engines as their exhaust contains high concentration of harmful diesel particulate matter (DPM) and other pollutant gases. Currently the DPFs have to be replaced frequently because moisture output from the wet scrubber blocks the filter media and causes reduced capacity. This paper presents experimental and theoretical studies on the heat and mass transfer mechanisms of the exhaust flow both under and above the water surface, aiming at finding the cause and effects of the moisture reaching the filters and employing a solution to reduce the humidity and DPM output, and to prolong the change-out period of the DPF. By assuming a steady flow condition, heat transfer from the inlet exhaust gas balances energy required for the water evaporation. Hence the exit humidity will decrease with the increase of exit temperature. Experiments on a real scrubber are underway.

Hydraulic instability onset detection in Kaplan turbines by monitoring shaft vibrations

Design of hydraulic turbines has often to deal with hydraulic instability. It is well-known that Francis and Kaplan types present hydraulic instability in their design power range. Even if modern CFD tools may help to define these dangerous operating conditions and optimize runner design, hydraulic instabilities may fortuitously arise during the turbine life and should be timely detected in order to assure a long-lasting operating life. In a previous paper, the authors have considered the phenomenon of helical vortex rope, which happens at low flow rates when a swirling flow, in the draft tube conical inlet, occupies a large portion of the inlet. In this condition, a strong helical vortex rope appears. The vortex rope causes mechanical effects on the runner, on the whole turbine and on the draft tube, which may eventually produce severe damages on the turbine unit and whose most evident symptoms are vibrations. The authors have already shown that vibration analysis is suitable for detecting vortex rope onset, thanks to an experimental test campaign performed during the commissioning of a 23 MW Kaplan hydraulic turbine unit. In this paper, the authors propose a sophisticated data driven approach to detect vortex rope onset at different power load, based on the analysis of the vibration signals in the order domain and introducing the so-called "residual order spectrogram", i.e. an order-rotation representation of the vibration signal. Some experimental test runs are presented and the possibility to detect instability onset, especially in real-time, is discussed.

Influence of tilt and surface roughness on the outflow wind field of an impinging jet

A physical and numerical steady flow impinging jet has been used to simulate the bulk characteristics of a downburst-like wind field. The influence of downdraft tilt and surface roughness on the ensuing wall jet flow has been investigated. It was found that a simulated downdraft impinging the surface at a non-normal angle has the potential for causing larger structural loads than the normal impingement case. It was also found that for the current impinging jet simulations, surface roughness played a minor role in determining the storm maximum wind structure, but this influence increased as the wall jet diverged. However, through comparison with previous research it was found that the influence of surface roughness is Reynolds number dependent and therefore may differ from that reported herein for full-scale downburst cases. Using the current experimental results an empirical model has been developed for laboratory-scale impinging jet velocity structure that includes the influence of both jet tilt and surface roughness.

Combustion enhancement in a scramjet engine using oxygen enrichment and porous fuel injection

This paper reports on the experimental testing of oxygen-enriched porous fuel injection in a scramjet engine. Fuel was injected via inlet mounted, oxide-based ceramic matrix composite (CMC) injectors on both flow path surfaces that covered a total of 9.2 % of the intake surface area. All experiments were performed at an enthalpy of 3.93−4.25±3.2% MJ kg−1, flight Mach number 9.2–9.6 and an equivalence ratio of 0.493±3%. At this condition, the engine was shown to be on the verge of achieving appreciable combustion. Oxygen was then added to the fuel prior to injection such that two distinct enrichment levels were achieved. Combustion was found to increase, by as much as 40 % in terms of combustion-induced pressure rise, over the fuel-only case with increasing oxygen enrichment. Further, the onset of combustion was found to move upstream with increasing levels of oxygen enrichment. Thrust, both uninstalled and specific, and specific impulse were found to be improved with oxygen enrichment. Enhanced fuel–air mixing due to the pre-mixing of oxygen with the fuel together with the porous fuel injection are believed to be the main contributors to the observed enhanced performance of the tested engine.

Impaction of spray droplets on leaves: Influence of formulation and leaf character on shatter, bounce and adhesion

This paper combines experimental data with simple mathematical models to investigate the influence of spray formulation type and leaf character (wettability) on shatter, bounce and adhesion of droplets impacting with cotton, rice and wheat leaves. Impaction criteria that allow for different angles of the leaf surface and the droplet impact trajectory are presented; their predictions are based on whether combinations of droplet size and velocity lie above or below bounce and shatter boundaries. In the experimental component, real leaves are used, with all their inherent natural variability. Further, commercial agricultural spray nozzles are employed, resulting in a range of droplet characteristics. Given this natural variability, there is broad agreement between the data and predictions. As predicted, the shatter of droplets was found to increase as droplet size and velocity increased, and the surface became harder to wet. Bouncing of droplets occurred most frequently on hard to wet surfaces with high surface tension mixtures. On the other hand, a number of small droplets with low impact velocity were observed to bounce when predicted to lie well within the adhering regime. We believe this discrepancy between the predictions and experimental data could be due to air layer effects that were not taken into account in the current bounce equations. Other discrepancies between experiment and theory are thought to be due to the current assumption of a dry impact surface, whereas, in practice, the leaf surfaces became increasingly covered with fluid throughout the spray test runs.

A nonconservative Lagrangian framework for statistical fluid registration - SAFIRA

In this paper, we used a nonconservative Lagrangian mechanics approach to formulate a new statistical algorithm for fluid registration of 3-D brain images. This algorithm is named SAFIRA, acronym for statistically-assisted fluid image registration algorithm. A nonstatistical version of this algorithm was implemented, where the deformation was regularized by penalizing deviations from a zero rate of strain. In, the terms regularizing the deformation included the covariance of the deformation matrices Σ and the vector fields (q). Here, we used a Lagrangian framework to reformulate this algorithm, showing that the regularizing terms essentially allow nonconservative work to occur during the flow. Given 3-D brain images from a group of subjects, vector fields and their corresponding deformation matrices are computed in a first round of registrations using the nonstatistical implementation. Covariance matrices for both the deformation matrices and the vector fields are then obtained and incorporated (separately or jointly) in the nonconservative terms, creating four versions of SAFIRA. We evaluated and compared our algorithms' performance on 92 3-D brain scans from healthy monozygotic and dizygotic twins; 2-D validations are also shown for corpus callosum shapes delineated at midline in the same subjects. After preliminary tests to demonstrate each method, we compared their detection power using tensor-based morphometry (TBM), a technique to analyze local volumetric differences in brain structure. We compared the accuracy of each algorithm variant using various statistical metrics derived from the images and deformation fields. All these tests were also run with a traditional fluid method, which has been quite widely used in TBM studies. The versions incorporating vector-based empirical statistics on brain variation were consistently more accurate than their counterparts, when used for automated volumetric quantification in new brain images. This suggests the advantages of this approach for large-scale neuroimaging studies.

Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter = 9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 10(8). Under these conditions the convection is turbulent, and the time-averaged velocity at any point is `zero'. The Reynolds number based on the Taylor microscale, Re-lambda, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6-7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as (RaSc1/2)-Sc-1/2, and the Reynolds number would scale as (RaSc-1/2)-Sc-1/2. The velocity and the flux measurements appear to be consistent with the Ra-1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are Much higher compared to what would be obtained in Rayleigh-Benard (R-B) convection for similar density differences.

On the origin of the inflectional instability of a laminar separation bubble

This is an experimental and theoretical Study of a laminar separation bubble and the associated linear stability mechanisms. Experiments were performed over a flat plate kept in a wind tunnel, with an imposed pressure gradient typical of an aerofoil that would involve a laminar separation bubble. The separation bubble was characterized by measurement of surface-pressure distribution and streamwise velocity using hot-wire anemometry. Single component hot-wire anemometry was also used for a detailed study of the transition dynamics. It was foundthat the so-called dead-air region in the front portion of the bubble corresponded to a region of small disturbance amplitudes, with the amplitude reaching a maximum value close to the reattachment point. An exponential growth rate of the disturbance was seen in the region upstream of the mean maximum height of the bubble, and this was indicative of a linear instability mechanism at work. An infinitesimal disturbance was impulsively introduced into the boundary layer upstream of separation location, and the wave packet was tracked (in an ensemble-averaged sense) while it was getting advected downstream. The disturbance was found to be convective in nature. Linear stability analyses (both the Orr-Sommerfeld and Rayleigh calculations) were performed for mean velocity profiles, starting from an attached adverse-pressure-gradient boundary layer all the way up to the front portion of the separation-bubble region (i.e. up to the end of the dead-air region in which linear evolution of the disturbance could be expected). The conclusion from the present work is that the primary instability mechanism in a separation bubble is inflectional in nature, and its origin can be traced back to upstream of the separation location. In other words, the inviscid inflectional instability of the separated shear layer should be logically seen as an extension of the instability of the upstream attached adverse-pressure-gradient boundary layer. This modifies the traditional view that pegs the origin of the instability in a separation bubble to the detached shear layer Outside the bubble, with its associated Kelvin-Helmholtz mechanism. We contendthat only when the separated shear layer has moved considerably away from the wall (and this happens near the maximum-height location of the mean bubble), a description by the Kelvin-Helmholtz instability paradigm, with its associated scaling principles, Could become relevant. We also propose a new scaling for the most amplified frequency for a wall-bounded shear layer in terms of the inflection-point height and the vorticity thickness and show it to be universal.