Search results

In this work, a new two-phase version of the finite element-based Artificial Compressibility (AC) Characteristic-Based Split (CBS) algorithm is developed and applied for the first time to heat and mass transfer phenomena in porous media with associated phase change. The purpose of this study is to provide an alternative for the theoretical analysis and numerical simulation of multiphase transport phenomena in porous media. Traditionally, the more complex Separate Flow Model was used in which the vapour and liquid phases were considered as distinct fluids and mathematically described by the conservation laws for each phase separately, resulting in a large number of governing equations.

Even though the adopted mathematical model presents analogies with the conventional multicomponent mixture flow model, it is characterized by a considerable reduction in the number of the differential equations for the primary variables. The fixed-grid numerical formulation can be applied to the resolution of general problems that may simultaneously include a superheated vapour region, a two-phase zone and a sub-cooled liquid region in a single physical domain with irregular and moving phase interfaces in between. The local thermal non-equilibrium model is introduced to consider the heat exchange between fluid and solid within the porous matrix.

The numerical model is verified considering the transport phenomena in a homogenous and isotropic porous medium in which water is injected from one side and heated from the other side, where it leaves the computational domain in a superheated vapour state. Dominant forces are represented by capillary interactions and two-phase heat conduction. The obtained results have been compared with the numerical data available in the scientific literature.

The present algorithm provides a powerful routine tool for the numerical modelling of complex two-phase transport processes in porous media.

For the first time, the stabilized AC-CBS scheme is applied to the resolution of compressible viscous flow transport in porous materials with associated phase change. A properly stabilized matrix inversion-free procedure employs an adaptive local time step that allows acceleration of the solution process even in the presence of large source terms and low diffusion coefficients values (near the phase change point).

The purpose of this paper is to describe two‐ and three‐dimensional numerical modelling of solid oxide fuel cells (SOFCs) by employing an accurate and stable fully matrix inversion free finite element algorithm.

A general and detailed mathematical model has been developed for the description of the coupled complex phenomena occurring in fuel cells. A fully matrix inversion free algorithm, based on the artificial compressibility (AC) version of the characteristic‐based split (CBS) scheme and single domain approach have been successfully employed for the accurate and efficient simulation of high temperature SOFCs.

For the first time, a stable fully explicit algorithm has been applied to detailed multi‐dimensional simulation transport phenomena, coupled to chemical and electrochemical reactions, in fluid, porous and solid parts of a SOFC. The accuracy of the present results has been verified via comparison with experimental and numerical data available in the literature.

For the first time, thanks to a stabilization analysis conducted, the AC‐CBS algorithm has been successfully used to numerically solve the generalized model, applied in this paper to describe transport phenomena through free fluid channels and porous electrodes of SOFCs, without the need of further conditions at the fluid‐electrode interface.

The purpose of this paper is to introduce a robust mathematical model and finite element‐based numerical approach to solve solid oxide fuel cell (SOFC) problems.

A robust mathematical model is constructed by studying pros and cons of different SOFC and other fuel cell models. The finite element‐based numerical approach presented is a unified approach to solve multi‐disciplinary aspects arising from SOFC problems. The characteristic‐based split approach employed here is an efficient way of solving various flow, heat and mass transfer regimes in SOFCs.

The results presented show that both the model and numerical algorithm proposed are robust. Furthermore, the approaches proposed are general and can be easily extended to other similar problems of practical interest.

The model proposed is the first of this kind and the unified approach for solving flow, heat and mass transfer within a fuel cell is also novel.

The purpose of this paper is to focus on the numerical analysis of transient free convection heat transfer in partially porous cylindrical domains. The authors analyze the dependence of velocity and temperature fields on the geometry, by analyzing transient flow behavior for different values of cavity aspect ratio and radii ratio; both inner and outer radius are assumed variable in order to not change the difference r_o-r_i. Moreover, several Darcy numbers have been considered.

A dual time-stepping procedure based on the transient artificial compressibility version of the characteristic-based split algorithm has been adopted in order to solve the transient equations of the generalized model for heat and fluid flow through porous media. The present model has been validated against experimental data available in the scientific literature for two different problems, steady-state free convection in a porous annulus and transient natural convection in a porous cylinder, showing an excellent agreement.

For vertically divided half porous cavities, with Rayleigh numbers equal to 3.4×106 for the 4:1 cavity and 3.4×105 for the 8:1 cavity, the numerical results show that transient oscillations tend to disappear in presence of cylindrical geometry, differently from what happens for rectangular one. The magnitude of this phenomenon increases with radii ratio; the porous layer also affects the stability of velocity and temperature fields, as oscillations tend to decrease in presence of a porous matrix with lower value of the Darcy number.

A proper analysis of partially porous annular cavities is fundamental for the correct estimation of Nusselt numbers, as the formulas provided for rectangular domains are not able to describe these problems.

The proposed model represents a useful tool for the study of transient natural convection problems in porous and partially porous cylindrical and annular cavities, typical of many engineering applications. Moreover, a fully explicit scheme reduces the computational costs and ensures flexibility.

This is the first time that a fully explicit finite element scheme is employed for the solution of transient natural convection in partially porous tall annular cavities.

The purpose of the paper is to numerically simulate steady‐state thermo‐solutal convection in rectangular cavities with different aspect ratios, subject to horizontal temperature and concentration gradients, and validate the results against numerical and experimental data available from literature.

The fully explicit Artificial Compressibility (AC) version of the Characteristic Based Split (CBS) scheme is adopted to solve double diffusion (DD) problems. A stabilization analysis is carried out to efficiently solve the problems considered in the present work. The thermal and solutal buoyancy forces acting on the fluid have been taken into account in case of aiding and opposing flow conditions.

The stability limits derived by the authors for the thermo‐solutal convection assume a fundamental role to efficiently solve the DD problems considered. In the cases characterized by higher Rayleigh number the convergent solution is obtained only by employing the new stability conditions. The efficient matrix free procedure employed is a powerful tool to study complex DD problems.

In this paper, the authors extend the stabilization analysis for the AC‐CBS scheme to the solution of DD, fundamental to efficiently solve the present problems, and apply the present fully explicit matrix free scheme, based on finite elements, to the solution of DD natural convection in cavities.

The optimal design of a small dc torque motor, acting as a position actuator, is considered. The magnetic field is simulated by means of a 2‐D finite element model, taking into account the magnetization curve of iron and permanent magnets. The optimum shape of the motor, exhibiting the minimum cost, subject to a set of behavioural constraints, is found by a stochastic optimizer of lowest order.

The purpose of this paper is to compare the fluid dynamic performance of two Aqueous Humor (AH) ocular drainage devices, the SOLX® Gold Micro Shunt (GMS) and the novel Silicon Shunt Device (SSD), implanted by surgeons in human eyes to reduce the IntraOcular Pressure towards physiological values, by draining the AH from the Anterior Chamber to the Suprachoroidal Space, to cure eyes with glaucoma.

The generalized porous medium model is solved to simulate the AH flow through the two ocular drainage devices and the surrounding porous tissues of the eye.

In the GMS, probable stagnation regions have been found, due to the very small AH velocity values inside the device and to the surrounding tissues, creating possible blockage and malfunction of the device. The simple microtubular geometry of the novel SSD allows to have a regular AH flow and to choose shunts with different diameters and/or with the presence of radial holes, based on patient needs, with consequent reduction of post-operative complications.

The present model will be further developed taking into account the insertion of the present drainage devices inside the anterior section of the eye. The present results show the comparative fluid dynamic performance of the two shunts considered, and can be useful for surgeons to choose the adequate shunt, based on the required AH flow rate for a specific patient.

The present numerical approach, employing the generalized porous medium model, represents a useful tool to study the fluid dynamics of ocular drainage devices and to design these shunts, to reduce post-operative complications.

The generalized porous medium model is here applied for the first time to simulate the interaction of ocular drainage devices with the surrounding porous tissues of the eye.

This paper aims to provide a comprehensive literature review on modelling electro-osmotic flow in porous media.

Modelling electro-osmosis in fluid systems without solid particles has been first introduced. Then, after a brief description of the existing approaches for porous media modelling, electro-osmotic flow in porous media has been considered by analysing the main contributions to the development of this topic.

The analysis of literature has highlighted the absence of a universal model to analyse electro-osmosis in porous media, whereas many different methods and assumptions are used.

For the first time, the existing approaches for modelling electro-osmotic flow in porous have been collected and analysed to provide detailed indications for future works concerning this topic.

Modelling accurately the transient behaviour of natural convection flow in enclosures been a challenging task because of a variety of numerical errors which have limited achieving the higher order temporal accuracy. A fourth-order accurate finite difference method in both space and time is proposed to overcome these numerical errors and accurately model the transient behaviour of natural convection flow in enclosures using vorticity–streamfunction formulation.

Fourth-order wide stencil formula with appropriate one-sided difference extrapolation technique near the boundary is used for spatial discretisation, and classical fourth-order Runge–Kutta scheme is applied for transient term discretisation. The proposed method is applied on two transient case studies, i.e. convection–diffusion of a Gaussian Pulse and Taylor Vortex flow having analytical solution.

Error magnitude comparison and rate of convergence analysis of the proposed method with these analytical solutions establish fourth-order accuracy and prove the ability of the proposed method to truly capture the transient behaviour of incompressible flow. Also, to test the transient natural convection flow behaviour, the algorithm is tested on differentially heated square cavity at high Rayleigh number in the range of 10³-10⁸, followed by studying the transient periodic behaviour in a differentially heated vertical cavity of aspect ratio 8:1. An excellent comparison is obtained with standard benchmark results.

The developed method is applied on 2D enclosures; however, the present methodology can be extended to 3D enclosures using velocity–vorticity formulations which shall be explored in future.

The proposed methodology to achieve fourth-order accurate transient simulation of natural convection flows is novel, to the best of the authors’ knowledge. Stable fourth-order vorticity boundary conditions are derived for boundary and external boundary regions. The selected case studies for comparison demonstrate not only the fourth-order accuracy but also the considerable reduction in error magnitude by increasing the temporal accuracy. Also, this study provides novel benchmark results at five different locations within the differentially heated vertical cavity of aspect ratio 8:1 for future comparison studies.

The purpose of this paper is to numerically investigate steady, laminar natural and mixed convection heat transfer in a two-dimensional cavity by using a finite volume method with a fourth-order approximation of convective terms, with and without the presence of nanoparticles. Highly accurate benchmark results are also provided.

A finite volume method on a non-uniform staggered grid is used for the solution of two-dimensional momentum and energy conservation equations. Diffusion terms, in the momentum and energy equations, are approximated using second-order central differences, whereas a non-uniform four-point fourth-order interpolation (FPFOI) scheme is developed for the convective terms. Coupled mass and momentum conservation equations are solved iteratively using a semi-implicit method for pressure-linked equation method.

For the case of natural convection problem at high-Rayleigh numbers, grid density must be sufficiently high in order to obtain grid-independent results and capture reality of the physics. Heat transfer enhancement for natural convection is observed up to a certain value of the nanoparticle volume fraction. After that value, heat transfer deterioration is found with increasing nanoparticle volume fraction.

Developed a non-uniform FPFOI scheme. Highly accurate benchmark results for the heat transfer of Al₂O₃-water nanofluid in a cavity are provided.