A level-set method for thermal motion of bubbles and droplets

Published under licence in Journal of Physics: Conference Series by IOP Publishing Ltd. Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.A conservative level-set model for direct simulation of two-phase flows with thermocapillary effects at dynamically deformable interface is presented. The Navier-Stokes equations coupled with the energy conservation equation are solved by means of a finite-volume/level-set method. Some numerical examples including thermocapillary motion of single and multiple fluid particles are computed by means of the present method. The results are compared with analytical solutions and numerical results from the literature as validations of the proposed model.Peer ReviewedPostprint (published version

Evaluation of the Lubrication Regime for Rotary Compressors. Influence of Thermal Expansions on the Minimum Film Thickness.

In this work, the lubrication regime of a rotary compressor is studied in the regions of the roller, cylinder, shaft, and bearings. A code is programmed so it can solve multiple cases with different suction-discharge pressures, angular velocity, viscosity, and geometrical conditions, among others. In order to study the lubrication regime, the balance of the different forces that act on the roller is computed. To calculate the oil pressure forces, the Reynolds equation coupled with an iterative method is used to find the minimal distance and pressure distribution that balances the external forces. This process is repeated for a discrete number of rotation angles to evaluate the overall rotation of a cycle. In addition, the effect that the thermal expansion of the different materials can have on the lubrication regime is studied. If the materials get too hot, the distances between surfaces might increase or reduce with respect to the distances at room temperature. These changes might lead to a point where the lubrication regime changes. The objective is to include this thermal expansion effect in the model and evaluate how it can affect the lubrication regime.The author Jordi Vera has been financially supported by the Ministerio de Educación y Ciencia (MEC), Spain, (FPI grant PRE2018-084017). The author E.Schillaci acknowledges the financial support of the Programa Torres Quevedo (PTQ2018-010060).Postprint (published version

A level-set model for thermocapillary motion of deformable fluid particles

A new level-set model is proposed for simulating immiscible thermocapillary flows with variable fluid-property ratios at dynamically deformable interfaces. The Navier–Stokes equations coupled with the energy conservation equation are solved by means of a finite-volume/level-set approach, adapted to a multiple marker methodology in order to avoid the numerical coalescence of the fluid particles. The temperature field is coupled to the surface tension through an equation of state. Some numerical examples including thermocapillary driven convection in two superimposed fluid layers, and thermocapillary motion of single and multiple fluid particles are computed using the present method. These results are compared against analytical solutions and numerical results from the literature as validations of the proposed model.Peer ReviewedPostprint (author's final draft

Analysis and design of a drain water heat recovery storage unit based on PCM plates

© 2016. This version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/This paper is focused on the detailed analysis of a PCM plate heat storage unit with a particular configuration, looking for the maximum contact area with the fluid (water) and the minimum volume to be used in a real household application. In that sense, a numerical study of the thermal and fluid dynamic behaviour of the water flow and the PCM melting-solidification processes, together with the thermal behaviour of the solid elements of the unit, has been carried out. On the other hand, an experimental set-up has been designed and built to validate the numerical model and characterise the performance of the heat storage unit. The purpose of the numerical and experimental study is to present a series of results to describe the heat storage unit performance in function of the time. Thus, after a preliminary design study three different cases have been simulated and tested. A 7.2% of discrepancy between numerical results and experimental data has been evaluated for the heat transfer. The PCM heat storage unit designed is capable to store approx. 75% of the thermal energy from the previous process wasted water heat, and recover part of it to supply around 50% of the thermal energy required to heat up the next process.Peer ReviewedPostprint (author's final draft

Numerical dynamic analysis of reciprocating compressor mechanism. Parametric studies for optimization purposes

© 2016. This version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/A complete numerical dynamic analysis of reciprocating compressor mechanism is presented, coupling the instantaneous pressure in the compression chamber, the electric motor torque and the hydrodynamic reactions, which arise from the piston and crankshaft secondary movements. Additionally, non-constant crankshaft angular velocity and the piston and crankshaft misalignment torques have also been considered. Two sensitivity analyses have been carried out to prove that neither the inertial forces in the directions of the secondary movements, nor the oscillations of the angular velocity produce significant differences in the compressor behaviour. Finally, a set of parametric studies has been developed to evaluate the influence of geometrical parameters in the stability of the secondary movements, the friction power losses and the compressor consumptionPeer ReviewedPostprint (author's final draft

On the Use of a Parallel Object-Oriented Code for Solving the Heat Transfer in Hermetic Reciprocating Compressors

The heat transport phenomenon in a hermetic reciprocating compressor is addressed in this work. This is far from straightforward. It involves several transient physical phenomena interacting to each other. The heat is exchanged between the refrigerant fluid and the solid parts of the compressor (suction muffler, cylinder head, crankcase, etc.). At the same time, the solid parts exchange heat to each other by means of conduction and radiation. Moreover, the phenomenon happens in non-symmetrical complex geometries and the solid parts are made of different materials. This is interesting from both the software engineering and the compressor design viewpoint. A parallel object-oriented software platform for the resolution of multiphysics problems is employed. This platform allows the use of partitioned strategies so that the compressor heat transport problem -a global problem- can be divided into several smaller parts -local problems-. This makes possible the use of multilevel modeling strategies for thermal systems analysis. Furthermore, in order to couple the several sub-problems in an integrated simulation, the platform provides data transfer tools -for matching and non-matching meshes- to exchange sub-domain state information. In particular, the work provides detailed information on the heat distribution and the temperature of the components of a test compressor. By means of comparative studies the thermal properties of some of its components are analyzed. This highlights the importance of choosing proper materials. For example, different suction muffler materials are tested to investigate their influence on the volumetric efficiency. Since the whole compressor is simulated, the consequences of altering specific component properties are also appreciated on the other components. In sum, the work presents illustrative numerical results of the three-dimensional heat transfer in a compressor that show the potential use of computer simulation to support design of components to attain feasibility and energy efficiency

Simulation of fluid-structure interaction and impact force on a reed valve

The cyclic impact force between a reed valve and the seat plate is the main reason of the valve failure in many thermo-technical devices as compressors, engines, etc. According to experimental observations the latter is due to fatigue and usually occurs in the leading part of the valve ‘neck’. In this work, a complex numerical analysis is presented aimed to studying the external forces and internal stresses suffered by the valve. In particular, the impact force between the valve and the seat is studied. The numerical analysis relies on the coupled synergy of two different simulation concepts. In order to do so, two codes are used: (1) first, the in-house Computational Fluid Dynamics (CFD) code presented in [1] is employed to simulate the Fluid-Structure Interaction (FSI) between gas and valve, extracting reference data for valve displacement and external gas pressures; (2) second, the analysis of the internal structure stresses, together with the impact forces with the plate is implemented in a Computational Solid Dynamics (CSD) code developed in FreeFEM++ [2]. The impact force representation is based on the formulation presented in [3] where a conserving algorithm for frictionless dynamic contact/impact is developed. Due to the importance of obtaining an adequate impact force, an exhaustive study is carried out on its characterization in terms of numerical parameters, such as the penalty stiffness. Under this framework, the valve displacement and impact velocities are verified. Hence, impact forces are analysed in different scenarios, obtaining interesting observations about stresses distribution, with a particular focus on the points where failure is experienced.The authors acknowledge Voestalpine Precision Strip AB company for the previous research collaboration project that allowed to validate experimentally the presented numerical methods. P. Castrillo gratefully acknowledges the Universitat Politecnica de Catalunya and Banco Santander for the financial ` support of his predoctoral grant FPI-UPC (109 FPI-UPC 2018). E. Schillaci acknowledges the financial support of the Programa Torres Quevedo (PTQ2018-010060). This work has also been financially supported by a competitive R+D project (ENE2017-88697-R) by the Spanish Research Agency.Postprint (published version

Analysis of high-order interpolation schemes for solving linear problems in unstructured meshes using the finite volume method

Finite-volume strategies in fluid-structure interaction problems would be of crucialvimportance in many engineering applications such as in the analysis of reed valves in reciprocating compressors. The efficient implementation of this strategy passes from the formulation of reliable high-order schemes on 3D unstructured meshes. The development of high-order models is essential in bending-dominant problems, where the phenomenon of shear blocking appears. In order to solve this problem, it is possible to either increase the number of elements or increase the interpolation order of the main variable. Increasing the number of elements does not always yield good results and implies a very high computational cost that, in real problems, is inadmissible. Using unstructured meshes is also vital because they are necessary for real problems where the geometries are complex and depart from canonical rectangular or regular shapes. This work presents a series of tests to demonstrate the feasibility of a high-order model using finite volumes for linear elasticity on unstructured and structured meshes. The high-order interpolation will be performed using two different schemes such as the Moving Least Squares (MLS) and the Local Regression Estimators (LRE). The reliability of the method for solving 2D and 3D problems will be verified by solving some known test cases with an analytical solution such as a thin beam or problems where stress concentrations appear.P. Castrillo gratefully acknowledges the Universitat Politècnica de Catalunya and Banco Santander for the financial support of his predoctoral grant FPI-UPC (109 FPI-UPC 2018). The authors are supported by the Ministerio de Economía y Competitividad, Spain, RETOtwin project (PDC2021-120970-I00).Peer ReviewedPostprint (published version