Automatic Compilation of MATLAB Programs for Synergistic Execution on Heterogeneous Processors

MATLAB is an array language, initially popular for rapid prototyping, but is now being increasingly used to develop production code for numerical and scientific applications. Typical MATLAB programs have abundant data parallelism. These programs also have control flow dominated scalar regions that have an impact on the program's execution time. Today's computer systems have tremendous computing power in the form of traditional CPU cores and throughput oriented accelerators such as graphics processing units(GPUs). Thus, an approach that maps the control flow dominated regions to the CPU and the data parallel regions to the GPU can significantly improve program performance. In this paper, we present the design and implementation of MEGHA, a compiler that automatically compiles MATLAB programs to enable synergistic execution on heterogeneous processors. Our solution is fully automated and does not require programmer input for identifying data parallel regions. We propose a set of compiler optimizations tailored for MATLAB. Our compiler identifies data parallel regions of the program and composes them into kernels. The problem of combining statements into kernels is formulated as a constrained graph clustering problem. Heuristics are presented to map identified kernels to either the CPU or GPU so that kernel execution on the CPU and the GPU happens synergistically and the amount of data transfer needed is minimized. In order to ensure required data movement for dependencies across basic blocks, we propose a data flow analysis and edge splitting strategy. Thus our compiler automatically handles composition of kernels, mapping of kernels to CPU and GPU, scheduling and insertion of required data transfer. The proposed compiler was implemented and experimental evaluation using a set of MATLAB benchmarks shows that our approach achieves a geometric mean speedup of 19.8X for data parallel benchmarks over native execution of MATLAB.

Dispersion and vaporization of biofuels and conventional fuels in a crossflow pre-mixer

In lean premixed pre-vaporized (LPP) combustion, controlled atomization, dispersion and vaporization of different types of liquid fuel in the premixer are the key factors required to stabilize the combustion process and improve the efficiency. The dispersion and vaporization process for biofuels and conventional fuels sprayed into a crossflow pre-mixer have been simulated and analyzed with respect to vaporization rate, degree of mixedness and homogeneity. Two major biofuels under investigation are Ethanol and Rapeseed Methyl Esters (RME), while conventional fuels are gasoline and jet-A. First, the numerical code is validated by comparing with the experimental data of single n-heptane and decane droplet evaporating under both moderate and high temperature convective air now. Next, the spray simulations were conducted with monodispersed droplets with an initial diameter of 80 mu m injected into a turbulent crossflow of air with a typical velocity of 10 m/s and temperature of around 800K. Vaporization time scales of different fuels are found to be very different. The droplet diameter reduction and surface temperature rise were found to be strongly dependent on the fuel properties. Gasoline droplet exhibited a much faster vaporization due a combination of higher vapor pressure and smaller latent heat of vaporization compared to other fuels. Mono-dispersed spray was adopted with the expectation of achieving more homogeneous fuel droplet size than poly-dispersed spray. However, the diameter histogram in the zone near the pre-mixer exit shows a large range of droplet diameter distributions for all the fuels. In order to improve the vaporization performance, fuels were pre-heated before injection. Results show that the Sauter mean diameter of ethanol improved from 52.8% of the initial injection size to 48.2%, while jet-A improved from 48.4% to 18.6% and RME improved from 63.5% to 31.3%. The diameter histogram showed improved vaporization performance of jet-A. (C) 2011 Elsevier Ltd. All rights reserved.

Experimental Studies on High-Pressure Spray Structure of Biofuels

The high-pressure spray characteristics of biofuels, specifically, Pongamia oil and its blends with diesel are studied for various gas pressures. Two single-hole solenoid injectors with nozzle diameters of 200 and 260 mu m are used along with a high-pressure common-rail direct-injection system to inject fuel into a high-pressure spray visualization chamber. The spray structure is characterized using a high-speed laser-based shadowgraphy technique. The spray structure of Pongamia oil revealed the presence of an intact liquid core at low gas pressure. At high gas pressures, the spray atomization of the Pongamia oil showed marked improvement. The spray tip penetration of Pongamia oil and its blends with diesel is higher compared to that of diesel for all test conditions. The spray cone angle of Pongamia oil and 50% Pongamia oil blend with diesel is lower as compared to that of diesel. Both these observations are attributed to the presence of large droplets carrying higher momentum in oil and blend. The droplet size is measured at an injection pressure of 1000 bar and gas pressure of 30 bar at 25 mm below the nozzle tip using the particle/droplet image.analysis (PDIA) method. The droplet size measurements have shown that the Sauter mean diameter (SMD) in the spray core of Pongamia oil is more than twice that of diesel. The spray tip penetration of the 20% blend of Pongamia with diesel (P20) is similar to that of diesel but the SMD is 50% higher. Based on experimental data, appropriate spray tip penetration correlation is proposed for the vegetable oil fuels such as Pongamia.

Spray characterization of gasoline-ethanol blends from a multi-hole port fuel injector

This work reports the measured spray structure and droplet size distributions of ethanol-gasoline blends for a low-pressure, multi-hole, port fuel injector (PFI). This study presents previously unavailable data for this class of injectors which are widely used in automotive applications. Specifically, gasoline, ethanol, and gasoline-ethanol blends containing 10%, 20% and 50% ethanol were studied using laser backlight imaging, and particle/droplet image analysis (PDIA) techniques. The fuel mass injected, spray structure and tip penetrations, droplet size distributions, and Sauter mean diameter were determined for the blends, at two different injection pressures. Results indicate that the gasoline and ethanol sprays have similar characteristics in terms of spray progression and droplet sizes in spite of the large difference in viscosity. It appears that the complex mode of atomization utilized in these injectors involving interaction of multiple fuel jets is fairly insensitive to the fuel viscosity over a range of values. This result has interesting ramifications for existing gasoline fuel systems which need to handle blends and even pure ethanol, which is one of the renewable fuels of the future. (C) 2012 Elsevier Ltd. All rights reserved.

Bubble size measurement and error analysis in a gas liquid ejector

Bubble size in a gas liquid ejector has been measured using the image technique and analysed for estimation of Sauter mean diameter. The individual bubble diameter is estimated by considering the two dimensional contour of the ellipse, for the actual three dimensional ellipsoid in the system by equating the volume of the ellipsoid to that of the sphere. It is observed that the bubbles are of oblate and prolate shaped ellipsoid in this air water system. The bubble diameter is calculated based on this concept and the Sauter mean diameter is estimated. The error between these considerations is reported. The bubble size at different locations from the nozzle of the ejector is presented along with their percentage error which is around 18%.

AN EXPERIMENTAL STUDY ON EFFERVESCENT ATOMIZATION OF BIO-OIL FUELS

This paper presents the work on detailed characterization of effervescent spray of Jatropha and Pongamia pure plant oils. The spray characteristics of these biofuels are compared with those of diesel. Both macroscopic and microscopic spray characteristics at different injection pressures and gas-to-liquid ratio (GLR) have been studied. The particle/droplet imaging analysis (PDIA) technique along with direct imaging methods are used for the purpose of spray characterization. Due to their higher viscosity, pure plant oils showed poor atomization compared to diesel and a blend of diesel and pure plant oil at a given GLR. Pure plant oil sprays showed a lower spray cone angle when compared to diesel and blends at lower GLRs. However, the difference is not significant at higher GLRs. Droplet size measurements at 100 mm downstream of the exit orifice showed reduction in Sauter mean diameter (SMD) diameter with increase in GLR. A radial variation in the SMD is observed for the blend and pure plant oils. Pure oils showed a larger variation when compared to the blend. Spray unsteadiness has been characterized based on the image-to-image variation in the mean droplet diameter and fluctuations in the spray cone angle. Results showed that pure plant oil sprays are more unsteady at lower GLRs when compared to diesel and blend. A critical GLR is identified at which the spray becomes steady. The three regimes of spray operation, namely ``steady spray,'' ``pulsating spray,'' and ``spray and unbroken liquid jet'' are identified in the injection pressure-GLR parameter space for these pure plant oils. Two-phase flow imaging inside the exit orifice shows that for the pure plant oils, the flow is highly transient at low GLRs and the bubbly, slug, and annular two-phase flow regimes are all observed. However, at higher GLRs where the spray is steady, only the annular flow regime is observed.

Automatic compilation of MATLAB programs for synergistic execution on heterogeneous processors

MATLAB is an array language, initially popular for rapid prototyping, but is now being increasingly used to develop production code for numerical and scientific applications. Typical MATLAB programs have abundant data parallelism. These programs also have control flow dominated scalar regions that have an impact on the program's execution time. Today's computer systems have tremendous computing power in the form of traditional CPU cores and throughput oriented accelerators such as graphics processing units(GPUs). Thus, an approach that maps the control flow dominated regions to the CPU and the data parallel regions to the GPU can significantly improve program performance. In this paper, we present the design and implementation of MEGHA, a compiler that automatically compiles MATLAB programs to enable synergistic execution on heterogeneous processors. Our solution is fully automated and does not require programmer input for identifying data parallel regions. We propose a set of compiler optimizations tailored for MATLAB. Our compiler identifies data parallel regions of the program and composes them into kernels. The problem of combining statements into kernels is formulated as a constrained graph clustering problem. Heuristics are presented to map identified kernels to either the CPU or GPU so that kernel execution on the CPU and the GPU happens synergistically and the amount of data transfer needed is minimized. In order to ensure required data movement for dependencies across basic blocks, we propose a data flow analysis and edge splitting strategy. Thus our compiler automatically handles composition of kernels, mapping of kernels to CPU and GPU, scheduling and insertion of required data transfer. The proposed compiler was implemented and experimental evaluation using a set of MATLAB benchmarks shows that our approach achieves a geometric mean speedup of 19.8X for data parallel benchmarks over native execution of MATLAB.

A method for measurement of planar liquid volume fraction in dense sprays

A methodology for measurement of planar liquid volume fraction in dense sprays using a combination of Planar Laser-Induced Fluorescence (PLIF) and Particle/Droplet Imaging Analysis (PDIA) is presented in this work. The PLIF images are corrected for loss of signal intensity due to laser sheet scattering, absorption and auto-absorption. The key aspect of this work pertains to simultaneously solving the equations involving the corrected PLIF signal and liquid volume fraction. From this, a quantitative estimate of the planar liquid volume fraction is obtained. The corrected PLIF signal and the corrected planar Mie scattering can be also used together to obtain the Sauter Mean Diameter (SMD) distribution by using data from the PDIA technique at a particular location for calibration. This methodology is applied to non-evaporating sprays of diesel and a more viscous pure plant oil at an injection pressure of 1000 bar and a gas pressure of 30 bar in a high pressure chamber. These two fuels are selected since their viscosity values are very different with a consequently very different spray structure. The spatial distribution of liquid volume fraction and SMD is obtained for two fuels. The proposed method is validated by comparing liquid volume fraction obtained by the current method with data from PDIA technique. (C) 2012 Elsevier Inc. All rights reserved.

Experimental studies on air-assisted impinging jet atomization

In the present study, a novel air-assisted impinging jet atomization is demonstrated. A configuration in which a gas jet is directed on to the impinging point of two liquid jets is used to improve the atomization. The effect of liquid properties such as viscosity and surface tension, angle between liquid jets and gas injection orifice diameter on spray characteristics has been experimentally studied. Backlit imaging and particle/droplet imaging and analysis techniques are utilized to characterize the sprays. The experimental results indicate that the effect of liquid viscosity is significant on the liquid sheet break up formed by the impinging jets. However, surface tension does not affect the spray structure significantly in this mode of atomization. At low liquid jet velocity, the prompt mode of atomization is observed where as atomization occurs in classical mode at higher liquid jet velocity. Results showed that variation in the angle between liquid jets do not affect the breakup phenomenon significantly. The spray angle is computed by finding the angle between the lines joining the impinging point and spray edge at an axial distance of 15 mm downstream of the impinging point from the ensemble-averaged data over 100 spray images. It was observed that effect of liquid jets impinging angle on the spray angle is higher at higher liquid velocity. Higher viscosity liquids exhibit lower spray angles. Droplet size measurements indicate a radial variation in the spray. An overall Sauter Mean Diameter (SMD) value is obtained by combining the droplet statistics at all radial locations at a fixed axial location. A very interesting trend is that the SMD is constant beyond a critical Gas to Liquid Ratio (GLR) and momentum ratio for a large variation in liquid viscosity and surface tension. This observation has important ramifications for fuel flexible systems. (C) 2013 Elsevier Ltd. All rights reserved.

Toward a Scalable Working Set Size Estimation Method and Its Application for Chip Multiprocessors

It is essential to accurately estimate the working set size (WSS) of an application for various optimizations such as to partition cache among virtual machines or reduce leakage power dissipated in an over-allocated cache by switching it OFF. However, the state-of-the-art heuristics such as average memory access latency (AMAL) or cache miss ratio (CMR) are poorly correlated to the WSS of an application due to 1) over-sized caches and 2) their dispersed nature. Past studies focus on estimating WSS of an application executing on a uniprocessor platform. Estimating the same for a chip multiprocessor (CMP) with a large dispersed cache is challenging due to the presence of concurrently executing threads/processes. Hence, we propose a scalable, highly accurate method to estimate WSS of an application. We call this method ``tagged WSS (TWSS)'' estimation method. We demonstrate the use of TWSS to switch-OFF the over-allocated cache ways in Static and Dynamic NonUniform Cache Architectures (SNUCA, DNUCA) on a tiled CMP. In our implementation of adaptable way SNUCA and DNUCA caches, decision of altering associativity is taken by each L2 controller. Hence, this approach scales better with the number of cores present on a CMP. It gives overall (geometric mean) 26% and 19% higher energy-delay product savings compared to AMAL and CMR heuristics on SNUCA, respectively.

Studies on microscopic structure of diesel sprays under atmospheric and high gas pressures

In the present work, the spray structure of diesel from a 200-mu m, single-hole solenoid injector is studied using microscopic imaging at injection pressures of 700, 1000 and 1400 bar for various gas pressures. A long-distance microscope with a high resolution camera is used for spray visualization with a direct imaging technique. This study shows that even at very high injection pressures, the spray structure in an ambient environment of atmospheric pressure reveals presence of entangled ligaments and non-spherical droplets during the injection period. With increase in the injection pressure, the ligaments tend to get smaller and spread radially. The spray structure studies are also conducted at high gas pressures in a specially designed high pressure chamber with optical access. The near nozzle spray structure at the end of the injection shows that the liquid jet breakup is improved with increase in gas density. The droplet size measurement is possible only late in the injection duration when the breakup appears to be complete and mostly spherical droplets are observed. Hence, droplet size measurements are performed after 1.3 ms from start of the injection pulse. Spatial and temporal variation in Sauter Mean `Diameter (SMD) is observed and reported for the case corresponding to an injection pressure of 700 bar. Overall, this study has highlighted the importance of verifying the extentof atomization and droplet shape even in dense sprays before using conventional dropsizing methods such as PDPA.

Correlation between Optical Properties and Nanomorphology of Fluoranthene-Based Conjugated Copolymer

Nanoparticles of conjugated polymers are receiving attention due to their interesting optical properties. Here we report nanoparticles of fluoranthene-based conjugated copolymer prepared by the Suzuki coupling reaction. The copolymer forms nanoparticles by the spontaneous self-assembly after evaporation of organic solvent. The mean diameter of the nanoparticles can be manipulated by varying solvent composition. We investigated the parameters that govern the nanostructured morphology of polymer by systematic variation of good and poor solvent. The UV vis and time-resolved fluorescence spectroscopy measurement reveal the use of poor solvent in the organization of nanostructures. Furthermore, transmission electron microscopy highlights the importance of rigidity of the polymer backbone in morphological development.

Airblast spray in crossflow - Structure, trajectory and droplet sizing

This study reports results of an experimental investigation of airblast spray of water and ethanol in crossflow. Laser shadowgraphy and Particle/Droplet Imaging Analysis (PDIA) are used to derive spray trajectory and drop size information while Particle Tracking Velocimetry (PTV) is used to measure droplet velocities. A new phenomenon of spray bifurcation is observed for low Gas to Liquid Ratio (GLR) cases. The reasons for the spatial bifurcation can be attributed to a combination of reasons. These are (a) presence of large ligaments and droplets in the near-nozzle region for low GLRs (b) secondary breakup experienced by ligaments/droplets leading to formation of a large number of small droplets, and (c) the crossflow causing differential dispersion of the small and large droplets. A novel correlation for spray trajectory is proposed incorporating the momentum ratio and liquid surface tension. This correlation is shown to be effective in predicting the non-linear spray trajectory over a large range of conditions for not only water but ethanol and Jet-A also. It is observed that the larger droplets penetrate further into the crossflow, in the direction of injection. Thus, with increase in height of the measurement location from the injection plane, the droplet Sauter Mean Diameter (SMD) is found to increase. Moreover, as the droplets travel downstream in the crossflow direction, the droplet SMD is observed to decrease. The effect of drag is assessed by comparing velocity of different sizes of droplets at various locations. Smaller droplets are entrained into the crossflow at much lower elevations, whereas larger droplets tend to penetrate further into the crossflow. (C) 2015 Elsevier Ltd. All rights reserved.

Analysis on the atomization characteristics of aviation biofuel discharging from simplex swirl atomizer

The atomization characteristics of aviation biofuel discharging from a simplex swirl atomizer into quiescent atmospheric air are studied. The aviation biofuel is a mixture of 90% commercially available camelina-derived biofuel and 10% VonSol-53 (aromatics). The experiments are conducted in a spray test facility at varying fuel flow rate conditions. The measured characteristics include atomizer flow number, spray cone angle, breakup length of liquid sheet, wavelength of undulations on liquid sheet, and spray droplet size. The characteristics of biofuel sheet breakup are deduced from the captured images of biofuel spray. The measurements of spray droplet size distribution are obtained using Spraytec. The experimentally measured characteristics of the biofuel sheet breakup are compared with the predictions obtained from the liquid film breakup model proposed by Senecal et al. (1999). The measurements of wavelength and breakup length of the biofuel sheet discharging from the simplex swirl atomizer agree well with the model predictions. The model-predicted droplet size for the biofuel spray is significantly higher than the experimentally measured Sauter mean diameter (SMD). The spray droplets formed from the liquid sheet breakup undergo secondary atomization until 35-45 mm from the atomizer exit and thereafter the SMD increases downstream due to the combined effect of fuel evaporation and droplet coalescence. A good comparison is observed between the experimentally measured SMD of the biofuel spray and the predictions obtained using the empirical correlation reported in literature for sprays discharging from simplex swirl atomizers. (C) 2015 Elsevier Ltd. All rights reserved.

EXPERIMENTAL STUDIES ON AIR-ASSISTED ATOMIZATION OF JATROPHA PURE PLANT OIL

The present study focuses on exploring air-assisted atomization strategies for effective atomization of high-viscosity biofuels, such as pure plant oils (PPOs). The first part of the study concerns application of a novel air-assisted impinging jet atomization for continuous spray applications, and the second part concerns transient spray applications. The particle/droplet imaging analysis (PDIA) technique along with direct imaging methods are used for the purpose of spray characterization. In the first part, effective atomization of Jatropha PPO is demonstrated at gas-to-liquid ratios (GLRs) on the order 0.1. The effect of liquid and gas flow rates on the spray characteristics is evaluated, and results indicate a Sauter mean diameter (SMD) of 50 mu m is achieved with GLRs as low as 0.05. In the second part of the study, a commercially available air-assisted transient atomizer is evaluated using Jatropha PPO. The effect of the pressure difference across the air injector and ambient gas pressure on liquid spray characteristics is studied. The results indicate that it is possible to achieve the same level of atomization of Jatropha as diesel fuel by operating the atomizer at a higher pressure difference. Specifically, a SMD of 44 mu m is obtained for the Jatropha oil using injection pressures of <1 MPa. A further interesting observation associated with this injector is the near constancy of a nondimensional spray penetration rate for the Jatropha oil spray.