Interaction between multiple bubbles in microchannel flow boiling and the effects on heat transfer

Flow boiling in microchannels have widely been studied in order to design more efficient cooling systems with numerical simulations forming a crucial part to deal with the areas that cannot be investigated experimentally. Previously published research largely focussed on the behaviour of a single bubble. Here, we focus on the behaviour of multiple bubbles. In this study, the influence of the distance between the bubbles (liquid slug length) is investigated in both an axisymmetric and a planar domain for two and three bubbles present. In this regard, an interface-tracking adaptive mesh refinement model was implemented to improve simulation time. Results show that the heat transfer was improved with sequential bubbles, and a 50% increase in heat transfer coefficient was observed for the cases investigated with three bubbles present. The heat transfer also improved the closer the bubbles were together.ThermaSMART Project of the European Commission, the Edinburgh Compute & Data Facility (ECDF) and Centre for High Performance Computing (CHPC), South Africa.https://www.elsevier.com/locate/ichmt2023-10-22hj2022Mechanical and Aeronautical Engineerin

Spanning Trees in Random Satisfiability Problems

Working with tree graphs is always easier than with loopy ones and spanning trees are the closest tree-like structures to a given graph. We find a correspondence between the solutions of random K-satisfiability problem and those of spanning trees in the associated factor graph. We introduce a modified survey propagation algorithm which returns null edges of the factor graph and helps us to find satisfiable spanning trees. This allows us to study organization of satisfiable spanning trees in the space spanned by spanning trees.Comment: 12 pages, 5 figures, published versio

Numerical and Experimental Investigations on Aerodynamic Behavior of the Ahmed Body Model with Different Diffuser Angles

Due to many restrictions applied by the necessity of fulfilling dimensional analysis in a numerical-experimental research and also the limits in experimental facilities a Low Reynolds Number simulation seems to be widespread. In this paper, effects of the diffuser angle on the aerodynamic behavior of the Ahmed body have been investigated for low Reynolds number flows. Numerical simulations were performed by solving the Reynolds Averaged Navier-Stokes (RANS) equations combined with different turbulence models. The Finite Volume Method (FVM) is used for simulations in Fluent 6.3.26 Software. The main objectives of the study are to improve the aerodynamic design of the body, analyzing the flow field to understand the nature of these improvements and reaching a suitable and reliable experimental-numerical setup for such a flow. Finally, it was concluded that the SST k-ω turbulence model with transitional flow corrections is the best choice. From the flow simulation and obtained experimental data, it was concluded that that drag coefficient is a function of three main phenomena. Results showed that the drag coefficient has its minimum value at a specific diffuser angle (8◦) and further increases in the angle lead to higher drag coefficient. On the other hand, the lift coefficient constantly decreases by increasing the diffuser angle. In order to show the validity of the numerical results, experimental data were obtained by measuring the drag and lift coefficients of scaled standard Ahmed body and a model with the diffuser angle of 8 degrees in a wind tunnel. Results confirmed that improvement of drag and lift coefficients occurs when diffuser region is considered for the Ahmed body. In addition, the flow field around the body was studied in detail to show the effects of the diffuser geometry on the aerodynamic characteristics of the body

Adaptive mesh refinement method for the reduction of computational costs while simulating slug flow

Microchannel flow boiling has been the focus of many experimental and numerical investigations due to the high heat transfer coefficients that it can induce. However, experimental research has been limited due to the small scales involved, leading researchers to employ computational fluid dynamics (CFD) simulations to resolve the dearth of research on microchannel flow boiling. Conventional CFD methods use a fine uniform mesh to capture the small scales and gradients, such as the liquid-vapour interface. This method has a large computational cost, and as a result, most research reported in the literature has been limited to two-dimensional axisymmetric domains. An interface-tracking adaptive mesh refinement model was created in this study to overcome the limitation of high computational costs without losing accuracy. This model dynamically refined the mesh only in the regions of interest and allowed a coarser mesh in the rest of the domain. This novel approach was able to recreate previously published results with a maximum error of 6.7%, while using less than 1.6% of the mesh elements. Several simulations were conducted in ANSYS Fluent 19.1 to determine the optimal settings for this new method to maintain accuracy and reduce cell count. These settings were determined as three levels of refinement (δL = 3), four refined cells on either side of the interface (δM = 4), and was implemented every five time steps (δT = 5). Finally, a case study was conducted to illustrate the possibility of simulating two-phase flow in microchannels in three dimensions with this method.ThermaSMART project of the European Commission, the Edinburgh Compute & Data Facility (ECDF) and the Centre for High Performance Computing (CHPC), South Africa.https://www.elsevier.com/locate/ichmt2023-10-23hj2022Mechanical and Aeronautical Engineerin

Chaos in Sandpile Models

We have investigated the "weak chaos" exponent to see if it can be considered as a classification parameter of different sandpile models. Simulation results show that "weak chaos" exponent may be one of the characteristic exponents of the attractor of \textit{deterministic} models. We have shown that the (abelian) BTW sandpile model and the (non abelian) Zhang model posses different "weak chaos" exponents, so they may belong to different universality classes. We have also shown that \textit{stochasticity} destroys "weak chaos" exponents' effectiveness so it slows down the divergence of nearby configurations. Finally we show that getting off the critical point destroys this behavior of deterministic models.Comment: 5 pages, 6 figure

Spatial Asymmetric Two dimensional Continuous Abelian Sandpile Model

We insert some asymmetries in the continuous Abelian sandpile models, such as directedness and ellipticity. We analyze probability distribution of different heights and also find the field theory corresponding to the models. Also we find the fields associated with some height variables.Comment: 14 Pages, 11 Figure

A novel computational approach to combine the optical and thermal modelling of Linear Fresnel Collectors using the finite volume method

A computational approach is presented, which uses the finite volume (FV) method in the Computational Fluid Dynamics (CFD) solver ANSYS Fluent to conduct the ray tracing required to quantify the optical performance of a line concentration Concentrated Solar Power (CSP) receiver, as well as the conjugate heat transfer modelling required to estimate the thermal efficiency of such a receiver. A Linear Fresnel Collector (LFC) implementation is used to illustrate the approach. It is shown that the Discrete Ordinates method can provide an accurate solution to the Radiative Transfer Equation (RTE) if the shortcomings of its solution are resolved appropriately in the FV CFD solver. The shortcomings are due to false scattering and the so-called ray effect inherent in the FV solution. The approach is first evaluated for a 2-D test case involving oblique collimated radiation and then for a more complex 2-D LFC optical domain based on the FRESDEMO project. For the latter, results are compared with and validated against those obtained with the Monte Carlo ray tracer, SolTrace. The outcome of the FV ray tracing in the LFC optical domain is mapped as a non-uniform heat flux distribution in the 3-D cavity receiver domain and this distribution is included in the FV conjugate heat transfer CFD model as a volumetric source. The result of this latter model is the determination of the heat transferred to the heat transfer fluid running in the collector tubes, thereby providing an estimation of the overall thermal efficiency. To evaluate the effectiveness of the phased approach in terms of accuracy and computational cost, the novel 2-D:3-D phased approach is compared with results of a fully integrated, but expensive 3-D optical and thermal model. It is shown that the less expensive model provides similar results and hence a large cost saving. The novel approach also provides the benefit of working in one simulation environment, i.e. ANSYS Workbench, where optimisation studies can be carried out to maximise the performance of linear CSP reflector layout and receiver configurations.University of Pretoria (South Africa) and the South African National Research Foundation (DST-NRF Solar Spoke).http://www.elsevier.com/locate/solener2016-06-30hb201

Annulus eccentricity optimisation of a phase-change material (PCM) horizontal double-pipe thermal energy store

The application of phase-change materials (PCMs) has received significant interest for use in thermal energy storage (TES) systems that can adjust the mismatch between the energy availability and demand. In the building sector, for example, PCMs can be used to reduce air-conditioning energy consumption by increasing the thermal capacity of the walls. However, as promising this technology may be, the poor thermal conductivity of PCMs has acted as a barrier to its commercialization, with many heat-transfer enhancement solutions proposed in the literature, such as microencapsulation or metal foam inserts, being either too costly and/or complex. The present study focuses on a low-cost and highly practical solution, in which natural-convective heat transfer is enhanced by placing the PCM in an eccentric annulus within a horizontal double-pipe TES heat exchanger. This paper presents an annulus-eccentricity optimisation study, whereby the optimal radial and tangential eccentricities are determined to minimize the charging time of a PCM thermal energy store. The storage performance of several geometrical configurations is predicted using a computational fluid dynamics (CFD) model based on the enthalpy-porosity formulation. The optimal geometrical configuration is then determined with response surface methods. The horizontal double-pipe heat exchanger studied considered here is an annulus filled with N-eicosane as the PCM for initial studies. In presence of N-eicosane, for the concentric configuration (which is the baseline case), the charging is completed at Fo = 0.64, while the charging of optimum eccentric geometries with the quickest and slowest charging is completed at Fo = 0.09 and Fo = 2.31, respectively. In addition, an investigation on the discharging performance of the studied configurations with N-eicosane shows the quickest discharge occurs with the concentric annulus case at Fo = 0.99, while the discharge time of the proposed optimum annuli is about three times this value. In other words, the proposed optimum geometry with the quickest charging time charges ~7.1 times faster but also discharges ~3 times slower, which is ideal for a TES, especially when used as passive thermal storage systems in nearly zero-emission buildings. Complementary studies demonstrate that the proposed optimum configuration improves the TES performance also when employing other PCM types as well as various shell-to-tube diameter ratios

Optimization of a trapezoidal cavity absorber for the Linear Fresnel Reflector

To increase the efficiency of Concentrated Solar Power (CSP) plants, the use of optimization methods is a current topic of research. This paper focuses on applying an integrated optimization technology to a solar thermal application, more specifically for the optimization of a trapezoidal cavity absorber of an LFR (Linear Fresnel Reflector), also called a Linear Fresnel Collector (LFC), CSP plant. LFR technology has been developed since the 1960s, and while large improvements in efficiencies have been made, there is still room for improvement. Once such area is in the receiver design where the optimal cavity shape, coatings, insulation thickness, absorber pipe selection, layout and spacing always need to be determined for a specific application. This paper uses a commercial tool to find an optimal design for a set of operating conditions. The objective functions that are used to judge the performance of a 2-D cavity are the combined heat loss through convection, conduction and radiation, as well as a wind resistance area. In this paper the effect of absorbed irradiation is introduced in the form of an outer surface of pipe temperature. Seven geometrical parameters are used as design variables. Based on a sample set requiring 79 CFD simulations, a global utopia point is found that minimizes both objectives. The most sensitive parameters were found to be the top insulation thickness and the cavity depth. Based on the results, the Multi-Objective Genetic Algorithm (MOGA) as contained in ANSYS DesignXplorer is shown to be effective in finding candidate optimal designs as well as the utopia point.University of Pretoria, South Africa, the South African National Research Foundation, as well as the Solar Spoke of the South African Department of Trade and Industry.http://www.elsevier.com/locate/solener2016-09-30hb201