A computer code for enclosed natural convection in the boundary layer regime

A computer code is developed for the numerical prediction of natural convection in rectangular two-dimensional cavities at high Rayleigh numbers. The governing equations are retained in the primitive variable form. The numerical method is based on finite differences and an ADI scheme. Convective terms may be approximated with either central or hybrid differencing for greater stability. A non-uniform grid distribution is possible for greater efficiency. The pressure is dealt with via a SIMPLE type algorithm and the use of a fast elliptic solver for the solenoidal velocity correction field significantly reduces computing times. Preliminary results indicate that the code is reasonably accurate, robust and fast compared with existing benchmarks and finite difference based codes, particularly at high Rayleigh numbers. Extension to three-dimensional problems and turbulence studies in similar geometries is readily possible and indicated.

The laminar boundary layer regime for natural convection of air in a square cavity

Numerical predictions are obtained for laminar natural convection of air in a square two dimensional cavity at high Rayleigh numbers. Proper resolution of the core reveals weak multi-cellular structure which varies in a complex manner as the effects of convection are increased. The end of the steady laminar regime is numerically estimated to occur at Ra=2.2x10^8.

Laminar and turbulent natural convection in rectangular cavities

A computer program has been developed for the prediction of buoyancy-driven laminar and turbulent flow in rectangular air-filled two-dimensional cavities with differentially heated side walls. Laminar flow predictions for a square cavity and Rayleigh numbers from Ra = 10^3 up to the onset of unsteady flow have been obtained. Accurate solutions for Ra = 5 x 10^6, 10^7, 5 x 10^7 and 10^8 are presented and an estimate for the critical Rayleigh number at which the steady laminar flow becomes unsteady is given for this geometry. Numerical predictions of turbulent flow have been obtained for RaH~0(10^9 -10^11 ) and compared with existing experimental data. A previously developed second moment closure model (Behnia et al. 1987) has been used to model the turbulence. Results indicate that a second moment closure model is capable of predicting the observed flow features.

Laminar and turbulent natural convection in rectangular cavities

"This chapter discusses laminar and turbulent natural convection in rectangular cavities. Natural convection in rectangular two-dimensional cavities has become a standard problem in numerical heat transfer because of its relevance in understanding a number of problems in engineering. Current research identified a number of difficulties with regard to the numerical methods and the turbulence modeling for this class of flows. Obtaining numerical predictions at high Rayleigh numbers proved computationally expensive such that results beyond Ra ∼ 1014 are rarely reported. The chapter discusses a study in which it was found that turbulent computations in square cavities can't be extended beyond Ra ∼ O (1012) despite having developed a code that proved very efficient for the high Ra laminar regime. As the Rayleigh number increased, thin boundary layers began to form next to the vertical walls, and the central region became progressively more stagnant and highly stratified. Results obtained for the high Ra laminar regime were in good agreement with existing studies. Turbulence computations, although of a preliminary nature, indicated that a second moment closure model was capable of predicting the experimentally observed flow features."--Publisher Summary

Tutoring with higher mathematics and the use of technology

We discuss three approaches to the use of technology as a teaching and learning tool that we are currently implementing for a target group of about one hundred second level engineering mathematics students. Central to these approaches is the underlying theme of motivating relatively poorly motivated students to learn, with the aim of improving learning outcomes. The approaches to be discussed have been used to replace, in part, more traditional mathematics tutorial sessions and lecture presentations. In brief, the first approach involves the application of constructivist thinking in the tertiary education arena, using technology as a computational and visual tool to create motivational knowledge conflicts or crises. The central idea is to model a realistic process of how scientific theory is actually developed, as proposed by Kuhn (1962), in contrast to more standard lecture and tutorial presentations. The second approach involves replacing procedural or algorithmic pencil-and-paper skills-consolidation exercises by software based tasks. Finally, the third approach aims at creating opportunities for higher order thinking via "on-line" exploratory or discovery mode tasks. The latter incorporates the incubation period method, as originally discussed by Rubinstein (1975) and others.

Flow field and traverse times for fan forced injection of fumigant via circular or annular inlet into stored grain

Fan forced injection of phosphine gas fumigant into stored grain is a common method to treat infestation by insects. For low injection velocities the transport of fumigant can be modelled as Darcy flow in a porous medium where the gas pressure satisfies Laplace's equation. Using this approach, a closed form series solution is derived for the pressure, velocity and streamlines in a cylindrically stored grain bed with either a circular or annular inlet, from which traverse times are numerically computed. A leading order closed form expression for the traverse time is also obtained and found to be reasonable for inlet configurations close to the central axis of the grain storage. Results are interpreted for the case of a representative 6m high farm wheat store, where the time to advect the phosphine to almost the entire grain bed is found to be approximately one hour.

Mathematical modelling and numerical simulation of phosphine flow during grain fumigation in leaky cylindrical silos

The phosphine distribution in a cylindrical silo containing grain is predicted. A three-dimensional mathematical model, which accounts for multicomponent gas phase transport and the sorption of phosphine into the grain kernel is developed. In addition, a simple model is presented to describe the death of insects within the grain as a function of their exposure to phosphine gas. The proposed model is solved using the commercially available computational fluid dynamics (CFD) software, FLUENT, together with our own C code to customize the solver in order to incorporate the models for sorption and insect extinction. Two types of fumigation delivery are studied, namely, fan- forced from the base of the silo and tablet from the top of the silo. An analysis of the predicted phosphine distribution shows that during fan forced fumigation, the position of the leaky area is very important to the development of the gas flow field and the phosphine distribution in the silo. If the leak is in the lower section of the silo, insects that exist near the top of the silo may not be eradicated. However, the position of a leak does not affect phosphine distribution during tablet fumigation. For such fumigation in a typical silo configuration, phosphine concentrations remain low near the base of the silo. Furthermore, we find that half-life pressure test readings are not an indicator of phosphine distribution during tablet fumigation.

Evaluation of k-epsilon and RNG turbulence models for confined swirling and non-swirling flow

In a very recent study [1] the Renormalisation Group (RNG) turbulence model was used to obtain flow predictions in a strongly swirling quarl burner, and was found to perform well in predicting certain features that are not well captured using less sophisticated models of turbulence. The implication is that the RNG approach should provide an economical and reliable tool for the prediction of swirling flows in combustor and furnace geometries commonly encountered in technological applications. To test this hypothesis the present work considers flow in a model furnace for which experimental data is available [2]. The essential features of the flow which differentiate it from the previous study [1] are that the annular air jet entry is relatively narrow and the base wall of the cylindrical furnace is at 90 degrees to the inlet pipe. For swirl numbers of order 1 the resulting flow is highly complex with significant inner and outer recirculation regions. The RNG and standard k-epsilon models are used to model the flow for both swirling and non-swirling entry jets and the results compared with experimental data [2]. Near wall viscous effects are accounted for in both models via the standard wall function formulation [3]. For the RNG model, additional computations with grid placement extending well inside the near wall viscous-affected sublayer are performed in order to assess the low Reynolds number capabilities of the model.

Computation of confined swirling flow using RNG and k-epsilon turbulence models

In this work we numerically model isothermal turbulent swirling flow in a cylindrical burner. Three versions of the RNG k-epsilon model are assessed against performance of the standard k-epsilon model. Sensitivity of numerical predictions to grid refinement, differing convective differencing schemes and choice of (unknown) inlet dissipation rate, were closely scrutinised to ensure accuracy. Particular attention is paid to modelling the inlet conditions to within the range of uncertainty of the experimental data, as model predictions proved to be significantly sensitive to relatively small changes in upstream flow conditions. We also examine the characteristics of the swirl--induced recirculation zone predicted by the models over an extended range of inlet conditions. Our main findings are: - (i) the standard k-epsilon model performed best compared with experiment; - (ii) no one inlet specification can simultaneously optimize the performance of the models considered; - (iii) the RNG models predict both single-cell and double-cell IRZ characteristics, the latter both with and without additional internal stagnation points. The first finding indicates that the examined RNG modifications to the standard k-e model do not result in an improved eddy viscosity based model for the prediction of swirl flows. The second finding suggests that tuning established models for optimal performance in swirl flows a priori is not straightforward. The third finding indicates that the RNG based models exhibit a greater variety of structural behaviour, despite being of the same level of complexity as the standard k-e model. The plausibility of the predicted IRZ features are discussed in terms of known vortex breakdown phenomena.

A heterogeneous computing approach to simulation of the Heston Stochastic Volatility Model

Stochastic volatility models are of fundamental importance to the pricing of derivatives. One of the most commonly used models of stochastic volatility is the Heston Model in which the price and volatility of an asset evolve as a pair of coupled stochastic differential equations. The computation of asset prices and volatilities involves the simulation of many sample trajectories with conditioning. The problem is treated using the method of particle filtering. While the simulation of a shower of particles is computationally expensive, each particle behaves independently making such simulations ideal for massively parallel heterogeneous computing platforms. In this paper, we present our portable Opencl implementation of the Heston model and discuss its performance and efficiency characteristics on a range of architectures including Intel cpus, Nvidia gpus, and Intel Many-Integrated-Core (mic) accelerators.

Computation of turbulent flow in an IFRF quarl burner

A computational model for isothermal axisymmetric turbulent flow in a quarl burner is set up using the CFD package FLUENT, and numerical solutions obtained from the model are compared with available experimental data. A standard k-e model and and two versions of the RNG k-e model are used to model the turbulence. One of the aims of the computational study is to investigate whether the RNG based k-e turbulence models are capable of yielding improved flow predictions compared with the standard k-e turbulence model. A difficulty is that the flow considered here features a confined vortex breakdown which can be highly sensitive to flow behaviour both upstream and downstream of the breakdown zone. Nevertheless, the relatively simple confining geometry allows us to undertake a systematic study so that both grid-independent and domain-independent results can be reported. The systematic study includes a detailed investigation of the effects of upstream and downstream conditions on the predictions, in addition to grid refinement and other tests to ensure that numerical error is not significant. Another important aim is to determine to what extent the turbulence model predictions can provide us with new insights into the physics of confined vortex breakdown flows. To this end, the computations are discussed in detail with reference to known vortex breakdown phenomena and existing theories. A major conclusion is that one of the RNG k-e models investigated here is able to correctly capture the complex forward flow region inside the recirculating breakdown zone. This apparently pathological result is in stark contrast to the findings of previous studies, most of which have concluded that either algebraic or differential Reynolds stress modelling is needed to correctly predict the observed flow features. Arguments are given as to why an isotropic eddy-viscosity turbulence model may well be able to capture the complex flow structure within the recirculating zone for this flow setup. With regard to the flow physics, a major finding is that the results obtained here are more consistent with the view that confined vortex breakdown is a type of axisymmetric boundary layer separation, rather than a manifestation of a subcritical flow state.

Analysis of turbulent confined highly swirling flow

We report here a CFD model of highly swirling flow in a quarl burner using three versions of the k-epsilon model. Results for the recirculating zone, the bounding shear layer and the downstream flow are presented. We discuss, with suitable qualifications, how the model predictions can inform our understanding of this class of flows.

Knowledge Exchange RDM Training & Skills Workshop, 9-10 Feb 2016, London

In today’s changing research environment, RDM is important in all stages of research. The skills and know-how in RDM that researchers and research support staff need, should be nurtured all though their career. At the end of 2015, KE initiated a project to compare approaches in RDM training within the partnership’s five member countries. The project was structured around two strands of activity: In the last months of 2015 a survey was conducted to collect information on current practice around RDM training, in order to provide an overview of the RDM training landscape. In February 2016 a workshop was held to share successful approaches to RDM training and capacity building provided within institutions and by infrastructure. The report describes the outputs of both the analysis of the survey and the outcomes of the workshop. The document provides an evidence base and informed suggestions to help improve RDM training practices in KE partner countries and beyond.

Programas de pagamento por serviços ambientais na gestão dos recursos hídricos: a experiência do FUNBOAS na bacia Lagos São João - RJ.

Os programas de Pagamento por Serviços Ambientais (PSA) correspondem a incentivos a quem ajuda a manter ou a produzir os serviços ecossistêmicos. O presente trabalho preocupou-se com a compreensão sobre a real contribuição destes programas de PSA para a gestão sustentável dos recursos hídricos, destacando o caso do programa FUNBOAS - Fundo de Boas Práticas Sócio Ambientais, criado pelo Comitê de Bacias Hidrográficas Lagos São João (CBHLSJ). Para atingir o objetivo proposto e entender os impactos do programa na microbacia, foram utilizadas informações fornecidas pelo CBHLSJ e foram realizadas visitas em campo, onde amostras de água foram coletadas para análise de parâmetros físico-químicos (pH, turbidez, NH3, P, OD, DQO e COT) e entrevistas com a população local foram realizadas. O FUNBOAS é um mecanismo de incentivo técnico e financeiro aos que ajudam a conservar os recursos hídricos e é alimentado com recursos oriundos da cobrança pelo uso da água da bacia. De uma forma geral, foi possível observar que o FUNBOAS possibilitou que os agricultores da microbacia beneficiada adotassem práticas de manejo mais sustentáveis, colaborando para a prestação dos serviços ecossistêmicos. Os resultados das análises de água indicaram que ainda existem pontos críticos de contaminação na microbacia, onde é necessário que medidas sejam tomadas para recuperar a qualidade da água do córrego. As entrevistas possibilitaram demonstrar que o FUNBOAS afetou o modo de vida dos produtores rurais da microbacia de forma positiva, melhorando a qualidade de vida das famílias beneficiadas. Também foi constatada a importância da preservação das florestas nativas na conservação da qualidade da água e dos trabalhos de educação ambiental na capacitação da população e na conscientização sobre a preservação dos recursos naturais. Além disso, o monitoramento ambiental da qualidade da água demonstrou ser uma ferramenta fundamental em programas de PSA em recursos hídricos, pois podem embasar a tomada de decisão, apontar pontos críticos de intervenção e evidenciar os impactos das ações, servindo para o acompanhamento dos resultados alcançados pelo programa. De uma forma geral, os resultados alcançados por esta pesquisa permitem concluir que programas de PSA são importantes instrumentos de gestão sustentável dos recursos hídricos, uma vez que possibilitam a preservação dos ecossistemas, colaborando com desenvolvimento socioeconômico das comunidades afetadas.

A global experiment on tilapia aquaculture: impacts of the Pond Dynamics/Aquaculture CRSP in Rwanda, Honduras, Philippines and Thailand

The Pond Dynamics/Aquaculture Collaborative Research Support Program (PDA/CRSP) is a global research network to generate basic science that may be used to advance aquaculture development. One of a family of research programs funded by the United States Agency for International Development (USAID), the CRSP focuses on improving the efficiency of aquaculture systems. The PDA/CRSP began work in 1982 in Thailand, and subsequently in the Philippines, Honduras, the US and, until recently, Rwanda. At all the sites, the goal is the same: to identify constraints to aquaculture production, and to design responses that are environmentally and culturally appropriate. The research network's global experiment has focused on tilapia (Oreochromis niloticus), although some sites have devoted attention to marine shrimp and other locally significant species. Impact of the network's investigations with tilapia is examined in this article.