R-parameter: A local truncation error based adaptive framework for finite volume compressible flow solvers

A residual-based strategy to estimate the local truncation error in a finite volume framework for steady compressible flows is proposed. This estimator, referred to as the -parameter, is derived from the imbalance arising from the use of an exact operator on the numerical solution for conservation laws. The behaviour of the residual estimator for linear and non-linear hyperbolic problems is systematically analysed. The relationship of the residual to the global error is also studied. The -parameter is used to derive a target length scale and consequently devise a suitable criterion for refinement/derefinement. This strategy, devoid of any user-defined parameters, is validated using two standard test cases involving smooth flows. A hybrid adaptive strategy based on both the error indicators and the -parameter, for flows involving shocks is also developed. Numerical studies on several compressible flow cases show that the adaptive algorithm performs excellently well in both two and three dimensions.

Pulsatile blood flow in a rigid pulmonary alveolar sheet with porous walls

This paper deals with the pulsatile blood flow in the lung alveolar sheets by idealizing each of them as a channel covered by porous media. As the blood flow in the lung is of low Reynolds number, a creeping flow is assumed in the channel. The analytical and numerical results for the velocity and pressure distribution in the porous medium are presented. The effect of an imposed slip condition is also studied. Comparisons with the corresponding results for the steady-state case are made at the end.

Dynamics of dense sheared granular flows. Part 1. Structure and diffusion

Shear flows of inelastic spheres in three dimensions in the Volume fraction range 0.4-0.64 are analysed using event-driven simulations.Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to the line Joining the centres). Here, we have considered both e(t) = +1 and e(t) = e(n) (rough particles) and e(t) =-1 (smooth particles), and the normal coefficient of restitution e(n) was varied in the range 0.6-0.98. Care was taken to avoid inelastic collapse and ensure there are no particle overlaps during the simulation. First, we studied the ordering in the system by examining the icosahedral order parameter Q(6) in three dimensions and the planar order parameter q(6) in the plane perpendicular to the gradient direction. It was found that for shear flows of sufficiently large size, the system Continues to be in the random state, with Q(6) and q(6) close to 0, even for volume fractions between phi = 0.5 and phi = 0.6; in contrast, for a system of elastic particles in the absence of shear, the system orders (crystallizes) at phi = 0.49. This indicates that the shear flow prevents ordering in a system of sufficiently large size. In a shear flow of inelastic particles, the strain rate and the temperature are related through the energy balance equation, and all time scales can be non-dimensionalized by the inverse of the strain rate. Therefore, the dynamics of the system are determined only by the volume fraction and the coefficients of restitution. The variation of the collision frequency with volume fraction and coefficient of estitution was examined. It was found, by plotting the inverse of the collision frequency as a function of volume fraction, that the collision frequency at constant strain rate diverges at a volume fraction phi(ad) (volume fraction for arrested dynamics) which is lower than the random close-packing Volume fraction 0.64 in the absence of shear. The volume fraction phi(ad) decreases as the coefficient of restitution is decreased from e(n) = 1; phi(ad) has a minimum of about 0.585 for coefficient of restitution e(n) in the range 0.6-0.8 for rough particles and is slightly larger for smooth particles. It is found that the dissipation rate and all components of the stress diverge proportional to the collision frequency in the close-packing limit. The qualitative behaviour of the increase in the stress and dissipation rate are well Captured by results derived from kinetic theory, but the quantitative agreement is lacking even if the collision frequency obtained from simulations is used to calculate the pair correlation function used In the theory.

Structure of shock waves at re-entry speeds

The unified structure of steady, one-dimensional shock waves in argon, in the absence of an external electric or magnetic field, is investigated. The analysis is based on a two-temperature, three-fluid continuum approach, using the Navier—Stokes equations as a model and including non-equilibrium collisional as well as radiative ionization phenomena. Quasi charge neutrality and zero velocity slip are assumed. The integral nature of the radiative terms is reduced to analytical forms through suitable spectral and directional approximations. The analysis is based on the method of matched asymptotic expansions. With respect to a suitably chosen small parameter, which is the ratio of atom-atom elastic collisional mean free-path to photon mean free-path, the following shock morphology emerges: within the radiation and electron thermal conduction dominated outer layer occurs an optically transparent discontinuity which consists of a chemically frozen heavy particle (atoms and ions) shock and a collisional ionization relaxation layer. Solutions are obtained for the first order with respect to the small parameter of the problem for two cases: (i) including electron thermal conduction and (ii) neglecting it in the analysis of the outer layer. It has been found that the influence of electron thermal conduction on the shock structure is substantial. Results for various free-stream conditions are presented in the form of tables and figures.

Flow of a micropolar fluid in a tube with sinusoidal deformation

The flow of a micropolar fluid for sinusoidally deforming boundaries is discussed in detail. The velocity and microrotation fields and the streamfunction are determined and plotted under the assumption of small deformations.

A Method of Characteristics for Solving Two-dimensional Unsteady Flow Problems

The aim of this investigation is to evolve a method of solving two-dimensional unsteady flow problems by the method of characteristics. This involves the reduction of the given system of equations to an equivalent system where only interior derivatives occur on a characteristic surface. From this system, four special bicharacteristic directional derivatives are chosen. A finite difference scheme is prescribed for solving the equations. General rectangular lattices are also considered. As an example, we investigate the propagation of an initial pressure distribution in a medium at rest.

Pulsatile flow in tubes of varying cross-sections

The pulsatile flow of an incompressible viscous fluid in a cylindrical tube of varying cross section is investigated for small Reynolds numbers. The solutions consist of a stedy and an oscillatory part. The shear stress distribution on the wall is evaluated and discussed in detail for special geometries like tapered tubes, locally constricted tubes and peristaltic tubes. The existence of separation in the flow field is noticed.

The non-planar peptide unit II. Comparison of theory with crystal structure datastar, open

The theoretical results derived in Part I (Ramachandran, G.N., Lakshminarayan, A.V. and Kolaskar, A.S. (1973) Biochim. Biophys. Acta 303, 8–13) that the three bonds of the peptide unit meeting at N can have a pyramidal structure is confirmed by an analysis of 14 published crystal structures of small peptides. It is shown that the dihedral angles θN and Δω are correlated, while θC, is small and is uncorrelated with Δω, showing that the non-planar distortion at C′ is generally small.

Theory of the non-planar peptide unitstar, open

By means of CNDO/2 calculations on N-methyl acetamide, it is shown that the state of minimum energy of the trans-peptide unit is a non-planar conformation, with the NH and NC2α bonds being significantly out of the plane formed by the atoms C1α, C′, O and N.

Uniformity of film thickness on rotating planetary planar substrates

For the purposes of obtaining a number of components with nearly identical thickness distributions over the substrate area and of minimizing the inhomogeneities of the film, it is logical to presume that a substrate rotating on its own axis and revolving around another axis will give more uniformity in film thickness than a substrate only revolving around one axis. In relation to the practical applications, an investigation has been undertaken to study the refinement that can be achieved by using a planar planetary substrate holder. It is shown theoretically that the use of the planetary substrate holder under ideal conditions of source and geometry does not offer any further improvement in uniformity of thickness over the conventional rotary work-holder. It is also shown that the geometrical parameters alone have little influence over the uniformity achieved on a planetary substrate, because of the complex cyclidal motion of any point on it. However, for any given geometry, a non-integral speed ratio of the planetary substrate and the work-holder shows considerably less variation in thickness over the substrate area.

Flow-graph techniques for epicyclic gear train analysis

Flow-graph techniques are applied in this article for the analysis of an epicyclic gear train. A gear system based on this is designed and constructed for use in Numerical Control Systems.

A parametric differentiation version with finite-difference scheme applicable to a class of problems in boundary layer flow with massive blowing

A numerical procedure, based on the parametric differentiation and implicit finite difference scheme, has been developed for a class of problems in the boundary-layer theory for saddle-point regions. Here, the results are presented for the case of a three-dimensional stagnation-point flow with massive blowing. The method compares very well with other methods for particular cases (zero or small mass blowing). Results emphasize that the present numerical procedure is well suited for the solution of saddle-point flows with massive blowing, which could not be solved by other methods.

Experimental and numerical investigation of jet injection in a wall bounded supersonic flow

Numerical and experimental studies of a supersonic jet (Helium) inclined at 45 degrees to a oncoming Mach 2 flow have been carried out. The numerical study has been used to arrive at a geometry that could reduce an oncoming Mach 5.75 flow to Mach 2 flow and in determining the jet parameters. Experiments are carried out in the IISc. hypersonic shock tunnel HST2 at similar conditions obtained from numerical studies. Flow visualization studies carried out using Schlieren technique clearly show the presence of the bow shock in front of the jet exposed to supersonic cross flow. The jet Mach number is experimentally found to be approximate to 3. Visual observations show that the jet has penetrated up to 60% of the total height of the chamber.

Experimental investigation of cowl shape and location on inlet characteristics at hypersonic Mach number

An experimental study is presented to show the effect of the cowl location and shape on the shock interaction phenomena in the inlet region for a 2D, planar scramjet inlet model. Investigations include schlieren visualization around the cowl region and heat transfer rate measurement inside the inlet chamber.Both regular and Mach reflections are observed when the forebody ramp shock reflects from the cowl plate. Mach stem heights of 3.3 mm and 4.1 mm are measured in 18.5 mm and 22.7 mm high inlet chambers respecively. Increased heat transfer rate is measured at the same location of chamber for cowls of longer lenghs is indicating additional mass flow recovery by the inlet.

Investigations of separated flow over backward facing steps in IISc hypersonic shock tunnel

Two backward facing step (2 mm and 3 mm step height) models are selected for surface heat transfer measurements. The platinum thin film gauges are deposited on the Macor inserts using both hand paint and vacuum sputtering technique. Using the Eckert reference temperature method the heating rates has been theoretically calculated along the flat plate portion of the model and the theoretical estimates are compared with experimentally determined surface heat transfer rate. Theoretical analysis of heat flux distribution down stream of the backward facing step model has been carried out using Gai’s non-dimensional analysis. Based on the measured surface heating rates on the backward facing step, the reattachment distance is estimated for 2 and 3 mm step height at nominal Mach number of 7.6. It has been found from the present study that for 2 and 3 mm step height, it approximately takes about 10 and 8 step heights downstream of the model respectively for the flow to re-attach.