Ayyaswamy, Portonovo S.

Email Address

ORCID

Disciplines

relationships.isProjectOf

relationships.isOrgUnitOf

Position

Introduction

Research Interests

Search Results

Now showing 1 - 10 of 22
  • Publication
    Numerical study of wall effects on buoyant gas-bubble rise in a liquid-filled finite cylinder
    (2007-09-01) Mukundakrishnan, Karthik; Ayyaswamy, Portonovo S; Eckmann, David M; Quan, Shaoping
    The wall effects on the axisymmetric rise and deformation of an initially spherical gas bubble released from rest in a liquid-filled, finite circular cylinder are numerically investigated. The bulk and gas phases are considered incompressible and immiscible. The bubble motion and deformation are characterized by the Morton number Mo, Eötvös number Eo, Reynolds number Re, Weber number We, density ratio, viscosity ratio, the ratios of the cylinder height and the cylinder radius to the diameter of the initially spherical bubble (H* =H/d0, R*=R/d0). Bubble rise in liquids described by Eo and Mo combinations ranging from (1,0.01) to (277.5,0.092), as appropriate to various terminal state Reynolds numbers (ReT) and shapes have been studied. The range of terminal state Reynolds numbers includes 0.02T<70. Bubble shapes at terminal states vary from spherical to intermediate spherical-cap–skirted. The numerical procedure employs a front tracking finite difference method coupled with a level contour reconstruction of the front. This procedure ensures a smooth distribution of the front points and conserves the bubble volume. For the wide range of Eo and Mo examined, bubble motion in cylinders of height H*=8 and R≥3, is noted to correspond to the rise in an infinite medium, both in terms of Reynolds number and shape at terminal state. In a thin cylindrical vessel (small R*) the motion of the bubble is retarded due to increased total drag and the bubble achieves terminal conditions within a short distance from release. The wake effects on bubble rise are reduced, and elongated bubbles may occur at appropriate conditions. For a fixed volume of the bubble, increasing the cylinder radius may result in the formation of well-defined rear recirculatory wakes that are associated with lateral bulging and skirt formation. The paper includes figures of bubble shape regimes for various values of R*, Eo, Mo, and ReT. Our predictions agree with existing results reported in the literature.
  • Publication
    Effect of Negative Ions on Electrical Breakdown in a Nonuniform Air Gap Between a Wire and a Plane
    (1994-01-24) Ramakrishna, K.; Cohen, Ira M.; Ayyaswamy, Portonovo S.
    Electrical breakdown of an axisymmetric, atmospheric pressure air gap between a wire and a plane has been investigated for a gap length of 0.5 mm. 0- and 02- have been identified as the negative ions affecting the discharge development in air, besides electrons and positive ions, and have been included in the electrical breakdown model. Five coupled two-dimensional transient partial differential equations describing the discharge evolution in the air gap have been solved using a finite difference algorithm developed earlier. Temporal development of the charged particle number densities, electrostatic potential, electric field, and current at both the electrodes is presented when the wire is negatively biased at 2500 V. The impact of negative ions on gap breakdown has been assessed by comparing the results of analyses with and without negative ions. It is concluded that the negative ions have negligible effect during the early stages of the discharge development. However, as the discharge evolves, the negative ions cause a net loss of electrons from the discharge. The effect is most pronounced away from the discharge axis, where peaks in the electron density occur as breakdown proceeds. Radial spread of discharge and current growth rate are relatively unaffected by the presence of negative ions, but the magnitude of total current at the electrodes has been found to decrease by a decade when the negative ions are present.
  • Publication
    Two-Dimenslonal Analysis of Electrical Breakdown in a Nonuniform Gap Between a Wire and a Plane
    (1989) Ramakrishna, K.; Cohen, Ira M.; Ayyaswamy, Portonovo S.
    Electrical breakdown of a gap between a wire (modeled as a hyperboloid) and a plane has been investigated numerically by solving the two-dimensional form of the diffusion flux equations for the charged particle number densities and Poisson's equation for the self-consistent electric field. Electron impact ionization, thermal ionization, and three-body recombination have been considered as the charged particle production and loss mechanisms. The electrode surfaces are considered to be absorbing and the initial density of the particles is small, but nonzero, A gap length of 0.5 mm is investigated and the gas medium is air or argon at atmospheric pressure. The temporal development of the profiles of ion and electron number densities, potential and electric field, and current growth on both the electrodes are presented when the applied voltage is 1500 and 2500 V for both positive and negative wires. When the wire is negatively biased, the peaks in the radial distribution of both of the charged particle densities near the wire occur off the axis except during the very early part of the breakdown. With positive polarity, the electron density maximum always occurs on the discharge axis, while for ions it moves away from the axis, later in the transient, due to the reverse particle drift in the electric field from the negative polarity case, The discharge spreads farther out into the ambient (almost two times the gap length) when the wire is negatively biased than with positive polarity. The effect of charge separation on the externally applied electric field is significant at voltages 2500 V and higher. Ionization is greater in argon than in air for a fixed potential difference between the electrodes.
  • Publication
    Breakdown of a Wire-to-Plane Discharge: Transient Effects
    (1991-08-12) Jog, Milind A.; Cohen, Ira M.; Ayyaswamy, Portonovo S.
    A wire-to-plane discharge during the early phases of breakdown has been studied. The discharge has been modeled in a prolate spheroidal coordinate system with the wire shape taken as a hyperboloid of revolution. Four simultaneous coupled, time-dependent, nonlinear partial differential equations describe the electrical discharge. These are the conservation equations for ion and electron densities, the energy equation for electron temperature, and Poisson’s equation for the self-consistent electric field. By solving this formulation subject to appropriate initial and boundary conditions, charged particle densities and temperature variations have been obtained as the ionization progresses in the discharge. The results show that both the electron temperature and the charged particle densities increase with the progress of ionization. The effect of different wire polarities is also examined. With a positive wire polarity, the increases in electron temperature and charged particle densities are confined to regions of the discharge in the vicinity of the wire tip. With a negative wire polarity, the breakdown occurs in the entire gap at a faster rate than with a positive wire polarity. The wire polarity affects the magnitude of energy transfer between the particles.
  • Publication
    Laminar Condensation on a Moving Drop. Part 2. Numerical Solutions
    (1984-02-01) Chung, J. N.; Ayyaswamy, Portonovo S.; Sadhal, Satwindar S.
    In this paper, we investigate the problem of transient laminar condensation on a moving drop by the semianalytical series-truncation method. The objectives are to assess the validity and the accuracy of the matched-asymptotic method employed in Part 1 . The fluid flow and thermodynamic variables are expanded as complete series of Legendre polynomials. The resulting transient momentum, energy and species equations are integrated numerically. The numerical scheme basically involves a three-point central difference for the spatial derivatives and a backward difference expression for the temporal derivatives. The finite-difference equations have been solved by the strongly implicit procedure. Good agreement of the fully transient numerical results with the singular perturbation approximation results of Part 1 lends credibility to a quasi-steady treatment of the continuous phase. The computational time requirements for the fully numerical solutions increase with decreasing non-condensable gas mass fraction in the bulk environment.
  • Publication
    Charged Particle Distributions and Heat Transfer in a Discharge Between Geometrically Dissimilar Electrodes: From Breakdown to Steady State
    (2000-02-01) Qin, Wei; Ayyaswamy, Portonovo S.; Cohen, Ira M.
    The low-current electric discharge from a fine wire anode to a planar cathode in atmospheric pressure air is numerically simulated from high-voltage prebreakdown through electron temperature growth, then ionization and consequent current growth to steady state, limited by a ballast resistor in the external circuit. Conservation of number ~mass! for ions and electrons, Gauss’ law for the self-consistent electric field, and energy conservation for electrons have been solved from breakdown to steady state in a body fitted coordinate system generated specifically for these two geometrically dissimilar electrodes. To facilitate the discussion of the results, the discharge has been categorized under ~a! electron acceleration period, ~b! charged particle generation period, ~c! current increase and voltage drop period, and ~d! current and voltage stabilization period. Results are given for transient electron, ion, and temperature distributions in the gap as well as current growth and voltage drop across the gap. Heat flux from the discharge to the wire is calculated. The numerical simulations were compared with experiments performed under the same conditions on a wire bonding machine with very close correspondence.
  • Publication
    Oscillatory enhancement of the squeezing flow of yield stress fluids: A novel experimental result
    (1997-05-25) Zwick, K. J.; Ayyaswamy, Portonovo S.; Cohen, Ira M.
    The extrusion of a yield stress fluid from the space between two parallel plates is investigated experimentally. Oscillating the magnitude of the squeezing force about a mean value (F = f[1+αcos(ωt)]) was observed to significantly enhance the flow rate of yield stress fluids, while having no effect on the flow rate of Newtonian fluids. This is a novel result. The enhancement depends on the magnitude of the force, the oscillatory frequency and amplitude, the fluid being squeezed, and the thickness of the fluid layer. Non-dimensional results for the various flow quantities have been presented by using the flow predicted for the constant-force squeezing of a Herschel-Bulkley yield stress fluid as the reference. In the limit of constant-force squeezing, the present experimental results compare very well with those of our earlier theoretical model for this situation (Zwick, Ayyaswamy & Cohen 1996). The results presented in this paper have significance, among many applications, for injection moulding, in the adhesive bonding of microelectronic chips, and in surgical procedures employed in health care.
  • Publication
    Thin-Flame Theory for the Combustion of a Moving Liquid Drop: Effects Due to Variable Density
    (1986-10-01) Gogos, George; Sadhal, Satwindar S.; Ayyaswamy, Portonovo S; Sundararajan, T.
    The combustion of a moving liquid fuel drop has been investigated. The drop experiences a strong evaporation-induced radial velocity while undergoing slow translation. In view of the high evaporation velocity, the flow field is not in the Stokes regime. The combustion process is modelled by an indefinitely fast chemical reaction rate. While the flow and the transport in the continuous phase and the drop internal circulation are treated as quasisteady, the drop heat-up is regarded as a transient process. The transport equations of the continuous phase require analysis by a singular perturbation technique. The transient heat-up of the drop interior is solved by a series-truncation numerical method. The solution for the total problem is obtained by coupling the results for the continuous and dispersed phases. The enhancement in the mass burning rate and the deformation of the flame shape due to drop translation have been predicted. The initial temperature of the drop and the subsequent heating influence the temporal variations of the flamefront standoff ratio and the flame distance. The friction drag, the pressure drag and the drag due to interfacial momentum flux are individually predicted, and the total drag behaviour is discussed. The circulation inside the drop decreases with evaporation rate. A sufficiently large non-uniform evaporation velocity causes the circulation to reverse.
  • Publication
    Thermal and electrical characteristics of a two‐dimensional tanh‐conductivity arc
    (1978) Ayyaswamy, Portonovo S.; Das, G. C.; Cohen, Ira M.
    The two-dimensional variable-property arc has been studied through the use of the tanh-conductivity model. Results that describe the thermal and electric arc characteristics for various values of the electrode temperatures and aspect ratios are given. The numerical evaluation is carried out by the use of a Galerkin technique. The results exhibit several novel and interesting features depending on the arc parameters. For large aspect ratios (ratio of the interelectrode distance to that between the bounding walls) and small electrode temperatures, the current---electric-field characteristics tend toward those of a slender arc. However, at a given aspect ratio with large enough electrode temperatures, the distinct minimum noted in the slender-arc characteristics does not occur. Also, for a given aspect ratio and large enough differences in electrode potential, the electric-field-current characteristic is nearly linear and is independent of the electrode temperature. The transverse electrostatic potential is found to have no significant variation in cross-sectional planes. The qualitative nature of the thermal characteristics are similar to those of a constant-property arc although significant differences in quantitative results exist. Wall and electrode heat transfer rates are provided.
  • Publication
    Finite-sized gas bubble motion in a blood vessel: Non-Newtonian effects
    (2008-09-01) Mukundakrishnan, Karthik; Ayyaswamy, Portonovo S; Eckmann, David M
    We have numerically investigated the axisymmetric motion of a finite-sized nearly occluding air bubble through a shear-thinning Casson fluid flowing in blood vessels of circular cross section. The numerical solution entails solving a two-layer fluid model - a cell-free layer and a non-Newtonian core together with the gas bubble. This problem is of interest to the field of rheology and for gas embolism studies in health sciences. The numerical method is based on a modified front-tracking method. The viscosity expression in the Casson model for blood (bulk fluid) includes the hematocrit [the volume fraction of red blood cells (RBCs)] as an explicit parameter. Three different flow Reynolds numbers, Reapp=ΡlUmaxd/µapp, in the neighborhood of 0.2, 2, and 200 are investigated. Here, Ρl is the density of blood, Umax is the centerline velocity of the inlet Casson profile, d is the diameter of the vessel, and µapp is the apparent viscosity of whole blood. Three different hematocrits have also been considered: 0.45, 0.4, and 0.335. The vessel sizes considered correspond to small arteries, and small and large arterioles in normal humans. The degree of bubble occlusion is characterized by the ratio of bubble to vessel radius (aspect ratio), λ, in the range 0.9 ≤ λ≤1.05. For arteriolar flow, where relevant, the Fahraeus-Lindqvist effects are taken into account. Both horizontal and vertical vessel geometries have been investigated. Many significant insights are revealed by our study: (i) bubble motion causes large temporal and spatial gradients of shear stress at the "endothelial cell" (EC) surface lining the blood vessel wall as the bubble approaches the cell, moves over it, and passes it by; (ii) rapid reversals occur in the sign of the shear stress (+ → - → +) imparted to the cell surface during bubble motion; (iii) large shear stress gradients together with sign reversals are ascribable to the development of a recirculation vortex at the rear of the bubble; (iv) computed magnitudes of shear stress gradients coupled with their sign reversals may correspond to levels that cause injury to the cell by membrane disruption through impulsive compression and stretching; and (v) for the vessel sizes and flow rates investigated, gravitational effects are negligible.