Vaporization and collision modeling of liquid fuel sprays in a co-axial fuel and air pre-mixer

Droplet collision occurs frequently in regions where the droplet number density is high. Even for Lean Premixed and Pre-vaporized (LPP) liquid sprays, the collision effects can be very high on the droplet size distributions, which will in turn affect the droplet vaporization process. Hence, in conjunction with vaporization modeling, collision modeling for such spray systems is also essential. The standard O'Rourke's collision model, usually implemented in CFD codes, tends to generate unphysical numerical artifact when simulations are performed on Cartesian grid and the results are not grid independent. Thus, a new collision modeling approach based on no-time-counter method (NTC) proposed by Schmidt and Rutland is implemented to replace O'Rourke's collision algorithm to solve a spray injection problem in a cylindrical coflow premixer. The so called ``four-leaf clover'' numerical artifacts are eliminated by the new collision algorithm and results from a diesel spray show very good grid independence. Next, the dispersion and vaporization processes for liquid fuel sprays are simulated in a coflow premixer. Two liquid fuels under investigation are jet-A and Rapeseed Methyl Esters (RME). Results show very good grid independence in terms of SMD distribution, droplet number distribution and fuel vapor mass flow rate. A baseline test is first established with a spray cone angle of 90 degrees and injection velocity of 3 m/s and jet-A achieves much better vaporization performance than RME due to its higher vapor pressure. To improve the vaporization performance for both fuels, a series of simulations have been done at several different combinations of spray cone angle and injection velocity. At relatively low spray cone angle and injection velocity, the collision effect on the average droplet size and the vaporization performance are very high due to relatively high coalescence rate induced by droplet collisions. Thus, at higher spray cone angle and injection velocity, the results expectedly show improvement in fuel vaporization performance since smaller droplet has a higher vaporization rate. The vaporization performance and the level of homogeneity of fuel-air mixture can be significantly improved when the dispersion level is high, which can be achieved by increasing the spray cone angle and injection velocity. (C) 2012 Elsevier Ltd. All rights reserved.

Rack Level Modeling of Air Flow Through Perforated Tile in a Data Center

Effective air flow distribution through perforated tiles is required to efficiently cool servers in a raised floor data center. We present detailed computational fluid dynamics (CFD) modeling of air flow through a perforated tile and its entrance to the adjacent server rack. The realistic geometrical details of the perforated tile, as well as of the rack are included in the model. Generally, models for air flow through perforated tiles specify a step pressure loss across the tile surface, or porous jump model based on the tile porosity. An improvement to this includes a momentum source specification above the tile to simulate the acceleration of the air flow through the pores, or body force model. In both of these models, geometrical details of tile such as pore locations and shapes are not included. More details increase the grid size as well as the computational time. However, the grid refinement can be controlled to achieve balance between the accuracy and computational time. We compared the results from CFD using geometrical resolution with the porous jump and body force model solution as well as with the measured flow field using particle image velocimetry (PIV) experiments. We observe that including tile geometrical details gives better results as compared to elimination of tile geometrical details and specifying physical models across and above the tile surface. A modification to the body force model is also suggested and improved results were achieved.

Smoldering combustion of "incense" sticks - Experiments and modeling

This paper is concerned with the experimental and modeling studies on the smoldering rates of incense sticks as a function of ambient oxygen fraction in air, the flow velocity and size. The experimental results are obtained both for forward and reverse smolder conditions. The results are explained on the basis of surface combustion due to diffusion of oxygen to the surface by both free and forced convection supporting the heat transfer into the solid by conduction, into the stream by convection and the radiant heat transfer from the surface. The heat release at the surface is controlled by the convective transport of the oxidizer to the surface. To obtain the diffusion rates particularly for the reverse smolder, CFD calculations of fluid flow with along with a passive scalar are needed; these calculations have been made both for forward and reverse smolder. The interesting aspect of the CFD calculations is that while the Nusselt umber for forward smolder shows a clear root( Re-u) dependence ( Re-u = Flow Reynolds Number), the result for reverse smolder shows a peak in the variation with Reynolds number with the values lower than for forward smolder and unsteadiness in the flow beyond a certain flow rate. The results of flow behavior and Nusselt number are used in a simple model for the heat transfer at the smoldering surface to obtain the dependence of the smoldering rate on the diameter of the incense stick, the flow rate of air and the oxygen fraction. The results are presented in terms of a correlation for the non-dimensional smoldering rate with radiant flux from the surface and heat generation rate at the surface. The correlations appear reasonable for both forward and reverse smolder cases.

Wave dispersion in a cellular composite with modulated microstructure

In this paper we report a modeling technique and analysis of wave dispersion in a cellular composite laminate with spatially modulated microstructure, which can be modeled by parameterization and homogenization in an appropriate length scale. Higher order beam theory is applied and the system of wave equations are derived. Homogenization of these equations are carried out in the scale of wavelength and frequency of the individual wave modes. Smaller scale scattering below the order of cell size are filtered out in the present approach. The longitudinal dispersion relations for different values of a modulation parameter are analyzed which indicates the existence of stop and pass band patterns. Dispersion relations for flexural-shear case are also analyzed which indicates a tendency toward forming the stop and pass bands for increasing values of a shear stiffness modulation parameter. The effect the phase angle (θ) of the incident wave indicates the existence more number of alternative stop bands and pass bands for θ = 45°.

Multi-Mode Phonon-Controlled Field Emission from Carbon Nanotubes: Modeling and Experiments

The main idea proposed in this paper is that in a vertically aligned array of short carbon nanotubes (CNTs) grown on a metal substrate, we consider a frequency dependent electric field, so that the mode-specific propagation of phonons, in correspondence with the strained band structure and the dispersion curves, take place. We perform theoretical calculations to validate this idea with a view of optimizing the field emission behavior of the CNT array. This is the first approach of its kind, and is in contrast to the the conventional approach where a DC bias voltage is applied in order to observe field emission. A first set of experimental results presented in this paper gives a clear indication that phonon-assisted control of field emission current in CNT based thin film diode is possible.

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.

Modeling stochastic hybrid systems

Stochastic hybrid systems arise in numerous applications of systems with multiple models; e.g., air traffc management, flexible manufacturing systems, fault tolerant control systems etc. In a typical hybrid system, the state space is hybrid in the sense that some components take values in a Euclidean space, while some other components are discrete. In this paper we propose two stochastic hybrid models, both of which permit diffusion and hybrid jump. Such models are essential for studying air traffic management in a stochastic framework.

Unified Dispersion Characteristics of Structural Acoustic Waveguides

In this article, we show with some formalism that infinite flexible structural acoustic waveguides have a general form for the dispersion equation. The dispersion equation of all such waveguides should conform to a generic form. This allows us to bring out the common features of structural acoustic waveguides. We take three examples to demonstrate this fact, namely, the rectangular, the circular cylindrical and the elliptical geometries. Where necessary, the equations are simplified for applicability to a particular frequency-regime before demonstrating the conformance to the generic form of the dispersion relation. It is then shown that the coupled wavenumber solutions of all these systems can be represented on a single schematic.

Monodispersity and stability: case of ultrafine aluminium nanoparticles (< 5 nm) synthesized by the solvated metal atom dispersion approach

The synthesis of THF coordinated aluminium nanoparticles by the solvated metal atom dispersion (SMAD) method is described. These colloids are not stable with respect to precipitation of aluminium nanoparticles. The precipitated aluminium nanopowder is highly pyrophoric. Highly monodisperse colloidal aluminium nanoparticles (3.1 +/- 0.6 nm) stabilized by a capping agent, hexadecyl amine (HDA), have also been prepared by the SMAD method. They are stable towards precipitation of particles for more than a week. The Al-HDA nanoparticles are not as pyrophoric as the Al-THF samples. Particles synthesized in this manner were characterized by high-resolution electron microscopy and powder X-ray diffraction. Annealing of the Al-HDA nanoparticles resulted in carbonization of the capping agent on the surface of the particles which imparts air stability to them. Carbonization of the capping agent was established using Raman spectroscopy and TEM. The annealed aluminium nanoparticles were found to be stable even upon their exposure to air for over a month which was evident from the powder XRD, TGA/DSC, and TEM studies. The successful passivation was further confirmed with the determination of high active aluminium content (95 wt%) upon exposure and storage under air.

Modeling and optimization of a natural convection-based agricultural dryer

This study presents development of a computational fluid dynamic (CFD) model to predict unsteady, two-dimensional temperature, moisture and velocity distributions inside a novel, biomass-fired, natural convection-type agricultural dryer. Results show that in initial stages of drying, when material surface is wet and moisture is easily available, moisture removal rate from surface depends upon the condition of drying air. Subsequently, material surface becomes dry and moisture removal rate is driven by diffusion of moisture from inside to the material surface. An optimum 9-tray configuration is found to be more efficient than for the same mass of material and volume of dryer. A new configuration of dryer, mainly to explore its potential to increasing uniformity in drying across all trays, is also analyzed. This configuration involves diverting a portion of hot air before it enters over the first tray and is supplied directly at an intermediate location in the dryer. Uniformity in drying across trays has increased for the kind of material simulated.

Premixed flame response to equivalence ratio fluctuations: Comparison between reduced order modeling and detailed computations

Combustion instability events in lean premixed combustion systems can cause spatio-temporal variations in unburnt mixture fuel/air ratio. This provides a driving mechanism for heat-release oscillations when they interact with the flame. Several Reduced Order Modelling (ROM) approaches to predict the characteristics of these oscillations have been developed in the past. The present paper compares results for flame describing function characteristics determined from a ROM approach based on the level-set method, with corresponding results from detailed, fully compressible reacting flow computations for the same two dimensional slot flame configuration. The comparison between these results is seen to be sensitive to small geometric differences in the shape of the nominally steady flame used in the two computations. When the results are corrected to account for these differences, describing function magnitudes are well predicted for frequencies lesser than and greater than a lower and upper cutoff respectively due to amplification of flame surface wrinkling by the convective Darrieus-Landau (DL) instability. However, good agreement in describing function phase predictions is seen as the ROM captures the transit time of wrinkles through the flame correctly. Also, good agreement is seen for both magnitude and phase of the flame response, for large forcing amplitudes, at frequencies where the DL instability has a minimal influence. Thus, the present ROM can predict flame response as long as the DL instability, caused by gas expansion at the flame front, does not significantly alter flame front perturbation amplitudes as they traverse the flame. (C) 2012 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Dispersion characteristics of paired waveguide modes in 2D photonic band gap structures

In this paper the case of a typical line defect in 2D Photonic crystal is analyzed. The 2D photonic crystals are of dielectric rods in air in square and triangular lattice configurations. This line defect serves as waveguide with a pair of modes having opposite dispersion characteristics.

Energy Hyperspace for Stacking Interaction in AU/AU Dinucleotide Step: Dispersion-Corrected Density Functional Theory Study

Double helical structures of DNA and RNA are mostly determined by base pair stacking interactions, which give them the base sequence-directed features, such as small roll values for the purine-pyrimidine steps. Earlier attempts to characterize stacking interactions were mostly restricted to calculations on fiber diffraction geometries or optimized structure using ab initio calculations lacking variation in geometry to comment on rather unusual large roll values observed in AU/AU base pair step in crystal structures of RNA double helices. We have generated stacking energy hyperspace by modeling geometries with variations along the important degrees of freedom, roll, and slide, which were chosen via statistical analysis as maximally sequence dependent. Corresponding energy contours were constructed by several quantum chemical methods including dispersion corrections. This analysis established the most suitable methods for stacked base pair systems despite the limitation imparted by number of atom in a base pair step to employ very high level of theory. All the methods predict negative roll value and near-zero slide to be most favorable for the purine-pyrimidine steps, in agreement with Calladine's steric clash based rule. Successive base pairs in RNA are always linked by sugar-phosphate backbone with C3-endo sugars and this demands C1-C1 distance of about 5.4 angstrom along the chains. Consideration of an energy penalty term for deviation of C1-C1 distance from the mean value, to the recent DFT-D functionals, specifically B97X-D appears to predict reliable energy contour for AU/AU step. Such distance-based penalty improves energy contours for the other purine-pyrimidine sequences also. (c) 2013 Wiley Periodicals, Inc. Biopolymers 101: 107-120, 2014.

Morphing fixed wing MAV modeling using VAM

The design and implementation of a morphing Micro Air Vehicle (MAV) wing using a smart composite is attempted in this research work. Control surfaces actuated by traditional servos are difficult to instrument and fabricate on thin composite-wings of MAVs. Piezoelectric Fiber Reinforced Composites (PFRCs) are the chosen smart structural materials in the current work for incorporation onto fixed-wing MAVs to simultaneously perform the dual functions of structural load-bearing and actuation of flexure, torsion and/or extension for morphing. Further, PFRC use can be extended towards shape control of a “fixed” wing MAV to meet changing performance requirements. Wings that can warp into desired shapes and/or have variable camber are well-known to exhibit improved efficiency in aerodynamic control. During an entire flight cycle, there are multiple optimal configurations, each of which suits a particular phase of the flight regime. Widely proposed methods of wing morphing include changes in camber, twist, sweep and span. However, camber change during flight is already established, in terms of its potential, as a major factor in improving the aerofoil efficiency and flow separation behavior. Hence, for this work, morphing with camber change is adopted with the goal to better tailor aerodynamic properties.