Search results

A calculation procedure for turbulent natural convection in enclosures is described. A two‐equation model based on the eddy diffusivity concept for the temperature field possessing a form similar to the k—ε model of flow is incorporated, thus, extending the applicability of the eddy diffusivity models by removing constraints of the Reynolds analogy between momentum and thermal transport processes. As a test problem, natural convection in a square cavity subjected to differential side‐wall heating is analysed. The vertical walls are divided into isothermal and constant heat‐flux surfaces and heated non‐uniformly. AtRa = 10¹⁰ and for an air—filled cavity (Pr = 0.71), variations of heating patterns are found to significantly alter the field characteristics. Numerical predictions demonstrate dissimilar features of the velocity and temperature fluctuations.

The purpose of this study is to determine the effect of ventilation openings and fire intensity on heat transfer and fluid flow within the microclimate between 3D human body and clothing.

On account of interaction effects of fire and ventilation openings on heat transfer process, a 3D transient computational fluid dynamics model considering the real shape of human body and clothing was developed. The model was validated by comparing heat flux history and distribution with experimental results. Heat transfer modes and fluid flow were investigated under three levels of fire intensity for the microclimate with ventilation openings and closures.

Temperature distribution on skin surface with open microclimate was heavily depended on the heat transfer through ventilation openings. Higher temperature for the clothing with confined microclimate was affected by the position and direction of flames injection. The presence of openings contributed to the greater velocity at forearms, shanks and around neck, which enhanced the convective heat transfer within microclimate. Thermal radiation was the dominant heat transfer mode within the microclimate for garment with closures. On the contrary, convective heat transfer within microclimate for clothing with openings cannot be neglected.

The findings provided fundamental supports for the ease and pattern design of the improved thermal protective systems, so as to realize the optimal thermal insulation of the microclimate on the garment level in the future.

The outcomes broaden the insights of results obtained from the mesoscale models. Different high skin temperature distribution and heat transfer modes caused by thermal environment and clothing structure provide basis for advanced thermal protective clothing design.

The purpose of this study is to present and demonstrate a numerical method for solving chemically reacting flows. These are important for energy conversion devices, which rely on chemical reactions as their operational mechanism, with heat generated from the combustion of the fuel, often gases, being converted to work.

The numerical study of such flows requires the set of Navier-Stokes equations to be extended to include multiple species and the chemical reactions between them. The numerical method implemented in this study also accounts for changes in the material properties because of temperature variations and the process to handle steep spatial fronts and stiff source terms without incurring any numerical instabilities. An all-speed numerical framework is used through simple low-dissipation advection upwind splitting (SLAU) convective scheme, and it has been extended in a multi-component species framework on the in-house density-based flow solver. The capability of solving turbulent combustion is also implemented using the Eddy Dissipation Concept (EDC) framework and the recent k-kl turbulence model.

The numerical implementation has been demonstrated for several stiff problems in laminar and turbulent combustion. The laminar combustion results are compared from the corresponding results from the Cantera library, and the turbulent combustion computations are found to be consistent with the experimental results.

This paper has extended the single gas density-based framework to handle multi-component gaseous mixtures. This paper has demonstrated the capability of the numerical framework for solving non-reacting/reacting laminar and turbulent flow problems. The all-speed SLAU convective scheme has been extended in the multi-component species framework, and the turbulent model k-kl is used for turbulent combustion, which has not been done previously. While the former method provides the capability of solving for low-speed flows using the density-based method, the later is a length-scale-based method that includes scale-adaptive simulation characteristics in the turbulence modeling. The SLAU scheme has proven to work well for unsteady flows while the k-kL model works well in non-stationary turbulent flows. As both these flow features are commonly found in industrially important reacting flows, the convection scheme and the turbulence model together will enhance the numerical predictions of such flows.

This paper aims to assess the sensitivity of numerical simulation results of turbulent reactive flow to the formulation of inlet boundary conditions. The analysis concerns the profiles of the mean velocity the turbulence kinetic energy k and its dissipation rate ϵ. It is intended to provide guidance to the determination of inlet conditions when only global flow data are available. This situation can be met both in simple laboratory experiments and in industrial full-scale applications, when measurements are either incomplete or infeasible, resulting in lack of detailed inlet data.

Two turbulence–chemistry interaction models were studied: eddy dissipation concept and partially stirred reactor. Three different velocity profiles and related turbulence statistics were applied to present feasible scenarios and their consequences. Simulations with the most appropriate inlet data were accompanied with profiles of turbulent quantities obtained with a proposed method. This method was contrasted to other approaches popular in the literature: the pre-inlet pipe and the separate cold flow simulations of a burner. The methodology was validated on two laboratory-scale jet flames: Delft Jet-in-Hot-Coflow and Sandia CHN B. The simulations were carried out with open source code OpenFOAM.

The proposed relations for turbulence kinetic energy and its dissipation rate at the inlet are found to provide results comparable to those obtained with the use of experimental data as inlet boundary conditions. Moreover, from a certain location downstream the jet, weakly dependent on the Reynolds number, the influence of inlet conditions on flow statistics was found to be negligible.

This work reveals the consequences of the use of rather crude assumptions made for inlet boundary conditions. Proposed formulas for the profiles for k and epsilon are attractive alternatives to other approaches aiming to determine the inlet boundary conditions for turbulent jet flows.

This paper aims to develop a multidomain boundary element method (BEM) for modeling 2D complex turbulent thermal flow using low Reynolds two‐equation turbulence models.

The integral boundary domain equations are discretised using mixed boundary elements and a multidomain method also known as a subdomain technique. The resulting system matrix is an overdetermined, sparse block banded and solved using a fast iterative linear least squares solver.

The simulation of a turbulent flow over a backward step is in excellent agreement with the finite volume method using the same turbulent model. A grid consisting of over 100,000 elements could be solved in the order of a few minutes using a 3.0 Ghz P4 and 1 GB memory indicating good efficiency.

The paper shows, for the first time, that the BEM is applicable to thermal flows using k‐ε.

Ultrahigh-speed projectile running in water with the velocity close to the speed of sound usually causes large supercavity. The computation of such transonic cavitating flows is usually difficult, thus high-speed model reflecting the compressibility of both the liquid and the vapor phases should be introduced to model such flow. The purpose of this paper is to achieve a model within an in-house developed solver to simulate the ultrahigh-speed subsonic supercavitating flows.

An improved TAIT equation adjusted by local temperature is adopted as the equation of state (EOS) for the liquid phase, and the Peng-Robinson EOS is used for the vapor phase. An all-speed variable coupling algorithm is used to unify the computations and regulate the convergence at arbitrary Mach number. The ultrahigh-speed (Ma=0.7) supercavitating flows around circular disk are investigated in contrast with the case of low subsonic (Ma=0.007) flow.

The characteristic physical variables are reasonably predicted, and the cavity profiles are compared to be close to the experimental empirical formula. An important conclusion in the compressible cavitating flow theory is verified by the numerical result that, at any specific cavitation number the cavity’s size and the drag coefficient both increase along with the rise of Mach number. On the contrary, it is found as well that the cavity’s slenderness ratio decreases when Mach number goes up. It indicates that the compressibility has different influences on the length and the radius of the supercavity.

A high-speed model reflecting the compressibility of both the liquid and the vapor phases was suggested to model the ultrahigh-speed supercavitating flows around underwater projectiles.

The purpose of this article is to explore the changing skill sets, operating environments and community engagement activities which can create a robust and valued future for libraries and librarians. It is easy to paint a doom-laden picture of the future of libraries against a background of library cuts, professional deskilling and the relentless advance of ever more powerful digital information systems. However, these self-same threats offer new challenges for information management.

Personal view of the future based on current literature.

It is better to seek forgiveness than permission. Seize the day.

Just one point of view among many.

Librarians need to act.

The library can remain core to civilised societies.

Others must judge that.

The purpose of this paper is to propose an efficient/robust numerical algorithm for solving the two‐dimensional laminar mixed‐convection in a lid‐driven cavity using the mixed finite element (FE) technique.

A numerical algorithm was based on the so‐called consistent splitting scheme, which improved the numerical accuracy of the primitive variables. In order to obtain a stable solution, two choices of mixed FEs, the Taylor‐Hood and Crouzeix‐Raviart types, were used. Two mesh layouts were considered; uniform and non‐uniform.

To verify that the proposed scheme had a second‐order accuracy, some numerical results are presented and compared with the known solution. The answer was confirmative. Numerically accurate solutions were obtained for a fixed Prandtl number, Pr=0.71, for a range of the Reynolds number, Re from 100 to 3,000, and for a range of the Richardson number, Ri from 0.001 to 100. The results from these calculations, using the mixed FE consistent splitting scheme, agreed with the existing ones.

Further extensions of this work could include the influence of various choices of Reynolds numbers, Prandtl numbers and Richardson numbers, and the effect of aspect ratio.

The present work was the first to apply a mixed FE in association with the consistent splitting scheme to the mixed convection problem.