Investigation of Entry Region Turbulent Flows of Non-Gray Gases

The entrance region turbulent flow between parallel planes has been investigated numerically. Variation of fluid properties with temperature has been considered. Monte Carlo simulation has been used with a narrow band non-gray model to investigate radiative transfer. The effects of radiative transfer and variable properties have been separately and systematically considered. Representative results are included for carbon monoxide and carbon dioxide. For the conditions considered, the results show that the variable properties exert a stronger influence compared to radiative transfer alone.     

Monte Carlo method of radiative transfer applied to a turbulent flame modeling with LES

Radiative transfer plays an important role in the numerical simulation of turbulent combustion. However, for the reason that combustion and radiation are characterized by different time scales and different spatial and chemical treatments, the radiation effect is often neglected or roughly modelled. The coupling of a large eddy simulation combustion solver and a radiation solver through a dedicated language, CORBA, is investigated. Two formulations of Monte Carlo method (Forward Method and Emission Reciprocity Method) employed to resolve RTE have been compared in a one-dimensional flame test case using three-dimensional calculation grids with absorbing and emitting media in order to validate the Monte Carlo radiative solver and to choose the most efficient model for coupling. Then the results obtained using two different RTE solvers (Reciprocity Monte Carlo method and Discrete Ordinate Method) applied on a three-dimensional flame holder set-up with a correlated-k distribution model describing the real gas medium spectral radiative properties are compared not only in terms of the physical behavior of the flame, but also in computational performance (storage requirement, CPU time and parallelization efficiency). To cite this article: J. Zhang et al., C. R. Mecanique 337 (2009).     

Radiative transfer effects of an axisymmetric particulate jet into a cylindrical pipe

This paper discusses the radiative transfer effects of an axisymmetric gas-particles jet into a cylindrical pipe. The medium is gray and it participates to radiation by emission, absorption and scattering. The two-phase flow problem is solved numerically by the finite volume method. We investigate the radiative transfer through a sensitivity analysis which considers the effects of the particle radiative properties and the particle number density on the temperature field and on the radiative heat fluxes of the two-phase flow domain. Analysis of the temperature profile in the cylinder, without and then with particle radiation effects, shows a decrease in the medium temperatures and thus an important role of the radiative transfer. These results also show that the presence of scattering makes the medium temperature more uniform. Finally, analysis of the particle number density, through the variation of the injection velocity, shows that a decrease in the injection velocity decreases the temperatures of the gas and particles and rapidly equilibrates the gas and particle temperatures.     

Coupled conductive radiative transfer in CO2 with narrow band statistical model and Monte-Carlo simulation

This paper describes an approach in order to determine the radiative source term and the temperature fields in a coupled conductive and radiative transfer. We consider an emitting and absorbing gas (carbon dioxide) enclosed between two horizontal and parallel reflecting plates. The medium is in local thermodynamic equilibrium (LTE). Only radiative and conductive exchanges are considered. An implicit finite difference technique is used to solve the energy conservation equation, and the narrow band statistical model (NBSM) with the Monte Carlo method (MCM) are used for simulation of radiative transfer equation in a coupled manner. This study of coupled conductive-radiative shows many phenomenons: The interaction radiation conduction in an emitting and absorbing gas is sensitive to the molar fraction of gas and the emissivities of walls. Received on 6 July 1999     

Soot and radiation modeling in laminar ethylene flames with tabulated detailed chemistry

A strategy has been developed in order to compute unsteady convective and radiative heat transfers in an industrial combustion device. This strategy involves a tabulation method to describe gas-phase chemistry, coupled with a semi-empirical soot model. A Monte Carlo method is used to evaluate gas and soot radiative transfer. This paper presents the first validation step of this strategy, in which four laminar premixed ethylene flames have been simulated. The tabulation method well predicts gas-phase species concentrations, including acetylene, considered as the main soot precursor. The soot model gives results in the experimental uncertainty range of measurements, whereas radiative powers highlight the dominating role of soot particles.     

Continuum and kinetic approaches to the simulation of the hypersonic flow past a flat plate

The data of a mathematical simulation of hypersonic flow past a flat plate at zero incidence obtained on the basis of a numerical solution of the complete Navier-Stokes equations and the Monte Carlo statistical modeling method are compared. The effect of the slip and temperature jump conditions imposed on the body surface is examined for various values of the temperature factor. The behavior of the gasdynamic variables on the body surface and in the flowfield is analyzed. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 133–145, January–February, 1997.     

Spacecraft aerodynamics and trajectory simulation during aerobraking

This paper uses a direct simulation Monte Carlo （DSMC） approach to simulate rarefied aerodynamic characteristics during the aerobraking process of the NASA Mars Global Surveyor （MGS） spacecraft. The research focuses on the flowfield and aerodynamic characteristics distribution under various free stream densities. The vari- ation regularity of aerodynamic coefficients is analyzed. The paper also develops an aerodynamics-aeroheating-trajectory integrative simulation model to preliminarily calculate the aerobraking orbit transfer by combining the DSMC technique and the classical kinematics theory. The results show that the effect of the planetary atmospheric density, the spacecraft yaw, and the pitch attitudes on the spacecraft aerodynamics is significant. The numerical results are in good agreement with the existing results reported in the literature. The aerodynamics-aeroheating-trajectory integrative simulation model can simulate the orbit transfer in the complete aerobraking mission. The current results of the spacecraft trajectory show that the aerobraking maneuvers have good performance of attitude control.     

Asymptotic solutions of the thin viscous shock layer equations near the symmetry plane of blunt bodies in hypersonic rarefied gas flow

The hypersonic rarefied gas flow over blunt bodies near the symmetry plane is investigated for the regime transitional from continuum to free-molecular. Three rarefied gas flow regimes are considered depending on the relationship between the determining parameters of the problem. For all regimes, at small Reynolds numbers, asymptotic solutions of the thin viscous shock layer equations near the symmetry plane of blunt bodies are obtained in the form of simple analytical expressions for the heat transfer, skin friction and pressure coefficients as functions of the gas-dynamic parameters of the free-stream flow and the geometric parameters and temperature of the body. With decrease in the Reynolds number these coefficients approach their values in free-molecular flow (with the accommodation coefficient equal to unity). From comparison with the data calculated using the direct simulation Monte Carlo method, the accuracy and applicability limits of the analytical solution are estimated.     

Asymptotic Investigation of the Thin Viscous Shock Layer Equations in the Neighborhood of the Stagnation Point for Low Reynolds Numbers

Hypersonic three-dimensional viscous rarefied gas flow past blunt bodies is investigated in the neighborhood of the stagnation point. The problem of applicability of the model of a thin viscous shock layer to the regime of transition from continuum to free-molecular flow is considered. In [1], it was shown that at low Reynolds numbers three hypersonic flow regimes can be distinguished and one of those regimes was investigated. In the present study an asymptotic solution of the thin viscous shock layer equations is obtained for another flow regime. With decrease in the Reynolds number the heat transfer coefficient determined by the solution obtained approaches its free-molecular value and the friction coefficient approaches its free-molecular limit, provided that the shock layer thickness is small. The analytical solution is compared with a numerical solution and the results of calculations based on direct Monte Carlo simulation.     

Relaxation processes in hypersonic rarefied binary gas mixture flow around a cylinder

Using the direct simulation Monte Carlo method, the hypersonic flow of a binary gas mixture around a cylinder is investigated over a wide rarefaction range: from an almost continuum regime (at the Knudsen number Kn = 0.01) to free-molecular flow. The effect of a small admixture of heavy diatomic particles in a light gas flow on the relaxation processes near the cylinder and the heat flux is studied.     

Multiphysics Simulation Combining Large-Eddy Simulation,Wall Heat Conduction and Radiative Energy Transfer to Predict Wall Temperature Induced by a Confined Premixed Swirling Flame

A multi-physics simulation combining large-eddy simulation, conjugate heat transfer and radiative heat transfer is used to predict the wall temperature field of a confined premixed swirling flame operating under atmospheric pressure. The combustion model accounts for the effect of enthalpy defect on the flame structure whose stabilization is here sensitive to the wall heat losses. The conjugate heat transfer is accounted for by solving the heat conduction within the combustor walls and with the Hybrid-Cell Neumann-Dirichlet coupling method, enabling to dynamically adapt the coupling period. The latter coupling procedure is enhanced to determine statistics (mean, RMS, \(\ldots \)) in a permanent regime accurately and efficiently thanks to an acceleration technique which is derived and validated. The exact radiative heat transfer equation is solved with an advanced Monte Carlo method with a local control of the statistical error. The coupled simulation is carried out with or without accounting for radiation. Excellent results for the wall temperature are achieved by the fully coupled simulation which are then further analyzed in terms of radiative effects, global energy budget and fluctuations of wall heat flux and temperature.     

Application of Model Kinetic Equations to Calculations of Super- and Hypersonic Molecular Gas Flows

For the purpose of taking the internal degrees of freedom into account, threetemperature approximating model equations, which are a generalization of the R- and ES–BGKmodels, are proposed for a diatomic gas. The surface pressure, friction, and heat transfer coefficients are compared with the direct simulation Monte Carlo (DSMC) solution in the problem of flow past a cylinder in the super- and hypersonic flow regimes. The dependence of the surface coefficients on the rotational collision number is analyzed.     

Monte Carlo simulation of self-avoiding lattice chains subject to oscillatory shear flow in polymer solution

 This paper has introduced a pseudo-potential in bond-fluctuation model to simulate oscillatory shear flow of multiple self-avoiding chains in three dimensions following our previous work under simple shear flow. The oscillatory flow field was reasonably reproduced by lattice Monte Carlo simulation using this pseudo-potential neglecting hydrodynamic interaction. By sampling the configuration distribution functions, the macroscopic viscoelasticity of semi-concentrated polymer solution was determined. Both Newtonian and non-Newtonian regimes were studied. The complex modulus and dynamic viscosity exhibit a reasonable power relation with oscillatory frequency, which is consistent with present theories and experiments. Consequently, lattice Monte Carlo simulation has been extended to model free-draining self-avoiding multi chains subject to oscillatory shear flow and to investigate associated viscoelasticity on the molecular level. Received: 1 October 1999 Accepted: 19 October 1999     

Supersonic flow of a rarefied gas around a wire grid

The direct simulation Monte Carlo method is used to study a plane-parallel supersonic gas flow through a grid formed by a number of infinite parallel wires. Characteristic features of the shock disturbance formation of in the interaction of a supersonic flow with a permeable grid are revealed. Particular attention is paid to studying the influence of geometrical parameters of the wire target on the number of particles colliding with the surface of the wires.     

Scattering of nonspherical particles rebounding from a smooth and a rough surface in a high-speed gas-particle flow

Scattering indicatrices of nonspherical particles rebounding from a smooth and a rough surface are obtained by direct Monte Carlo simulations. Particles shaped as ellipsoids of revolution, rectangular prisms, and prisms with truncated vertices are considered. Surface roughness is defined as a two-dimensional profile whose scattering characteristics are close to those of real roughness induced by abrasive erosion of the surface in a high-speed gas-particle flow. Impact interaction of an individual particle with the surface is considered in a three-dimensional formulation. The scattering indicatrices of reflected particles are found to depend substantially on the particle shape in the case of rebound from a smooth surface and to be almost independent of the particle shape if the particles rebound from a rough surface. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 2, pp. 79–88, March–April, 2008.     

Rarefied planar jets into vacuum

This paper presents a fundamental gaskinetic study on high speed rarefied jets expanding into vacuum from a cluster of planar exits. Based on the corresponding exact expressions for one planar jet, this paper straightforwardly derives the combined multiple jet flowfield solutions of density and velocity components, however, for the combined temperature and pressure solutions, extra attention shall be practiced. Several direct simulation Monte Carlo simulation results are provided and they validate these analytical solutions of rarefied planar jet flows.     

Assessment of an unstructured exponential scheme discrete ordinates radiation model for non-gray media

For radiative transfer in complex geometries, Sakami and his co-workers have developed a discrete ordinates method (DOM) exponential scheme for unstructured meshes which was mainly applied to gray media. The present study investigates the application of the unstructured exponential scheme to a wider range of non-gray scenarios found in fire and combustion applications, with the goal to implement it in an in-house Computational Fluid Dynamics (CFD) code for fire simulations. The original unstructured gray exponential scheme is adapted to non-gray applications by employing a statistical narrow-band/correlated-k (SNB-CK) gas model and meshes generated using the authors’ own mesh generator. Different non-gray scenarios involving spectral gas absorption by H₂O and CO₂ are investigated and a comparative analysis is carried out between heat flux and radiative source terms predicted and literature data based on ray-tracing and Monte Carlo methods. The maximum discrepancies for total radiative heat flux do not typically exceed 5%.     

Effects of wall conduction and rarefaction on shock propagation in a micro-channel

A Direct Simulation Monte Carlo code is coupled with a continuum unsteady conduction finite volume code to study the effects of wall conduction on shock propagation in a micro-channel. In addition, the effects of rarefaction and wall accommodation are also investigated. Wall accommodation is seen to influence the shock attenuation more significantly than wall conduction. The results indicate that both the shock front Mach number and shock strength are comparatively higher in case of adiabatic wall condition as there is no loss of energy from the shocked region. At higher rarefaction, the effect of wall accommodation as well as that of wall conduction are seen to become more noticeable. The simulation results show that it is not necessary to exactly model the wall conduction heat transfer for studying shock attenuation, and the two limiting situations of isothermal and adiabatic approximations are satisfactory as the thermal wall boundary conditions.     

极高超声速稀薄气体原子辐射效应的p-DSMC方法

极高超声速流动激波层内的高温导致内能模态的激发并伴随热辐射发生, 过高的温度使得空气分子完全解离, 原子组分对辐射热的贡献将达到80%以上. 本文基于优化的原子辐射模型, 提出追踪光子?直接模拟蒙特卡罗(p-DSMC)方法, 研究了稀薄流区不同马赫数下的高超声速二维圆柱绕流的壁面辐射加热, 获得了有无激发辐射效应的壁面压力和热流以及沿驻点线变化的平动、振动和转动温度. 在不考虑激发辐射效应的情况下, 得到的壁面压力和热流与已有的模拟结果符合的非常好, 误差均在5%以内, 尤其是在驻点位置, 误差在1%以内; 获得的平动、振动以及转动温度均与文献结果符合的很好. 在相同的来流条件下, 考虑辐射效应后发现, 来流速度低于10 km/s时, 辐射加热不明显, 在驻点区域, 辐射加热占对流加热比重在7%左右; 来流速度大于10 km/s时, 在驻点区域, 辐射加热占对流加热比重将超过30%. 考虑辐射效应后, 对非平衡区的平动、转动和振动温度的最大值影响不大. 此外, 另一个重要结论是, 流场中原子的浓度是影响壁面辐射热流大小的一个重要因素.      

Interaction of a centered isentropic compression wave with the bow shock ahead of a body in a rarefied gas

Rarefied gas flow with a centered isentropic compression wave is investigated using direct Monte Carlo simulation of the solution of the Boltzmann equation. For monatomic gas flow the pattern of formation of a suspended compression shock near the geometric center of the compression wave is considered. The flow pattern is compared with the results obtained within the framework of gas dynamics. For a diatomic gas the interference of a centered compression wave with the bow shock ahead of a cylinder is investigated. The dependence of the pressure and the heat transfer to the surface on the Reynolds number and the wave center position relative to the cylinder center is analyzed. The results are compared with those of numerical simulation of the Euler and boundary-layer equations.