Modified kinetic flux vector splitting (m-KFVS) method for compressible flows

To resolve many flow features accurately, like accurate capture of suction peak in subsonic flows and crisp shocks in flows with discontinuities, to minimise the loss in stagnation pressure in isentropic flows or even flow separation in viscous flows require an accurate and low dissipative numerical scheme. The first order kinetic flux vector splitting (KFVS) method has been found to be very robust but suffers from the problem of having much more numerical diffusion than required, resulting in inaccurate computation of the above flow features. However, numerical dissipation can be reduced by refining the grid or by using higher order kinetic schemes. In flows with strong shock waves, the higher order schemes require limiters, which reduce the local order of accuracy to first order, resulting in degradation of flow features in many cases. Further, these schemes require more points in the stencil and hence consume more computational time and memory. In this paper, we present a low dissipative modified KFVS (m-KFVS) method which leads to improved splitting of inviscid fluxes. The m-KFVS method captures the above flow features more accurately compared to first order KFVS and the results are comparable to second order accurate KFVS method, by still using the first order stencil. (C) 2011 Elsevier Ltd. All rights reserved.

A note on the damping and oscillations of a fluid drop moving in another fluid

The oscillations of a drop moving in another fluid medium have been studied at low values of Reynolds number and Weber number by taking into consideration the shape of the drop and the viscosities of the two phases in addition to the interfacial tension. The deformation of the drop modifies the Lamb's expression for frequency by including a correction term while the viscous effects split the frequency into a pair of frequencies—one lower and the other higher than Lamb's. The lower frequency mode has ample experimental support while the higher frequency mode has also been observed. The two modes almost merge with Lamb's frequency for the asymptotic cases of a drop in free space or a bubble in a dense viscous fluid but the splitting becomes large when the two fluids have similar properties. Instead of oscillations, aperiodic damping modes are found to occur in drops with sizes smaller than a critical size ($\sim\hat{\rho}\hat{\nu}^2/T $). With the help of these calculations, many of the available experimental results are analyzed and discussed.

Stability of the flow of a fluid through a flexible tube at high Reynolds number

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (ρVR / η), the ratio of the viscosities of the wall and fluid ηr = (ηs/η), the ratio of radii H and the dimensionless velocity Γ = (ρV2/G)1/2. Here ρ is the density of the fluid, G is the coefficient of elasticity of the wall and Vis the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter ε = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate do), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctruations due to the Reynolds strees. There is an O(ε1/2) correction to the growth rate, s(1), due to the presence of a wall layer of thickness ε1/2R where the viscous stresses are O(ε1/2) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Γ and wavenumber k where s(l) = 0. At these points, the wail layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(ε) correction to the growth rate s(2) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s(2) increases [is proportional to] (H − 1)−2 for (H − 1) [double less-than sign] 1 (thickness of wall much less than the tube radius), and decreases [is proportional to] (H−4 for H [dbl greater-than sign] 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube.

The collective dynamics of self-propelled particles

We propose a method for the dynamic simulation of a collection of self-propelled particles in a viscous Newtonian fluid. We restrict attention to particles whose size and velocity are small enough that the fluid motion is in the creeping flow regime. We propose a simple model for a self-propelled particle, and extended the Stokesian Dynamics method to conduct dynamic simulations of a collection of such particles. In our description, each particle is treated as a sphere with an orientation vector p, whose locomotion is driven by the action of a force dipole Sp of constant magnitude S0 at a point slightly displaced from its centre. To simplify the calculation, we place the dipole at the centre of the particle, and introduce a virtual propulsion force Fp to effect propulsion. The magnitude F0 of this force is proportional to S0. The directions of Sp and Fp are determined by p. In isolation, a self-propelled particle moves at a constant velocity u0 p, with the speed u0 determined by S0. When it coexists with many such particles, its hydrodynamic interaction with the other particles alters its velocity and, more importantly, its orientation. As a result, the motion of the particle is chaotic. Our simulations are not restricted to low particle concentration, as we implement the full hydrodynamic interactions between the particles, but we restrict the motion of particles to two dimensions to reduce computation. We have studied the statistical properties of a suspension of self-propelled particles for a range of the particle concentration, quantified by the area fraction φa. We find several interesting features in the microstructure and statistics. We find that particles tend to swim in clusters wherein they are in close proximity. Consequently, incorporating the finite size of the particles and the near-field hydrodynamic interactions is of the essence. There is a continuous process of breakage and formation of the clusters. We find that the distributions of particle velocity at low and high φa are qualitatively different; it is close to the normal distribution at high φa, in agreement with experimental measurements. The motion of the particles is diffusive at long time, and the self-diffusivity decreases with increasing φa. The pair correlation function shows a large anisotropic build-up near contact, which decays rapidly with separation. There is also an anisotropic orientation correlation near contact, which decays more slowly with separation. Movies are available with the online version of the paper.

A Note On The Effect Of Short And Long Laminar Separation-Bubbles On Desinent Cavitation

Earlier desinent cavitation studies on a 1/8 caliber ogive by one of the authors (J. W. H.) showed a sudden change in the magnitude of the desinent cavitation number at a critical velocity. In the present work it is shown by means of oil-film flow visualization that below the critical velocity a long laminar separation bubble exists whereas above the critical velocity the laminar separation bubble is short. Thus the desinent cavitation characteristics of a 1/8 caliber ogive are governed by the nature of the viscous flow around the body.

Regimes of spray formation in gas-centered swirl coaxial atomizers

Spray formation in ambient atmosphere from gas-centered swirl coaxial atomizers is described by carrying out experiments in a spray test facility. The atomizer discharges a circular air jet and an axisymmetric swirling water sheet from its coaxially arranged inner and outer orifices. A high-speed digital imaging system along with a backlight illumination arrangement is employed to record the details of liquid sheet breakup and spray development. Spray regimes exhibiting different sheet breakup mechanisms are identified and their characteristic features presented. The identified spray regimes are wave-assisted sheet breakup, perforated sheet breakup, segmented sheet breakup, and pulsation spray regime. In the regime of wave-assisted sheet breakup, the sheet breakup shows features similar to the breakup of two-dimensional planar air-blasted liquid sheets. At high air-to-liquid momentum ratios, the interaction process between the axisymmetric swirling liquid sheet and the circular air jet develops spray processes which are more specific to the atomizer studied here. The spray exhibits a periodic ejection of liquid masses whose features are dominantly controlled by the central air jet.

Inception of Cavitation From a Backward Facing Step

Cavitation inception measurements are reported for flow past a downstream facing step with the height of the step varying from about 0.4 to 5 percent of the forebody diameter. The forebody was a 49 mm hemispherical nose and sigmai values were found to be very strong function of the height of the step. In addition, sigmai values were found to depend on whether the boundary layer approaching the step was laminar or turbulent. Generally sigmai values for turbulent case were lower.

Explicit filtering LES for turbulent non-premixed combustion

Two models for large eddy simulation of turbulent reacting flow in homogeneous turbulence were studied. The sub-grid stress arising out of non-linearities of the Navier-Stokes equations were modeled using an explicit filtering approach. A filtered mass density function (FMDF) approach was used for closure of the sub-grid scalar fluctuations. A posteriori calculations, when compared with the results from the direct numerical simulation, indicate that the explicit filtering is adequate in representing the effect of sub-grid stress on the filtered velocity field in the absence of reaction. Discrepancies arise when reactions occur, but the FMDF approach suffices to account for sub-grid scale fluctuations of the reacting scalars, accurately.

Filtered density function for large eddy simulation of turbulent reacting flows

A methodology termed the “filtered density function” (FDF) is developed and implemented for large eddy simulation (LES) of chemically reacting turbulent flows. In this methodology, the effects of the unresolved scalar fluctuations are taken into account by considering the probability density function (PDF) of subgrid scale (SGS) scalar quantities. A transport equation is derived for the FDF in which the effect of chemical reactions appears in a closed form. The influences of scalar mixing and convection within the subgrid are modeled. The FDF transport equation is solved numerically via a Lagrangian Monte Carlo scheme in which the solutions of the equivalent stochastic differential equations (SDEs) are obtained. These solutions preserve the Itô-Gikhman nature of the SDEs. The consistency of the FDF approach, the convergence of its Monte Carlo solution and the performance of the closures employed in the FDF transport equation are assessed by comparisons with results obtained by direct numerical simulation (DNS) and by conventional LES procedures in which the first two SGS scalar moments are obtained by a finite difference method (LES-FD). These comparative assessments are conducted by implementations of all three schemes (FDF, DNS and LES-FD) in a temporally developing mixing layer and a spatially developing planar jet under both non-reacting and reacting conditions. In non-reacting flows, the Monte Carlo solution of the FDF yields results similar to those via LES-FD. The advantage of the FDF is demonstrated by its use in reacting flows. In the absence of a closure for the SGS scalar fluctuations, the LES-FD results are significantly different from those based on DNS. The FDF results show a much closer agreement with filtered DNS results. © 1998 American Institute of Physics.

Influence of working fluid on the performance of a standing-wave thermoacoustic prime mover

The thermoacoustic prime mover (TAPM) is an attractive alternative as a pressure wave generator to drive Pulse Tube Cryocoolers (PTCs), by the absence of moving parts, construction simplicity, reasonable efficiency, and environmental friendly. Decreasing the resonance frequency and improving the efficiency of the TAPM are important to drive the PTCs. These are controlled by the working gas parameters other than the dimensions of TAPM. In this technical note, the experimental studies carried out to evaluate the influence of different working fluids on the performances of a twin standing wave TAPM at various operating pressures have been compared with the simulation studies of the same system using DeltaEc wherever possible. The reasonably good agreement between them indicates the utility of DeltaEc for the optimal design of TAPM with the right working fluids for practical applications. (C) 2011 Elsevier Ltd. All rights reserved.

Enhanced Survival of Probiotic Lactobacillus acidophilus by Encapsulation with Nanostructured Polyelectrolyte Layers through Layer-by-Layer Approach

The encapsulation of probiotic Lactobacillus acidophilus through layer-by-layer self-assembly of polyelectrolytes (PE) chitosan (CHI) and carboxymethyl cellulose (CMC) has been investigated,to enhance its survival m adverse conditions encountered in the GI tract The survival of encapsulated cells in simulated gastric (SGF) and intestinal fluids (SIF) is significant when compared to nonencapsulated cells On sequential exposure to SGF and SIF fox 120 nun, almost complete death of free cells is observed However, for cells coated with three nanolayers of PEs (CHI/CMC/CHI) about 33 log % of the cells (6 log cfu/500 mg) survived under the same conditions The enhanced survival rate of encapsulated L acidophilus can be attributed to the impermeability of polyelectrolyte nanolayers to large enzyme molecules like pepsin, and pancreatin that cause proteolysis and to the stability of the polyelectrolyte nanolayers in gastric and intestinal pH The PE coating also serves to reduce viability losses during freezing and freeze- drying About 73 and 92 log % of uncoated and coated cells survived after freeze:drying, and the losses occurring between freezing and freeze-drying were found to be lower for coated cells

Mode-coupling glass transition in a fluid confined by a periodic potential

We show that a fluid under strong spatially periodic confinement displays a glass transition within mode-coupling theory at a much lower density than the corresponding bulk system. We use fluctuating hydrodynamics, with confinement imposed through a periodic potential whose wavelength plays an important role in our treatment. To make the calculation tractable we implement a detailed calculation in one dimension. Although we do not expect simple 1d fluids to show a glass transition, our results are indicative of the behavior expected in higher dimensions. In a certain region of parameter space we observe a three-step relaxation reported recently in computer simulations [S. H. Krishnan, Ph.D. thesis, Indian Institute of Science (2005); Kim et al., Eur. Phys. J. Special Topics 189, 135 (2010)] and a glass-glass transition. We compare our results to those of Krakoviack [Phys. Rev. E 75, 031503 (2007)] and Lang et al. [Phys. Rev. Lett. 105, 125701 (2010)].

Viscous effects in the inception of cavitation

The inception of cavitation in the steady flow of liquids around bodies is seen to depend upon the real fluid flow around the bodies as well as the supply of nucleating cavitation sources—or nuclei—within the fluid. A primary distinction is made between bodies having a laminar separation or not having a laminar separation. The former group is relatively insensitive to the nuclei concentration whereas the latter is much more sensitive. Except for the case of fully separated wake flows and for gaseous cavitation by diffusion the cavitation inception index tends always to be less than the magnitude of the minimum pressure coefficient and only approaches that value for high Reynolds numbers in flows well supplied with nuclei.

2.5PN linear momentum flux from inspiralling compact binaries in quasicircular orbits and associated recoil: Nonspinning case

Anisotropic emission of gravitational waves (GWs) from inspiralling compact binaries leads to the loss of linear momentum and hence gravitational recoil of the system. The loss rate of linear momentum in the far-zone of the source (a nonspinning binary system of black holes in quasicircular orbit) is investigated at the 2.5 post-Newtonian (PN) order and used to provide an analytical expression in harmonic coordinates for the 2.5PN accurate recoil velocity of the binary accumulated in the inspiral phase. The maximum recoil velocity of the binary system at the end of its inspiral phase (i.e at the innermost stable circular orbit (ISCO)) estimated by the 2.5PN formula is of the order of 4 km s(-1) which is smaller than the 2PN estimate of 22 km s(-1). Going beyond inspiral, we also provide an estimate of the more important contribution to the recoil velocity from the plunge phase. The maximum recoil velocity at the end of the plunge, involving contributions both from inspiral and plunge phase, for a binary with symmetric mass ratio nu = 0.2 is of the order of 182 km s(-1).

A new 6-component accelerometer force balance for short duration ground testing facilities

A new six-component accelerometer force balance is developed and used in the HST2 shock tunnel of Indian Institute of Science. Aerodynamic forces and moments for a hypersonic slender body measured using this balance system at a free stream Mach number of 5.75 and Reynolds number of 1.5 million and stagnation enthalpy of 1.5 and 2 MJ/kg are presented. These measured values compare well with the theoretical values estimated using modified Newtonian theory.