Shock-Induced Dust Cloud Over a Deposit Layer of Fine Particles

The present paper investigates particle density pro les of a dust cloud induced by a normal shock wave moving at a constant speed along a at surface deposited with ne particles. In shock-fixxed coordinates, numerical simulation of ow structures of the carrier- and dispersed- phases was performed for the M = 2 case. The neness and non-uniformity of the particle size are taken into account and their effcts on the dust cloud are discussed in detail.

Computation of flow over a rotating body on unstructured Chimera mesh

Flow around moving boundary is ubiquitous in engineering applications. To increse the efficienly of the algorithm to handle moving boundaries is still a major challenge in Computational Fluid Dynamics (CFD). The Chimera grid method is one type of method to handle moving boundaries. A concept of domain de-composition has been proposed in this paper. In this method, sub-domains are meshed independently and governing equations are also solved separately on them. The Chimera grid method was originally used only on structured (curvilinear) meshes. However, in a problem which involves both moving boundary and complex geometry, the number of sub-domains required in a traditional (structured) Chimera method becomes fairly large. Thus the time required in the interior boundary locating, link-building and data exchanging also increases. The use of unstructured Chimera grid can reduce the time consumption significantly by the reduction of domain(block) number. Generally speaking, unstructured Chimera grid method has not been developed. In this paper, a well-known pressure correction scheme - SIMPLEC is modified and implemented on unstructured Chimera mesh. A new interpolation scheme regarding the pressure correction is proposed to prevent the possible decoupling of pressure. A moving-mesh finite volume approach is implemented in an inertial reference frame. This approach is then used to compute incompressible flow around a rotating circular and elliptic cylinder. These numerical examples demonstrate the capability of the proposed scheme in handling moving boundaries. The numerical results are in good agreement with other experimental and computational data in literature. The method proposed in this paper can be efficiently applied to more challenge cases such as free-falling objects or heavy particles in fluid.

An Atomistic-Continuum Hybrid Simulation Of Fluid Flows Over Superhydrophobic Surfaces

Recent experiments have found that slip length could be as large as on the order of 1 mu m for fluid flows over superhydrophobic surfaces. Superhydrophobic surfaces can be achieved by patterning roughness on hydrophobic surfaces. In the present paper, an atomistic-continuum hybrid approach is developed to simulate the Couette flows over superhydrophobic surfaces, in which a molecular dynamics simulation is used in a small region near the superhydrophobic surface where the continuum assumption is not valid and the Navier-Stokes equations are used in a large region for bulk flows where the continuum assumption does hold. These two descriptions are coupled using the dynamic coupling model in the overlap region to ensure momentum continuity. The hybrid simulation predicts a superhydrophobic state with large slip lengths, which cannot be obtained by molecular dynamics simulation alone.

Effects of wall temperature on boundary layer stability over a blunt cone at Mach 7.99

Effects of wall temperature on stabilities of hypersonic boundary layer over a 7-degree half-cone-angle blunt cone are studied by using both direct numerical simulation (DNS) and linear stability theory (LST) analysis. Four isothermal wall cases with Tw/T0= 0.5, 0.7, 0.8 and 0.9, as well as an adiabatic wall case are considered. Results of both DNS and LST indicate that wall temperature has significant effects on the growth of disturbance waves. Cooling the surface accelerates unstable Mack II mode waves and decelerates the first mode (Tollmien–Schlichting mode) waves. LST results show that growth rate of the most unstable Mack II mode waves for the cases of cold wall Tw/T0=0.5 and 0.7 are about 45% and 25% larger than that for the adiabatic wall, respectively. Numerical results show that surface cooling modifies the profiles of rdut/dyn and temperature in the boundary layers, and thus changes the stability haracteristic of the boundary layers, and then effects on the growth of unstable waves. The results of DNS indicate that the disturbances with the frequency range from about 119.4 to 179.1 kHz, including the most unstable Mack modes, produce strong mode competition in the downstream region from about 11 to 100 nose radii. And adiabatic wall enhances the amplitudes of disturbance according to the results of DNS, although the LST indicates that the growth rate of the disturbance of cold wall is larger. That because the growth of the disturbance does not only depend on the development of the second unstable mode.

An atomistic-continuum hybrid simulation of fluid flows over superhydrophobic surfaces

Recent experiments have found that slip length could be as large as on the order of 1 mu m for fluid flows over superhydrophobic surfaces. Superhydrophobic surfaces can be achieved by patterning roughness on hydrophobic surfaces. In the present paper an atomistic-continuum hybrid approach is developed to simulate the Couette flows over superhydrophobic surfaces in which a molecular dynamics simulation is used in a small region near the superhydrophobic surface where the continuum assumption is not valid and the Navier-Stokes equations are used in a large region for bulk flows where the continuum assumption does hold. These two descriptions are coupled using the dynamic coupling model in the overlap region to ensure momentum continuity. The hybrid simulation predicts a superhydrophobic state with large slip lengths which cannot be obtained by molecular dynamics simulation alone.

Direct numerical simulation of hypersonic boundary layer transition over a blunt cone with a small angle of attack

The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [Li et al., AIAA J. 46, 2899(2008)], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 0° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves VLFWs� are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [ Li et al., AIAA J. 46, 2899 (2008) ], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 20° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves (VLFWs) are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.

Simulation of gas flow over a micro cylinder

Gas flow over a micro cylinder is simulated using both a compressible Navier-Stokes solver and a hybrid continuum /particle approach. The micro cylinder flow has low Reynolds number because of the small length scale and the low speed, which also indicates that the rarefied gas effect exists in the flow. A cylinder having a diameter of 20 microns is simulated under several flow conditions where the Reynolds number ranges from 2 to 50 and the Mach number varies from 0.1 to 0.8. It is found that the low Reynolds number flow can be compressible even when the Mach number is less than 0.3, and the drag coefficient of the cylinder increases when the Reynolds number decreases. The compressible effect will increase the pressure drag coefficient although the friction coefficient remains nearly unchanged. The rarefied gas effect will reduce both the friction and pressure drag coefficients, and the vortex in the flow may be shrunk or even disappear.

Effect of coal-over-coal reburn on furnace temperature and heat flux distributions in 1 MW tangentially fired furnace

The furnace temperature and heat flux distributions of 1 MW tangentially fired furnace were studied during coal-over-coal reburn, and the influences of the position of reburn nozzle and reburn fuel fraction on furnace temperature and heat flux distributions were investigated. Compared with the baseline, the flue gas temperature is 70–90 C lower in primary combustion and 130–150 C higher at furnace exit, and the variations of the flue gas temperature distributions along furnace height are slower. The temperature distribution along the width of furnace wall decreases with the increase of the relative furnace height. In the primary combustion zone and the reburn zone, the temperature and heat flux distributions of furnace wall are much non-uniform and asymmetric along the width of furnace wall, those of furnace wall in the burnout zone are relatively uniform, and the temperature non-uniformity coefficients of the primary combustion zone, the reburn zone and the burnout zone are 0.290, 0.100 and 0.031, respectively.