ASYMPTOTIC PRESERVING SCHEME FOR A KINETIC MODEL DESCRIBING IN COMPRESSIBLE FLUIDS

The kinetic theory of fluid turbulence modeling developed by Degond and Lemou in 7] is considered for further study, analysis and simulation. Starting with the Boltzmann like equation representation for turbulence modeling, a relaxation type collision term is introduced for isotropic turbulence. In order to describe some important turbulence phenomenology, the relaxation time incorporates a dependency on the turbulent microscopic energy and this makes difficult the construction of efficient numerical methods. To investigate this problem, we focus here on a multi-dimensional prototype model and first propose an appropriate change of frame that makes the numerical study simpler. Then, a numerical strategy to tackle the stiff relaxation source term is introduced in the spirit of Asymptotic Preserving Schemes. Numerical tests are performed in a one-dimensional framework on the basis of the developed strategy to confirm its efficiency.

Instabilities of a compressible stratified fluid in horizontal sheared motion

The stability characteristics of a Helmholtz velocity profile in a stably stratified, compressible atmosphere in the presence of a lower boundary are studied. A jump in the Brunt–Väisälä frequency is introduced and the level at which this jump occurs is assumed to be different from the shear zone, to simulate sharp temperature discontinuities in the atmosphere. The results are compared with those of Pellacani, Tebaldi, and Tosi and Lindzen and Rosenthal. In the present configuration, new unstable modes with larger growth rates are found. The wavelengths of the most unstable gravity waves for the parameters pertaining to observed cases of clear air turbulence agree quite closely with the experimental values. Physics of Fluids is copyrighted by The American Institute of Physics

Instabilities induced by variation of Brunt-Vaisala frequency in compressible stratified shear flows

The stability characteristics of a Helmholtz velocity profile in a stably stratified, compressible fluid in the presence of a lower rigid boundary are studied. A jump in the Brunt-Vaisala frequency at a level different from the shear zone is introduced and the variation of the Brunt-Vaisala frequency with respect to the vertical coordinate in the middle layer of the three-layered model is considered. An analytic solution in each of the layers is obtained, and the dispersion relation is solved numerically for parameters relevant to the model. The effect of shear in the lowermost layer of the three-layered model for a Boussinesq fluid is discussed. The results are compared with the earlier studies of Lindzen and Rosenthal, and Sachdev and Satya Narayanan. In the present model, new unstable modes with larger growth rates are obtained and the most unstable gravity wave modes are found to agree closely with the observed ones at various heights. Physics of Fluids is copyrighted by The American Institute of Physics.

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.

On the computation of two-dimensional compressible flow fields by the modified local stability scheme

The modified local stability scheme is applied to several two-dimensional problems—blunt body flow, regular reflection of a shock and lambda shock. The resolution of the flow features obtained by the modified local stability scheme is found to be better than that achieved by the other first order schemes and almost identical to that achieved by the second order schemes incorporating artificial viscosity. The scheme is easy for coding, consumes moderate amount of computer storage and time. The scheme can be advantageously used in place of second order schemes.

Long waves in inviscid compressible atmospheres: Isothermal atmosphere

Solitary waves and cnoidal waves have been found in an adiabatic compressible atmosphere which, under ambient conditions, has winds, and is isothermal. The theory is illustrated with an example for which the background wind is linearly increasing. It is found that the number of possible critical speeds of the flow depends crucially on whether the Richardson number is greater or less than one‐fourth.

Long waves in inviscid compressible atmospheres: Isothermal atmosphere

Solitary waves and cnoidal waves have been found in an adiabatic compressible atmosphere which, under ambient conditions, has winds, and is isothermal. The theory is illustrated with an example for which the background wind is linearly increasing. It is found that the number of possible critical speeds of the flow depends crucially on whether the Richardson number is greater or less than one‐fourth.

A finite element method for compressible viscous fluid flows

This work presents a mixed three-dimensional finite element formulation for analyzing compressible viscous flows. The formulation is based on the primitive variables velocity, density, temperature and pressure. The goal of this work is to present a `stable' numerical formulation, and, thus, the interpolation functions for the field variables are chosen so as to satisfy the inf-sup conditions. An exact tangent stiffness matrix is derived for the formulation, which ensures a quadratic rate of convergence. The good performance of the proposed strategy is shown in a number of steady-state and transient problems where compressibility effects are important such as high Mach number flows, natural convection, Riemann problems, etc., and also on problems where the fluid can be treated as almost incompressible. Copyright (C) 2010 John Wiley & Sons, Ltd.

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.

Effect of side confining walls on the growth rate of compressible mixing layers

Model free simulations are performed to study the effect of the presence of side wall in compressible mixing of two parallel dissimilar gaseous streams with significant temperature difference. The turbulence statistics shows the three dimensional nature of the flow with and without the presence of side walls. The presence of side wall neither makes the flow field two dimensional, nor suppresses three dimensional disturbances. However, the comparison of shear layer growth rate and wall pressures reveal a better match with the two dimensional simulation results. This better match is explained on the basis of formation of oblique structures due to the presence of side walls which also suppress the distribution of momentum in third direction making the pressures to be higher as compared with the case without side walls. (C) 2013 Elsevier Ltd. All rights reserved.

Unsteady Compressible Boundary-Layer Flow On The Stagnation Line Of An Infinite Swept Cylinder

Numerical solutions of flow and heat transfer process on the unsteady flow of a compressible viscous fluid with variable gas properties in the vicinity of the stagnation line of an infinite swept cylinder are presented. Results are given for the case where the unsteady temperature field is produced by (i) a sudden change in the wall temperature (enthalpy) as the impulsive motion is started and (ii) a sudden change in the free-stream velocity. Solutions for the simultaneous development of the thermal and momentum boundary layers are obtained by using quasilinearization technique with an implicit finite difference scheme. Attention is given to the transient phenomenon from the initial flow to the final steady-state distribution. Results are presented for the skin friction and heat transfer coefficients as well as for the velocity and enthalpy profiles. The effects of wail enthalpy parameter, sweep parameter, fluid properties and transpiration cooling on the heat transfer and skin friction are considered.

Temperature dependent vibrational lifetimes in supercritical fluids near the critical point

Vibrational relaxation measurements on the CO asymmetric stretching mode (similar to 1980 cm(-1)) of tungsten hexacarbonyl (W(CO)(6)) as a function of temperature at constant density in several supercritical solvents in the vicinity of the critical point are presented. In supercritical ethane, at the critical density, there is a region above the critical temperature (Tc) in which the lifetime increases with increasing temperature. When the temperature is raised sufficiently (similar to T-c + 70 degrees C), the lifetime decreases with further increase in temperature. A recent hydrodynamic/thermodynamic theory of vibrational relaxation in supercritical fluids reproduces this behavior semiquantitatively. The temperature dependent data for fixed densities somewhat above and below the critical density is in better agreement with the theory. In fluoroform solvent at the critical density, the vibrational lifetime also initially increases with increasing temperature. However, in supercritical CO2 at the critical density, the temperature dependent vibrational lifetime decreases approximately linearly with temperature beginning almost immediately above T-c. The theory does not reproduce this behavior. A comparison between the absolute lifetimes in the three solvents and the temperature trends is made.

Mixed convection in non-Newtonian fluids along a vertical plate in a porous medium

The problem of mixed convection from vertical surfaces in a porous medium saturated with a power-law type non-Newtonian fluid is investigated. The transformed conservation laws are solved numerically for the case of variable wall temperature conditions. Results for the details of the velocity and temperature fields as well as the Nusselt number have been presented. The viscosity index ranged from 0.5-2.0.

First donor-acceptor interaction promoted gelation of organic fluids

In recent years there has been considerable interest in developing new types of gelators of organic solvents.1 Despite the recent advances, a priori design of a gelator for gelling a given solvent has remained a challenging task. Various noncovalent interactions like hydrogen-bonding,2 metal coordination3 etc. have been used as the driving force for the gelation process. A special class of cholesterol-based gelators were reported by Weiss,4 and by Shinkai.5 Gels derived from these molecules have been used for chiral recognition/sensing,6 for studying photo- and metal-responsive functions,7 and as templates to make hollow fiber silica.8 Other types of organogels have been used for designing polymerized 9 and reverse aerogels,10 and in molecular imprinting.11 Hanabusa’s group has recently reported organogels with a bile acid derivative.12 This has prompted us to disclose our results on a novel electron donor–acceptor (EDA) interaction mediated two-component13 gelator system based on the bile acid14 backbone.