Numerical Analysis of the Turbulent Flow Field Applied to Thermal Spallation Drilling

Thermal spallation is a possible drilling technique which consists of using hot supersonic jets as heat source to perforate hard rocks at high rates. This work presents a numerical analysis of a typical spallation drilling configuration, by the finite volume method. The time-averaged conservation equations of mass, momentum and energy are solved to determine the turbulent compressible gas phase flow field. Turbulence is predicted by the classical high Reynolds number κ-ε model, as well as with a low Reynolds number κ-ε model. The influence of the jet Reynolds number is investigated. Special attention is given to the rock surface temperature, since its accurate determination is required to predict spallation rates under field-drilling conditions.

Stefan Number-Insensitive Numerical Simulation of the Enthalpy Method for Stefan Problems Using the Space-Time Conservation Element and Solution Element Method

An improvement is introduced to the conservation element and solution element (CE/SE) phase change scheme presented previously. The improvement addresses a well known weakness in numerical simulations of the enthalpy method when the Stefan number, (the ratio of sensible to latent heat) is small (less than 0.1). Behavior of the improved scheme, at the limit of small Stefan numbers, is studied and compared with that of the original scheme. It is shown that high dissipative errors, associated with small Stefan numbers, do not occur using the new scheme.

Vehicle Thermal Management Simulation Using a Rapid Omni-Tree Based Adaptive Cartesian Mesh Generation Methodology

CFD simulation of vehicle under-hood and under-body poses several challenges. Specifically, the complexity of the geometry involved makes the use of traditional mesh generation approaches, based on the boundary-to-interior methodology, impractical and time consuming. The current work presents the use of an interior-to-boundary method wherein the need for creating a ‘water-tight’ surface mesh is not a pre-requisite for volume mesh generation. The application of the new method is demonstrated for an actual passenger vehicle under-hood model with nearly a hundred components. Coupled radiation/convection simulations are performed to obtain the complete airflow and thermal map of the engine compartment. Results are validated with test data. The new method results in significant gains in efficiency over traditional approaches allowing the simulation tool to be used effectively in the vehicle development process.

Simulation of Flow Around a Static and Oscillating in Pitch NACA 0015 Airfoil Using URANS and DES

This paper presents numerical simulations of the flow around a NACA 0015 airfoil at static and dynamic stall. The treatment of these configurations is a very challenging task for CFD applications. The turbulent flow around the static and in pitch oscillation airfoil is computed using different approaches: 2D RANS, 3D RANS and DES methodologies and with finer and finer meshes in order to try to reach a space converged solution. The main conclusion of the paper is that the prediction of static and all the more dynamic stall is not mature with present modeling capabilities.

Rheological Aspects of Drops Deforming in Finite Reynolds Number Oscillatory Extensional Flows

The evolving morphology of droplets in a flowing emulsion determines its rheological properties. A two-way interaction between drops and the flow governs the rheological stresses arising from drop deformation. In this paper, the rheology of droplet emulsions under oscillatory extensional flow is investigated using direct numerical simulation (DNS). The deformation of a three dimensional drop is simulated. The rheological responses are related with the interface morphology using Bachelor’s stress formulation [6]. Detailed investigation of the variation of parameters such as interfacial tension, flow frequency and inertia displayed complex non-Newtonian response of the emulsion that will have broad implication in industrial applications. The results are explained and discussed with a simple model for the drop dynamics.

Simulation of Contaminant Dispersal in an Apartment Building

Bio-terrorism events (like the 2001 anthrax attacks) accentuate the importance of countering these incidents. In order to develop reliable countermeasures for these events, it is essential to understand the associated transport processes. The transport processes involved pose challenges as they occur over wide ranges of spatial and temporal scales. CONTAMW, a multi zone indoor air quality and ventilation analysis program is used to predict the contaminant dispersal in an apartment building. Detailed simulation results and analysis of controlled release of propylene within a generic apartment building is presented. A zonal analysis is carried out for the entire apartment building (using CONTAMW) to obtain time histories of propylene concentration in different zones. The simulations provide the dispersion, transport and contaminant concentration within each zone of the apartment. This study also considers the effect of flow obstructions and ventilation rates on contaminant dispersal. The results are validated with the experimental results reported in Cybyk et al. (1999). We have also simulated propylene transport in the apartment with FDS, a large eddy simulation model.

Simulation of Acoustic Streaming in a Resonator

Formation of vortical flow structures in a rectangular enclosure due to acoustic streaming is investigated numerically. The oscillatory flow field in the enclosure is created by the vibration of a vertical side wall of the enclosure. The frequency of the wall vibration is chosen such that a standing wave forms in the enclosure. The interaction of this standing wave with the horizontal solid walls leads to the production of Rayleigh type acoustic streaming flow patterns in the enclosure. All four walls of the enclosure considered are thermally insulated. The fully compressible form of the Navier-Stokes equations is considered and an explicit time-marching algorithm is used to explicitly track the acoustic waves. Numerical solutions are obtained by employing a highly accurate flux corrected transport (FCT) algorithm for the convection terms. A time-splitting technique is used to couple the viscous and diffusion terms of the full Navier-Stokes equations. Non-uniform grid structure is employed in the computations. The simulation of the primary oscillatory flow and the secondary (steady) streaming flows in the enclosure is performed. Streaming flow patterns are obtained by time averaging the primary oscillatory flow velocity distributions. The effect of the amount of wall displacement on the formation of the oscillatory flow field and the streaming structures are studied. Computations indicate that the nonlinearity of the acoustic field increases with increasing amount of the vibration amplitude. The form and the strength of the secondary flow associated with the oscillatory flow field and viscous effects are found to be strongly correlated to the maximum displacement of the vibrating wall. Total number of acoustic streaming cells per wavelength is also determined by the strength and the level of the nonlinearity of the sound field in the resonator.

Fluid Flow and Heat Transfer Around Two Circular Cylinders in Side-by-Side Arrangement

Incompressible flows passing through two circular cylinders in side-by-side arrangement are investigated numerically. The calculations are carried out with pitch ratios from 1.1 to 2.0 at Reynolds number of 1000. The flow and temperature fields, flow interference, and the local and the mean Nusselt numbers are studied in this research. It is observed that for the pitch ratios in the range of 2.0 and 1.5, the emerging jet between cylinders deflects and one wide and one narrow wakes behind the cylinders are formed. The gap flow velocity increases as the pitch ratio decreases and consequently increases the mean Nusselt number of the cylinders. As the pitch ratio decreases and is less than 1.5, the jet deflection is more severe and the gap flow velocity starts to decrease slowly, which results in reducing the mean Nusselt number of the cylinders. Due to the rapid reduction of the narrow wake size, the mean Nusselt number of the cylinder with narrow wake shows an uprising tendency for the decreasing pitch ratio less than 1.2.

LES of Turbulent Separated Flow and Heat Transfer in a Symmetric Expansion Plane Channel

LES method is applied to simulate numerically a turbulent separated and reattached flow and heat transfer in a symmetric expansion plane channel of expansion ratio 2.0. Smagorinsky model is used in the analysis and fundamental equations are discretized by means of the finite difference method, and their resulting finite difference equations are solved using SMAC method. The calculations are conducted for Re = 15000. It is found that the present numerical results, in general, agree well with the previous experimental ones. The complicated vortical flow structures in the channel and their correlations with heat transfer characteristics are visualized through various fields of flow quantities.

Three-Dimensional Simulation of Turbulent Rectangular and Square Impinging Jet Flows

Flow characteristics of turbulent impinging jets issuing, respectively, from a rectangular and a square nozzles have been investigated numerically through the solution of three-dimensional Navier-Strokes equations in steady state. Two geometries with two nozzle-to-plate spacings of four and eight times of hydraulic diameters of the jet pipes, and two Reynolds numbers of 20000 and 30000 have been considered with fully developed inlet boundary conditions. An RNG based k–ε turbulence model and a deferred correction QUICK scheme in conjunction with the wall function method have been applied to the prediction of the flow fields within semi-confined spaces. A common feature revealed by the computational results is the presence of a toroidal recirculation zone around the jet. An adverse pressure gradient is found at the impingement surface downstream the stagnation point. Boundary layer separation will occur if the gradient is strong enough, and the separation manifests itself as a secondary recirculation zone at the surface. In addition, three-dimensional simulations reveal the existence of two and four pronounced streamwise velocity off-center peaks at the cross-planes near to the impingement plate, respectively, in the rectangular and square impinging jet flows. These peaks are found forming at the horizontal planes where the wall jets start forming accompanied by two or four pairs of counter-rotating vortex rings. It is believed that the formation of the off-center velocity peaks is due to the vorticity diffusion along the wall jet as the jet impinges on the target plate.