Numerical analysis of the stability of nonisothermal couette flow

The stability of convective motion of a liquid between two rotating heated cylinders is investigated in the absence of external forces. The mathematical model for describing the convection is obtained from the general equations [1, 2] on the assumption that the density of the liquid, the thermal conductivity, the specific heat and the viscosity coefficients depend only on temperature, and that the work done by the pressure forces and the viscous dissipation are negligibly small. The thermal expansion coefficient of the liquid is not assumed to be small, which distinguishes the models in question from the classical Oberbeck-Boussinesq model [1, 3, 4]. Rostov-on-Don. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 70–76, September–October, 1988.     

Measurement of the concentration diffusion coefficient for Ne-Ar,Ne-Xe,Ne-H2, Xe-H2, H2-N2 and H2-O2 gas systems

The concentration diffusion coefficient, D ₁₂, is measured for the equimolar mixtures of Ne-Ar, Ne-Xe, Ne-H₂, Xe-H₂, H₂-N₂ and H₂-O₂ binary gas systems in a two-bulb metal apparatus in the temperature range 0 C to 100 C. These values are compared with the existing data on these systems and with the predictions of the kinetic theory in conjunction with the modified Buckingham exp-six potential. Unlike the thermal diffusion coefficient, with the simple theory it is possible to predict D ₁₂ within a few percent even for systems involving polyatomic gases. The smoothed experimental D ₁₂ values are also used to obtain data for the coefficients of viscosity and thermal conductivity at round temperatures and compositions for these systems.     

Chapman–Enskog solutions to arbitrary order in Sonine polynomials III: Diffusion,thermal diffusion,and thermal conductivity in a binary,rigid-sphere,gas mixture

The Chapman–Enskog solutions of the Boltzmann equation provide a basis for the computation of important transport coefficients for both simple gases and gas mixtures. These coefficients include the viscosity, the thermal conductivity, and the diffusion coefficient. In a preceding paper (I), for simple, rigid-sphere gases (i.e. single-component, monatomic gases) we have shown that the use of higher-order Sonine polynomial expansions enables one to obtain results of arbitrary precision that are free of numerical error and, in a second paper (II), we have extended our initial simple gas work to modeling the viscosity in a binary, rigid-sphere, gas mixture. In this latter paper we reported an extensive set of order 60 results which are believed to constitute the best currently available benchmark viscosity values for binary, rigid-sphere, gas mixtures. It is our purpose in this paper to similarly report the results of our investigation of relatively high-order (order 70), standard, Sonine polynomial expansions for the diffusion- and thermal conductivity-related Chapman–Enskog solutions for binary gas mixtures of rigid-sphere molecules. We note that in this work, as in our previous work, we have retained the full dependence of the solution on the molecular masses, the molecular sizes, the mole fractions, and the intermolecular potential model via the omega integrals. For rigid-sphere gases, all of the relevant omega integrals needed for these solutions are analytically evaluated and, thus, results to any desired precision can be obtained. The values of the transport coefficients obtained using Sonine polynomial expansions for the Chapman–Enskog solutions converge and, therefore, the exact diffusion and thermal conductivity solutions to a given degree of convergence can be determined with certainty by expanding to sufficiently high an order. We have used Mathematica® for its versatility in permitting both symbolic and high-precision computations. Our results also establish confidence in the results reported recently by other authors who used direct numerical techniques to solve the relevant Chapman–Enskog equations. While in all of the direct numerical methods more-or-less full calculations need to be carried out with each variation in molecular parameters, our work has utilized explicit, general expressions for the necessary matrix elements that retain the complete parametric dependence of the problem and, thus, only a matrix inversion at the final step is needed as a parameter is varied. This work also indicates how similar results may be obtained for more realistic intermolecular potential models and how other gas-mixture problems may also be addressed with some additional effort.     

Study of the aerodynamic stability of a blunt cone with a side flap in a hypersonic flow

The stationary and time-dependent aerodynamic coefficients of a slender blunt cone with a flap located near the base section of the model are experimentally investigated. The freestream parameters (M_∞ = 6, Re _L = 0.88 × 10⁷, and γ = 1.4) ensured a turbulent regime of flow over the conical surface and the flap. At high angles of attack (α ～ 10°) laminar-turbulent transition is observable in the separation zone on the leeward side of the body. Emphasis is placed on the determination of the trimming angles of attack for different positions of the center of rotation and the static and dynamic stability coefficients (the model oscillation damping coefficient).     

Measurement of the concentration diffusion coefficient for He-Ar and Ne-Kr by a two bulb method

Summary A two bulb glass apparatus was used to measure the concentration diffusion coefficient of the binary gas systems He-Ar and Ne-Kr. The coefficients were determined for equimolar mixtures at temperatures between 0°C and 70°C. The diffusion was followed as a function of time by withdrawing samples and analyzing them in a specially designed thermal conductivity analyzer with high accuracy. The diffusion coefficients agree with earlier reported experimental values and with those obtained on the basis of the Chapman-Enskog theory in conjunction with the modified Buckingham exp-six and Lennard-Jones (12-6) intermolecular potentials. The smoothed values were used to predict viscosity and thermal conductivity of these mixtures as a function of composition and temperature.     

The uniqueness of the solution of boundary-value problems of a thermal model of discharge

The paper is devoted to a nonlinear analysis of superheating [1, 2] instability of an electric discharge stabilized by electrodes [3] in the framework of a thermal model [4] where the stability of the discharge relative to the long-wave and short-wave perturbations is proved in a linear approximation. Similar boundary-value problems arise in the theories of chemically and biologically reacting mixtures [5–7], thermal breakdown of dielectrics [8], thermal explosion [9], in the investigation of nonlinear waves in semiconductors and superconductors [10, 11], and in the investigation of Couette flow with variable viscosity [12]. The uniqueness of the one-dimensional steady solutions of the thermal model of discharge and the stability relative to the small spatial perturbations, respectively, for the exponential and step dependence of the electrical conductivity on the temperature are proved in [3, 13]. The uniqueness of the solutions in the one-dimensional case for the same electrode temperature and arbitrary dependences of the electrical and thermal conductivity on the temperature is established in paper [14]. In the present paper, the existence and uniqueness of steady solutions of the thermal model of discharge in a three-dimensional formulation for arbitrary fairly smooth electrical and thermal conductivity functions of the temperature in the case of isothermal isopotential electrodes are proved analytically.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 140–145, January–February, 1986.The author expresses his gratitude to A. G. Kulikovskii and A. A. Barmin for the formulation of the problem and their discussions.     

Chapman–Enskog solutions to arbitrary order in Sonine polynomials II: Viscosity in a binary,rigid-sphere,gas mixture

The Chapman–Enskog solutions of the Boltzmann equations provide a basis for the computation of important transport coefficients for both simple gases and gas mixtures. These coefficients include the viscosity, the thermal conductivity, and the diffusion coefficient. In a preceding paper (I), for simple, rigid-sphere gases (i.e. single-component, monatomic gases) we have shown that the use of higher-order Sonine polynomial expansions enables one to obtain results of arbitrary precision that are error free. It is our purpose in this paper to report the results of our investigation of relatively high-order, standard, Sonine polynomial expansions for the viscosity-related Chapman–Enskog solutions for binary gas mixtures of rigid-sphere molecules. We note that in this work we have retained the full dependence of the solution on the molecular masses, the molecular sizes, the mole fractions, and the intermolecular potential model via the omega integrals. For rigid-sphere gases, all of the relevant omega integrals needed for these solutions are analytically evaluated and, thus, results to any desired precision can be obtained. The values of viscosity obtained using Sonine polynomial expansions for the Chapman–Enskog solutions converge monotonically from below and, therefore, the exact viscosity solution to a given degree of convergence can be determined with certainty by expanding to sufficiently high an order. We have used Mathematica® for its versatility in permitting both symbolic and high precision computations. Our results also establish confidence in the results reported recently by other authors who used direct numerical techniques to solve the relevant Chapman–Enskog equations. While in all of the direct numerical methods more-or-less full calculations need to be carried out with each variation in molecular parameters, our work utilizes explicit, general expressions for the necessary matrix elements that retain the complete parametric dependence of the problem and, thus, only a matrix inversion at the final step is needed as a parameter is varied. This work also indicates how similar results may be obtained for more realistic intermolecular potential models and how other gas-mixture problems may also be addressed with some additional effort.     

Blowing coefficient with arbitrary variation of the mass flow rate of the blown gas past a body around which flows a stream with an exponential change in the velocity

Formulas are derived permitting calculation of the linear corrections to the friction and heat-transfer coefficients with the blowing into the boundary layer of different gases, in small amounts but with a mass flow rate varying arbitrarily along the body. The case of a Mach number equal to zero and a temperature factor equal to unity was studied. Here it is postulated that bringing the relative heat-transfer coefficient down to a dependence on the dimensionless blowing renders possible, as with blowing which permits a self-similar solution, the use of the results obtained for arbitrary values of these parameters [1]. The proposed method of solution is based on the application, in the linear approximation, of a Duhamel integral for an arbitrary law of change in the mass flow rate along the body, if a solution is known with a discontinuous change in the mass flow rate. For a discontinuous change in the mass flow rate, the solution is sought using a Laplace transform; in this sense, the proposed method is similar to the method of [5].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 55–63, July–August, 1970.     

Chapman–Enskog solutions to arbitrary order in Sonine polynomials V: Summational expressions for the viscosity-related bracket integrals

The Chapman–Enskog solutions of the Boltzmann equations provide a basis for the computation of important transport coefficients for both simple gases and gas mixtures. These coefficients include the viscosity, the thermal conductivity, and the diffusion coefficient. In a preceding paper on simple gases (I), we have shown that the use of higher-order Sonine polynomial expansions enables one to obtain results of arbitrary precision that are free of numerical error. In two subsequent papers (II–III), we extended our original simple gas work to encompass binary gas mixture computations of the viscosity, thermal conductivity, diffusion, and thermal diffusion coefficients to high-order. In a fourth paper (IV) we derived general summational representations for the diffusion- and thermal conductivity-related bracket integrals and provided compact, explicit expressions for all of these bracket integrals needed to compute the diffusion- and thermal conductivity-related transport coefficients up to order 5 in the Sonine polynomial expansions used. In all of this previous work we retained the full dependence of our solutions on the molecular masses, the molecular sizes, the mole fractions, and the intermolecular potential model via the omega integrals up to the final point of solution via matrix inversion. The elements of the matrices to be inverted are, in each case, determined by appropriate combinations of bracket integrals which contain, in general form, all of the various dependencies. Since accurate expressions for the needed bracket integrals have not previously been available in the literature beyond orders 2 or 3, and since such expressions are necessary for any extensive program of computations of the transport coefficients involving Sonine polynomial expansions to higher orders, we have investigated alternative methods of constructing appropriately general bracket integral expressions that do not rely on the term-by-term, expansion and pattern matching techniques that we developed for our previous work. It is our purpose in this paper to report the results of our efforts to obtain useful, alternative, general expressions for the bracket integrals associated with the viscosity-related Chapman–Enskog solutions for gas mixtures. Specifically, we have obtained such expressions in summational form that are conducive to use in high-order viscosity coefficient computations for arbitrary gas mixtures and have computed and reported explicit expressions for all of the orders up to 5.     

Chapman–Enskog solutions to arbitrary order in Sonine polynomials IV: Summational expressions for the diffusion- and thermal conductivity-related bracket integrals

The Chapman–Enskog solutions of the Boltzmann equations provide a basis for the computation of important transport coefficients for both simple gases and gas mixtures. These coefficients include the viscosity, the thermal conductivity, and the diffusion coefficient. In a preceding paper on simple gases, we have shown that the use of higher-order Sonine polynomial expansions enables one to obtain results of arbitrary precision that are free of numerical error. In two subsequent papers, we have extended our original simple gas work to encompass binary gas mixture computations of the viscosity, thermal conductivity, diffusion, and thermal diffusion coefficients to high-order. In all of this previous work we retained the full dependence of our solutions on the molecular masses, the molecular sizes, the mole fractions, and the intermolecular potential model via the omega integrals up to the final point of solution via matrix inversion. The elements of the matrices to be inverted are, in each case, determined by appropriate combinations of bracket integrals which contain, in general form, all of the various dependencies. Since accurate, explicit, general expressions for bracket integrals are not available in the literature beyond order 3, and since such general expressions are necessary for any extensive program of computations of the transport coefficients involving Sonine polynomial expansions to higher orders, we have investigated alternative methods of constructing appropriately general bracket integral expressions that do not rely on the term-by-term, expansion and pattern matching techniques that we developed for our previous work. It is our purpose in this paper to report the results of our efforts to obtain useful, alternative, general expressions for the bracket integrals associated with the diffusion- and thermal conductivity-related Chapman–Enskog solutions for gas mixtures. Specifically, we have obtained such expressions in summational form that are conducive to use in high-order transport coefficient computations for arbitrary gas mixtures and have computed and reported explicit expressions for all of the orders up to 5.     

Extended thermodynamics of molecular ideal gases

The objective of extended thermodynamics of molecular ideal gases is the determination of the 17 fields ofmass density, velocity, energy density, pressure deviator, heat flux, intrinsic energy density and intrinsic heat flux. The intrinsic energy represents the rotational or the vibrational energy of the molecules. The necessary field equations are based upon balance laws and the system of equations is closed by constitutive relations which are characteristic for the gas under consideration. The generality of the constitutive relations is restricted by theprinciple of material frame indifference, and by the entropy principle. These principles reduce the constitutive coefficients of all fluxes to the thermal and caloric equation of state of the gas and provide inequalities for the transport coefficients. The transport coefficients can be related to the shear viscosity, the heat conductivity, and the coefficients of self-diffusion and attenuation of sound waves, so that the field equations become quite specific. The theory is in perfect agreement with the kinetic theory of molecular gases. It is shown that in non-equilibrium the temperature is discontinuous at thermometric walls. The dynamic pressure and the volume viscosity, are discussed and it is shown how extended thermodynamics and ordinary thermodynamics are related.     

Transport phenomena in dissociating gases

On the basis of the Chapman—Enskog method, a mathematical formalism of kinetic theory is developed for investigating transport phenomena in rarefied gases with allowance for an equilibrium dissociation reaction. The contribution of the dissociation collision integrals to the transport coefficients is estimated. The results of the estimate show that the presence of dissociation collisions in the gas does not have an appreciable influence on the coefficients of shear viscosity, thermal conductivity, and diffusion. On the other hand, the influence of the equilibrium dissociation reaction on the coefficient of bulk viscosity must be appreciable.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 168–173, July–August, 1981.     

Effect of inelastic collisions on the first-order slip coefficients for a diatomic gas with rotational degrees of freedom

When solving problems of inhomogeneous gas dynamics in the slip regime, it is necessary to know the boundary conditions for the velocity, temperature, heat fluxes, etc., that is, the boundary conditions for the gas macroparameters. In particular, such problems arise in developing the theory of thermophoresis of moderately large aerosol particles [1].The problem of monatomic and molecular (di- and polyatomic) gas slip along a boundary surface is considered in many publications (see, for example, [2–8]). The first-order effects include the isothermal and thermal gas slips characterized by the coefficients C_m and K_TS, respectively.In contrast to a monatomic gas, the molecules of diatomic and polyatomic gases have internal degrees of freedom, which considerably complicates the kinetic equation [9]. The lack of reliable models for the intermolecular interaction potential predetermines the need to construct model kinetic equations [10].In this study, for a diatomic gas whose molecules have rotational degrees of freedom, we propose a model kinetic equation obtained by developing the approach described in [6]. With the use of this model equation, the problem of diatomic gas slip along a plane surface is solved. As a result, for diatomic gases the coefficients C_m and K_TS, which depend on the thermophysical gas parameters and the intensity of inelastic collisions, are obtained.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, 2004, pp. 176–182. Original Russian Text Copyright © 2004 by Poddoskin.     

Experimental investigations and theoretical determination of thermal conductivity and viscosity of Al2O3/water nanofluid

Experimental investigations and theoretical determination of effective thermal conductivity and viscosity of Al₂O₃/H₂O nanofluid are reported in this paper. The nanofluid was prepared by synthesizing Al₂O₃ nanoparticles using microwave assisted chemical precipitation method, and then dispersing them in distilled water using a sonicator. Al₂O₃/water nanofluid with a nominal diameter of 43 nm at different volume concentrations (0.33–5%) at room temperature were used for the investigation. The thermal conductivity and viscosity of nanofluids are measured and it is found that the viscosity increase is substantially higher than the increase in thermal conductivity. Both the thermal conductivity and viscosity of nanofluids increase with the nanoparticle volume concentration. Theoretical models are developed to predict thermal conductivity and viscosity of nanofluids without resorting to the well established Maxwell and Einstein models, respectively. The proposed models show reasonably good agreement with our experimental results.     

Measurement of the concentration diffusion coefficient for Ne-Ar, Ne-Xe, Ne-H2, Xe-H2, H2-N2 and H2-O2 gas systems

The concentration diffusion coefficient, D ₁₂, is measured for the equimolar mixtures of Ne-Ar, Ne-Xe, Ne-H₂, Xe-H₂, H₂-N₂ and H₂-O₂ binary gas systems in a two-bulb metal apparatus in the temperature range 0 C to 100 C. These values are compared with the existing data on these systems and with the predictions of the kinetic theory in conjunction with the modified Buckingham exp-six potential. Unlike the thermal diffusion coefficient, with the simple theory it is possible to predict D ₁₂ within a few percent even for systems involving polyatomic gases. The smoothed experimental D ₁₂ values are also used to obtain data for the coefficients of viscosity and thermal conductivity at round temperatures and compositions for these systems.Nomenclature C ₂ ^t relative amount of a gas in the mixture in the bulb 2 at an instant t - C ₂ relative amount of the same gas in the mixture in the bulb 2 at equilibrium - D ₁₂ diffusion coefficient - X ₁ mole-fraction of the heavier component in the mixture - _mix viscosity coefficient - _mix thermal conductivity coefficient     

Sound absorption near the boundary dividing two liquids

It is well known that sound absorption in finite media is caused mainly by fluid viscosity and thermal conductivity. Kirchhoff [1] developed a general theory describing the mechanism of such absorption and applied it to the particular case of sound propagating in tubes. Rayleigh [2] used Kirchhoff's theory to study sound absorption by a porous wall with normal incidence of the sound wave. Konstantinov [3] also used Kirchhoff's theory to solve the problem of sound absorption by a rigid, isothermal (with infinite thermal conductivity) and a thermally insulating plane wall with arbitrary angle of sound-wave incidence. A natural extension of these efforts is a study of sound absorption on the boundary dividing two liquids. Aside from its scientific interest, such a problem is of practical significance, for example, in hydroacoustics or in creating methods for visualization of sound in gases and liquids [4]. The present study will attempt to solve this problem. The results can be applied to both liquid and solid (resinlike) materials.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 6–9, January–February, 1984.The author thanks T. P. Zhizhina for much assistance in the study.     

A modified STRUCT model for efficient engineering computations of turbulent flows in hydro-energy machinery

A modified STRUCT (MST) turbulence model for efficient engineering computations of turbulent flows in hydro-energy machinery is proposed in this paper. The MST model switches between URANS and LES-like modes using a new damping function to adjust the turbulent viscosity. Compared with the original STRUCT method, the modifications are as follows: (1) the BSL k-ω model with the Spalart-Shur correction is chosen as the new baseline to improve the sensitivity to rotation and curvature; (2) a new adaptive time-scale ratio is proposed to avoid the arbitrariness of geometric averaging operation in the original method; (3) the normalized helicity is introduced into the new damping function to detect the energy backscatter phenomenon. Five classical high Reynolds number flow cases are tested. The results show that the turbulent viscosity of the MST model is reasonably reduced in the massively separated regions and LES-like mode is activated, which captures more turbulent vortices and fluctuations on the URANS grids. With high efficiency and robustness, the MST model inherits the advantages of the original STRUCT method and improves the prediction accuracy of the turbulence with rotation and curvature, which enables efficient engineering computations of turbulent flows in hydro-energy machinery.     

Effect of rotation on ferromagnetic porous convection with a thermal non-equilibrium model

The effect of rotation on the onset of thermal convection in a horizontal layer of ferrofluid saturated Brinkman porous medium is investigated in the presence of a uniform vertical magnetic field using a local thermal non-equilibrium (LTNE) model. A two-field model for temperature representing the solid and fluid phases separately is used for energy equation. The condition for the occurrence of stationary and oscillatory convection is obtained analytically. The stability of the system has been analyzed when the magnetic and buoyancy forces are acting together as well as in isolation and the similarities as well as differences between the two are highlighted. In contrast to the non-rotating case, it is shown that decrease in the Darcy number Da and an increase in the ratio of effective viscosity to fluid viscosity Λ is to hasten the onset of stationary convection at high rotation rates and a coupling between these two parameters is identified in destabilizing the system. Asymptotic solutions for both small and large values of scaled interphase heat transfer coefficient H _t are presented and compared with those computed numerically. Besides, the influence of magnetic parameters and also parameters representing LTNE on the stability of the system is discussed and the veracity of LTNE model over the LTE model is also analyzed.     

Methods correlating data on the heat conductivity of dissociating gases

This article describes three methods of correlating experimental data on the heat conductivity of dissociating systems with a single chemical reaction. Criterial equations for the coefficients of thermal conductivity are given.In studies on dissociating gases, the problem of determining the physical constants of the transport phenomena and, in particular, the coefficient of thermal conductivity is of special interest. In such investigations, however, experimental methods in the high-temperature range are complicated and need further development, while in the case of polyatomic mixtures in which chemical reactions are taking place the practical computations are cumbersome and inaccurate. Therefore, the development of general methods of correlating the available experimental data on the heat conductivity of chemically reacting gases is of definite importance. In this article three methods of evaluating experimental data are discussed in relation to dissociating gases.     

Modeling rotation and curvature effects within scalar eddy viscosity model framework

Two approaches to incorporate the effects of rotation and curvature in scalar eddy viscosity models are explored. One is the “Modified coefficients approach” – to parameterize the model coefficients such that the growth rate of turbulent kinetic energy is suppressed or enhanced. The other is the “Bifurcation approach” – to parameterize the eddy viscosity coefficient such that the equilibrium solution bifurcates from healthy to decaying solution branches. Simple, yet, predictive models in each of these two approaches are proposed and validated on some benchmark test cases characterized by profound effects of system rotation and/or streamline curvature. The results obtained with both the models are encouraging.