Development of an efficient particle approach for micro-scale gas flow simulations

The micro-scale gas flows are usually low-speed flows and exhibit rarefied gas effects. It is challenging to simulate these flows because traditional CFD method is unable to capture the rarefied gas effects and the direct simulation Monte Carlo (DSMC) method is very inefficient for low-speed flows. In this study we combine two techniques to improve the efficiency of the DSMC method. The information preservation technique is used to reduce the statistical noise and the cell-size relaxed technique is employed to increase the effective cell size. The new cell-size relaxed IP method is found capable of simulating micro-scale gas flows as shown by the 2D lid-driven cavity flows.

NUMERICAL AND EXPERIMENTAL STUDY ON LIQUID-SOLID FLOW IN A HYDROCYCLONE

Hydrocyclones are widely used in industry, of which the geometrical design using CFD techniques is gaining more popularity in recent years. In this study, the Euler-Euler approach and the Reynolds stress model are applied to simulate the liquid-solid flowfield in a hydrocyclone. The methodology is validated by a good agreement between experimental data and numerical results. Within the research range, the simulation indicates that the liquid-solid separation mainly occurs in the conical segment, and increasing conical height or decreasing cylindrical height helps to improve the grade efficiencies of solid particles. Based on these results, two of the same hydrocyclones are designed and installed in series to establish a liquid-solid separation system. Many experiments are then conducted under different conditions, in which the effects of the water cut and the second hydrocyclone on the separation are investigated. The results also confirm that smaller solid particles are more susceptible to the inlet conditions, and the second hydrocyclone plays a more important role as the water cut reduces.

Dissipative particle dynamics simulation of wettability alternation phenomena in the chemical flooding process

Wettability alternation phenomena is considered one of the most important enhanced oil recovery (EOR) mechanisms in the chemical flooding process and induced by the adsorption of surfactant on the rock surface. These phenomena are studied by a mesoscopic method named as dissipative particle dynamics (DPD). Both the alteration phenomena of water-wet to oil-wet and that of oil-wet to water-wet are simulated based on reasonable definition of interaction parameters between beads. The wetting hysteresis phenomenon and the process of oil-drops detachment from rock surfaces with different wettability are simulated by adding long-range external forces on the fluid particles. The simulation results show that, the oil drop is liable to spread on the oil-wetting surface and move in the form of liquid film flow, whereas it is likely to move as a whole on the water-wetting surface. There are the same phenomena occuring in wettability-alternated cases. The results also show that DPD method provides a feasible approach to the problems of seepage flow with physicochemical phenomena and can be used to study the mechanism of EOR of chemical flooding.

Axial compressor response to inlet flow distortions by a CFD analysis

The usual approach to compressor design considers uniform inlet flow characteristics. Especially in aircraft applications, the inlet flow is quite often non uniform, and this can result in severe performance degradation. The magnitude of this phenomenon is amplified in military engines due to the complexity of inlet duct configurations and the extreme flight conditions. CFD simulation is an innovative and powerful tool for studying inlet distortions and can bring this inside the very early phases of the design process. This project attempts to study the effects of inlet flow distortions in an axial flow compressor trying to minimize the use computer resources and computational time. The first stage of a low bypass ratio compressor has been analyzed and its clean and distorted performance compared outlining the principal changes due to uneven flow distribution: drop in mass flow, increase in pressure and temperature ratios, decrease in surge margin. Three different studies have then been conducted to better understand the effects of the level, the type and the frequency of the distortion.

Vortex filament simulation of trailing vortex merger

Aircraft in high-lift configuration shed multiple vortices. These generally merge to form a downstream wake consisting of two counter-rotating vortices of equal strength. The understanding of the merger of two co-rotating trailing vortices is important in evaluating the separation criteria for different aircraft to prevent wake vortex hazards during landing and take-off. There is no existing theoretical method on the basis of which such norms can be set. The present study is aimed at gaining a better understanding of the behaviour of wake vortices behind the aircraft. Two dimensional studies are carried out using the vortex blob method and compared with Bertenyi's experiment. It is shown that inviscid two dimensional effects are insufficient to explain the observations. Three dimensional studies, using the vortex filament method, are applied to the same test case. Two Lamb-Oseen profile vortices of the same dimensions and initial separation as the experiment are allowed to evolve from a straight starting condition until a converged steady flow is achieved. The results obtained show good agreement with the experimental distance to merger. Core radius and separation behaviour is qualitatively similar to experiment, with the exception of rapid increases at first. This may be partially attributable to the choice of filament-element length, and recommended further work includes a convergence study for this parameter. Copyright © 2005 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Reynolds number influence on turbulent blocking effects of a rigid plane boundary (3rd report, direct numerical simulation on wall blocking effects for anisotropic turbulence)

The Reynolds number influence on turbulent blocking effects by a rigid plane boundary is studied using direct numerical simulation (DNS). A new forcing method proposed in the second report using Townsend's "simple model eddies" for DNS was extended to generate axisymmetric anisotropic turbulence. A force field is obtained in real space by sprinkling many space-filling "simple model eddies" whose centers are randomly but uniformly distributed in space. The axes of rotation are controlled in this study to generate axisymmetric anisotropic turbulence. The method is applied to a shear-free turbulent boundary layer over a rigid plane boundary and the blocking effects for anisotropic turbulence are investigated. The results show that stationary axisymmetric anisotropic turbulence is generated using the present method. Turbulence intensities near the wall showed good agreements with the rapid distortion theory (RDT) for small t (t ≪ TL), where TL. is the eddy turnover time. The splat effect (i. e. turbulence intensities of the components parallel to the surface are amplified) occurs near the boundary and the viscous effect attenuates the splat effect at the quasi steady state at low Reynolds number as for Isotropic turbulence. Prandtl's secondary flow of the second kind does not occur for low Reynolds number flows, which qualitatively agrees with previous observetion in a mixing-box.

Numerical simulation of the growth of ZnO nanostructures in a tube furnace by physical vapour deposition

Zinc oxide is a versatile II-VI naturally n-type semiconductor that exhibits piezoelectric properties. By controlling the growth kinetics during a simple carbothermal reduction process a wide range of 1D nanostructures such as nanowires, nanobelts, and nanotetrapods have been synthesized. The driving force: for the nanostructure growth is the Zn vapour supersaturation and supply rate which, if known, can be used to predict and explain the type of crystal structure that results. A model which attempts to determine the Zn vapour concentration as a function of position in the growth furnace is described. A numerical simulation package, COMSOL, was used to simultaneously model the effects of fluid flow, diffusion and heat transfer in a tube furnace made specifically for ZnO nanostructure growth. Parameters such as the temperature, pressure, and flow rate are used as inputs to the model to show the effect that each one has on the Zn concentration profile. An experimental parametric study of ZnO nanostructure growth was also conducted and compared to the model predictions for the Zn concentration in the tube. © 2008 Materials Research Society.

The time-resolved flow within an unsteady ejector

The flow field within an unsteady ejector has been investigated using experimental and computational techniques. The experimental results show a peak thrust augmentation of 1.4; numerical simulation gives a value of 1.37. It is shown that the vortex ring dominates the flow field. At optimal thrust augmentation the vortex ring acts like a fluid piston accelerating the fluid inside the ejector. A model is proposed for the operation of unsteady ejectors, based on the vortex ring acting like a fluid piston. Control volume analysis is presented showing that mass entrainment is responsible for thrust augmentation. It is proposed that the spacing of successive vortex rings determines the mass entrainment and therefore thrust augmentation. The efficiency of unsteady ejectors was found to vary between 28% and 32% depending on the L/D ratio of the unsteady jet source. Copyright © 2008 by J H Heffer.

Geometrical properties and turbulent flame speed measurements in stationary premixed V-flames using direct numerical simulation

Three dimensional, fully compressible direct numerical simulations (DNS) of premixed turbulent flames are carried out in a V-flame configuration. The governing equations and the numerical implementation are described in detail, including modifications made to the Navier-Stokes Characteristic Boundary Conditions (NSCBC) to accommodate the steep transverse velocity and composition gradients generated when the flame crosses the boundary. Three cases, at turbulence intensities, u′/sL, of 1, 2, and 6 are considered. The influence of the flame holder on downstream flame properties is assessed through the distributions of the surface-conditioned displacement speed, curvature and tangential strain rates, and compared to data from similarly processed planar flames. The distributions are found to be indistinguishable from planar flames for distances greater than about 17δth downstream of the flame holder, where δth is the laminar flame thermal thickness. Favre mean fields are constructed, and the growth of the mean flame brush is found to be well described by simple Taylor type diffusion. The turbulent flame speed, sT is evaluated from an expression describing the propagation speed of an isosurface of the mean reaction progress variable c̃ in terms of the imbalance between the mean reactive, diffusive, and turbulent fluxes within the flame brush. The results are compared to the consumption speed, sC, calculated from the integral of the mean reaction rate, and to the predictions of a recently developed flame speed model (Kolla et al., Combust Sci Technol 181(3):518-535, 2009). The model predictions are improved in all cases by including the effects of mean molecular diffusion, and the overall agreement is good for the higher turbulence intensity cases once the tangential convective flux of c̃ is taken into account. © 2010 Springer Science+Business Media B.V.

Computational fluid-structure interaction of DGB parachutes in compressible fluid flow

This paper describes large-scale simulations of compressible flows over a supersonic disk-gap-band parachute system. An adaptive mesh refinement method is used to resolve the coupled fluid-structure model. The fluid model employs large-eddy simulation to describe the turbulent wakes appearing upstream and downstream of the parachute canopy and the structural model employed a thin-shell finite element solver that allows large canopy deformations by using subdivision finite elements. The fluid-structure interaction is described by a variant of the Ghost-Fluid method. The simulation was carried out at Mach number 1.96 where strong nonlinear coupling between the system of bow shocks, turbulent wake and canopy is observed. It was found that the canopy oscillations were characterized by a breathing type motion due to the strong interaction of the turbulent wake and bow shock upstream of the flexible canopy. Copyright © 2010 by ASME.