Articles de revista
http://hdl.handle.net/2117/3191
Fri, 18 Aug 2017 14:33:05 GMT
20170818T14:33:05Z

Numerical study of Taylor bubbles rising in a stagnant liquid using a levelset/movingmesh method
http://hdl.handle.net/2117/105970
Numerical study of Taylor bubbles rising in a stagnant liquid using a levelset/movingmesh method
Gutiérrez González, Ernesto; Balcázar Arciniega, Néstor; Bartrons Casademont, Eduard; Rigola Serrano, Joaquim
An Arbitrary LagrangianEulerian formulation has been posed to solve the challenging problem of the threedimensional Taylor bubble, within a Conservative Level Set (CLS) framework. By employing a domain optimization method (i.e. the moving mesh method), smaller domains can be used to simulate rising bubbles, thus saving computational resources. As the method requires the use of open boundaries, a careful treatment of both inflow and outflow boundary conditions has been carried out. The coupled CLS  moving mesh method has been verified by means of extensive numerical tests. The challenging problem of the full threedimensional Taylor bubble has then been thoroughly addressed, providing a detailed description of its features. The study also includes a sensitivity analyses with respect to the initial shape of the bubble, the initial volume of the bubble, the flow regime and the inclination of the channel.
Thu, 29 Jun 2017 08:29:59 GMT
http://hdl.handle.net/2117/105970
20170629T08:29:59Z
Gutiérrez González, Ernesto
Balcázar Arciniega, Néstor
Bartrons Casademont, Eduard
Rigola Serrano, Joaquim
An Arbitrary LagrangianEulerian formulation has been posed to solve the challenging problem of the threedimensional Taylor bubble, within a Conservative Level Set (CLS) framework. By employing a domain optimization method (i.e. the moving mesh method), smaller domains can be used to simulate rising bubbles, thus saving computational resources. As the method requires the use of open boundaries, a careful treatment of both inflow and outflow boundary conditions has been carried out. The coupled CLS  moving mesh method has been verified by means of extensive numerical tests. The challenging problem of the full threedimensional Taylor bubble has then been thoroughly addressed, providing a detailed description of its features. The study also includes a sensitivity analyses with respect to the initial shape of the bubble, the initial volume of the bubble, the flow regime and the inclination of the channel.

A 3D volumeoffluid advection method based on cellvertex velocities for unstructured meshes
http://hdl.handle.net/2117/105736
A 3D volumeoffluid advection method based on cellvertex velocities for unstructured meshes
Jofre Cruanyes, Lluís; Lehmkuhl Barba, Oriol; Castro González, Jesús; Oliva Llena, Asensio
A new geometrical VolumeofFluid (VOF) method for capturing interfaces on threedimensional (3D) Cartesian and unstructured meshes is introduced. The method reconstructs interfaces as first and secondorder piecewise planar approximations (PLIC), and advects volumes in a single unsplit Lagrangian–Eulerian (LE) geometrical algorithm based on constructing flux polyhedrons by tracing back the Lagrangian trajectories of the cellvertex velocities. In this way, the situations of overlapping between flux polyhedrons are minimized, consequently, the accuracy in the solution of the advection equation is improved by minimizing the creation of overshoots (volume fractions over one), undershoots (volume fractions below zero) and wisps (fluid in void regions or vice versa). However, if not treated carefully, the use of cellvertex velocities may result in the construction of flux polyhedrons that contain nonplanar faces and that do not conserve volume. Therefore, this work explains in detail a set of geometric algorithms necessary to overcome these two drawbacks. In addition, the new VOF method is analyzed numerically on 3D Cartesian and unstructured meshes, first, by reconstructing the interface of spherical geometries and, second, by evaluating the final advection result of a sphere placed in a rotation, shear and deformation field.
Thu, 22 Jun 2017 14:39:33 GMT
http://hdl.handle.net/2117/105736
20170622T14:39:33Z
Jofre Cruanyes, Lluís
Lehmkuhl Barba, Oriol
Castro González, Jesús
Oliva Llena, Asensio
A new geometrical VolumeofFluid (VOF) method for capturing interfaces on threedimensional (3D) Cartesian and unstructured meshes is introduced. The method reconstructs interfaces as first and secondorder piecewise planar approximations (PLIC), and advects volumes in a single unsplit Lagrangian–Eulerian (LE) geometrical algorithm based on constructing flux polyhedrons by tracing back the Lagrangian trajectories of the cellvertex velocities. In this way, the situations of overlapping between flux polyhedrons are minimized, consequently, the accuracy in the solution of the advection equation is improved by minimizing the creation of overshoots (volume fractions over one), undershoots (volume fractions below zero) and wisps (fluid in void regions or vice versa). However, if not treated carefully, the use of cellvertex velocities may result in the construction of flux polyhedrons that contain nonplanar faces and that do not conserve volume. Therefore, this work explains in detail a set of geometric algorithms necessary to overcome these two drawbacks. In addition, the new VOF method is analyzed numerically on 3D Cartesian and unstructured meshes, first, by reconstructing the interface of spherical geometries and, second, by evaluating the final advection result of a sphere placed in a rotation, shear and deformation field.

A lowdissipation convection scheme for the stable discretization of turbulent interfacial flow
http://hdl.handle.net/2117/105415
A lowdissipation convection scheme for the stable discretization of turbulent interfacial flow
Schillaci, Eugenio; Jofre Cruanyes, Lluís; Balcázar Arciniega, Néstor; Antepara Zambrano, Óscar; Oliva Llena, Asensio
This paper analyzes a lowdissipation discretization for the resolution of immiscible, incompressible multiphase flow by means of interfacecapturing schemes. The discretization is built on a threedimensional, unstructured finitevolume framework and aims at minimizing the differences in kinetic energy preservation with respect to the continuous governing equations. This property plays a fundamental role in the case of flows presenting significant levels of turbulence. At the same time, the hybrid form of the convective operator proposed in this work incorporates localized lowdispersion characteristics to limit the growth of spurious flow solutions. The lowdissipation discrete framework is presented in detail and, in order to expose the advantages with respect to commonly used methodologies, its conservation properties and accuracy are extensively studied, both theoretically and numerically. Numerical tests are performed by considering a threedimensional vortex, an exact sinusoidal function, and a spherical drop subjected to surface tension forces in equilibrium and immersed in a swirling velocity field. Finally, the turbulent atomization of a liquidgas jet is numerically analyzed to further assess the capabilities of the method.
© 2017. This version is made available under the CCBYNCND 4.0 license http://creativecommons.org/licenses/byncnd/4.0/
Wed, 14 Jun 2017 12:46:47 GMT
http://hdl.handle.net/2117/105415
20170614T12:46:47Z
Schillaci, Eugenio
Jofre Cruanyes, Lluís
Balcázar Arciniega, Néstor
Antepara Zambrano, Óscar
Oliva Llena, Asensio
This paper analyzes a lowdissipation discretization for the resolution of immiscible, incompressible multiphase flow by means of interfacecapturing schemes. The discretization is built on a threedimensional, unstructured finitevolume framework and aims at minimizing the differences in kinetic energy preservation with respect to the continuous governing equations. This property plays a fundamental role in the case of flows presenting significant levels of turbulence. At the same time, the hybrid form of the convective operator proposed in this work incorporates localized lowdispersion characteristics to limit the growth of spurious flow solutions. The lowdissipation discrete framework is presented in detail and, in order to expose the advantages with respect to commonly used methodologies, its conservation properties and accuracy are extensively studied, both theoretically and numerically. Numerical tests are performed by considering a threedimensional vortex, an exact sinusoidal function, and a spherical drop subjected to surface tension forces in equilibrium and immersed in a swirling velocity field. Finally, the turbulent atomization of a liquidgas jet is numerically analyzed to further assess the capabilities of the method.

Large eddy simulation of a turbulent diffusion flame: some aspects of subgrid modelling consistency
http://hdl.handle.net/2117/104310
Large eddy simulation of a turbulent diffusion flame: some aspects of subgrid modelling consistency
Ventosa Molina, Jordi; Lehmkuhl, Oriol; Pérez Segarra, Carlos David; Oliva Llena, Asensio
In the context of Large Eddy Simulation (LES) solely for the momentum
transport equation there may be found several models for the turbulent subgrid fluxes. Furthermore, among those relying on the eddy diffusivity approach, each model may be based on different invariants of the strain rate. Besides, when heat and mass transfer are also considered, closures for the subgrid turbulent scalar fluxes are also required. Hence, different model combinations may be considered. Additionally, when other physical phenomena are included, such as combustion, further subgrid modelling is involved. Therefore, in the present study a LES simulation of a turbulent diffusion flame is performed and different combination of subgrid models are used in order to analyse the numerical effects in the simulations. Several models for the turbulent momentum subgrid fluxes are considered, both constant and dynamically evaluated Schmidt numbers. Regarding combustion, in the context of the Flamelet/ProgressVariable (FPV) model, with an assumed probability density function for the turbulentchemistry interactions and four different closures for the subgrid mixture fraction variance are considered. Hence, a large number of model combinations are possible. The present study highlights the need for a consistent closure of subgrid effects. It is shown that, in the context of an FPV modelling, incorrect capture of subgrid mixing results in a flame liftoff for the studied flame (DLR A diffusion flame), even though experimentally an attached flame was reported. It is found that a consistent formulation is required, that is, all subgrid closures should become active in the same regions of the domain to avoid modelling inconsistencies. In contrast, when the classical flamelet approach is used, no liftoff is observed. The reason is that the classical flamelet includes only a limited subset of the possible flame states, i.e. only includes burning flamelets and extinguished flamelets for scalar dissipation rates past the extinction one.
This is a copy of the author 's final draft version of an article published in the journal Flow turbulence and combustion. The final publication is available at Springer via http://dx.doi.org/10.1007/s1049401798132
Thu, 11 May 2017 10:22:17 GMT
http://hdl.handle.net/2117/104310
20170511T10:22:17Z
Ventosa Molina, Jordi
Lehmkuhl, Oriol
Pérez Segarra, Carlos David
Oliva Llena, Asensio
In the context of Large Eddy Simulation (LES) solely for the momentum
transport equation there may be found several models for the turbulent subgrid fluxes. Furthermore, among those relying on the eddy diffusivity approach, each model may be based on different invariants of the strain rate. Besides, when heat and mass transfer are also considered, closures for the subgrid turbulent scalar fluxes are also required. Hence, different model combinations may be considered. Additionally, when other physical phenomena are included, such as combustion, further subgrid modelling is involved. Therefore, in the present study a LES simulation of a turbulent diffusion flame is performed and different combination of subgrid models are used in order to analyse the numerical effects in the simulations. Several models for the turbulent momentum subgrid fluxes are considered, both constant and dynamically evaluated Schmidt numbers. Regarding combustion, in the context of the Flamelet/ProgressVariable (FPV) model, with an assumed probability density function for the turbulentchemistry interactions and four different closures for the subgrid mixture fraction variance are considered. Hence, a large number of model combinations are possible. The present study highlights the need for a consistent closure of subgrid effects. It is shown that, in the context of an FPV modelling, incorrect capture of subgrid mixing results in a flame liftoff for the studied flame (DLR A diffusion flame), even though experimentally an attached flame was reported. It is found that a consistent formulation is required, that is, all subgrid closures should become active in the same regions of the domain to avoid modelling inconsistencies. In contrast, when the classical flamelet approach is used, no liftoff is observed. The reason is that the classical flamelet includes only a limited subset of the possible flame states, i.e. only includes burning flamelets and extinguished flamelets for scalar dissipation rates past the extinction one.

DNS and regularization modeling of a turbulent differentially heated cavity of aspect ratio 5
http://hdl.handle.net/2117/103914
DNS and regularization modeling of a turbulent differentially heated cavity of aspect ratio 5
Trias Miquel, Francesc Xavier; Gorobets, Andrei; Oliva Llena, Asensio; Pérez Segarra, Carlos David
This work is devoted to the study of turbulent natural convection flows in
differentially heated cavities. The adopted configuration corresponds to an airfilled
(Pr = 0.7) cavity of aspect ratio 5 and Rayleigh number Ra = 4.5 × 1010
(based on the cavity height). Firstly, a complete direct numerical simulation
(DNS) has been performed. Then, the DNS results have been used as reference
solution to assess the performance of symmetrypreserving regularization as a
simulation shortcut: a novel class of regularization that restrain the convective
production of small scales of motion in an unconditionally stable manner. In
this way, the new set of equations is dynamically less complex than the original
NavierStokes equations, and therefore more amenable to be numerically solved.
Direct comparison with the DNS results shows fairly good agreement even for very coarse grids.
Tue, 02 May 2017 14:37:13 GMT
http://hdl.handle.net/2117/103914
20170502T14:37:13Z
Trias Miquel, Francesc Xavier
Gorobets, Andrei
Oliva Llena, Asensio
Pérez Segarra, Carlos David
This work is devoted to the study of turbulent natural convection flows in
differentially heated cavities. The adopted configuration corresponds to an airfilled
(Pr = 0.7) cavity of aspect ratio 5 and Rayleigh number Ra = 4.5 × 1010
(based on the cavity height). Firstly, a complete direct numerical simulation
(DNS) has been performed. Then, the DNS results have been used as reference
solution to assess the performance of symmetrypreserving regularization as a
simulation shortcut: a novel class of regularization that restrain the convective
production of small scales of motion in an unconditionally stable manner. In
this way, the new set of equations is dynamically less complex than the original
NavierStokes equations, and therefore more amenable to be numerically solved.
Direct comparison with the DNS results shows fairly good agreement even for very coarse grids.

Three dimensionality in the wake of the flow around a circular cylinder at Reynolds number 5000
http://hdl.handle.net/2117/103891
Three dimensionality in the wake of the flow around a circular cylinder at Reynolds number 5000
Aljure Osorio, David E.; Lehmkuhl, Oriol; Rodríguez Pérez, Ivette María; Oliva Llena, Asensio
The turbulent flow around a circular cylinder has been investigated at Re=5000Re=5000 using direct numerical simulations. Low frequency behavior, vortex undulation, vortex splitting, vortex dislocations and three dimensional flow within the wake were found to happen at this flow regime. In order to successfully capture the wake three dimensionality, different spanwise lengths were considered. It was found that a length LZ=2pDLZ=2pD was enough to capture this behavior, correctly predicting different aspects of the flow such as drag coefficient, Strouhal number and pressure and velocity distributions when compared to experimental values. Two instability mechanisms were found to coexist in the present case study: a global type instability originating in the shear layer, which shows a characteristic frequency, and a convective type instability that seems to be constantly present in the near wake. Characteristics of both types of instabilities are identified and discussed in detail. As suggested by Norberg, a resonancetype effect takes place in the vortex formation region, as the coexistence of both instability mechanisms result in distorted vortex tubes. However, vortex coherence is never lost within the wake.
Tue, 02 May 2017 10:41:23 GMT
http://hdl.handle.net/2117/103891
20170502T10:41:23Z
Aljure Osorio, David E.
Lehmkuhl, Oriol
Rodríguez Pérez, Ivette María
Oliva Llena, Asensio
The turbulent flow around a circular cylinder has been investigated at Re=5000Re=5000 using direct numerical simulations. Low frequency behavior, vortex undulation, vortex splitting, vortex dislocations and three dimensional flow within the wake were found to happen at this flow regime. In order to successfully capture the wake three dimensionality, different spanwise lengths were considered. It was found that a length LZ=2pDLZ=2pD was enough to capture this behavior, correctly predicting different aspects of the flow such as drag coefficient, Strouhal number and pressure and velocity distributions when compared to experimental values. Two instability mechanisms were found to coexist in the present case study: a global type instability originating in the shear layer, which shows a characteristic frequency, and a convective type instability that seems to be constantly present in the near wake. Characteristics of both types of instabilities are identified and discussed in detail. As suggested by Norberg, a resonancetype effect takes place in the vortex formation region, as the coexistence of both instability mechanisms result in distorted vortex tubes. However, vortex coherence is never lost within the wake.

Heat and moisture insulation by means of air curtains: application to refrigerated chambers
http://hdl.handle.net/2117/103512
Heat and moisture insulation by means of air curtains: application to refrigerated chambers
Giráldez García, Héctor; Pérez Segarra, Carlos David; Oliet Casasayas, Carles; Oliva Llena, Asensio
The present study is devoted to the determination of the efficiency of air curtain units (ACUs) applied to heat and moisture insulation of refrigerated chambers. A detailed study of the fluid dynamics and heat and mass transfer of the ACU in the refrigerated space and the external ambient is carried out by means of large eddy simulations (LES). The heat and moisture entrainment through the doorway and their transport inside the inner space are fully described. Three different configurations are studied: nonrecirculating, recirculating and twinjet air curtains. The condensation produced in the cool walls of the refrigerated space is evaluated considering the warm humid air from the ambient which penetrates inside the chamber through the doorway. The influence of both the discharge velocities and the discharge angles on the sealing capabilities of the three different tested ACU configurations is determined.
Tue, 18 Apr 2017 11:30:21 GMT
http://hdl.handle.net/2117/103512
20170418T11:30:21Z
Giráldez García, Héctor
Pérez Segarra, Carlos David
Oliet Casasayas, Carles
Oliva Llena, Asensio
The present study is devoted to the determination of the efficiency of air curtain units (ACUs) applied to heat and moisture insulation of refrigerated chambers. A detailed study of the fluid dynamics and heat and mass transfer of the ACU in the refrigerated space and the external ambient is carried out by means of large eddy simulations (LES). The heat and moisture entrainment through the doorway and their transport inside the inner space are fully described. Three different configurations are studied: nonrecirculating, recirculating and twinjet air curtains. The condensation produced in the cool walls of the refrigerated space is evaluated considering the warm humid air from the ambient which penetrates inside the chamber through the doorway. The influence of both the discharge velocities and the discharge angles on the sealing capabilities of the three different tested ACU configurations is determined.

Simulation of the twofluid model on incompressible flow with Fractional Step method for both resolved and unresolved scale interfaces
http://hdl.handle.net/2117/103410
Simulation of the twofluid model on incompressible flow with Fractional Step method for both resolved and unresolved scale interfaces
Hou, Xiaofei; Rigola Serrano, Joaquim; Lehmkuhl Barba, Oriol; Oliva Llena, Asensio
In the present paper, the Fractional Step method usually used in single fluid flow is here extended and applied for the twofluid model resolution using the finite volume discretization. The use of a projection method resolution instead of the usual pressurecorrection method for multifluid flow, successfully avoids iteration processes. On the other hand, the main weakness of the two fluid model used for simulations of free surface flows, which is the numerical diffusion of the interface, is also solved by means of the conservative Level Set method (interface sharpening) (Strubelj et al., 2009). Moreover, the use of the algorithm proposed has allowed presenting different freesurface cases with or without Level Set implementation even under coarse meshes under a wide range of density ratios. Thus, the numerical results presented, numerically verified, experimentally validated and converged under high density ratios, shows the capability and reliability of this resolution method for both mixed and unmixed flows.
Thu, 06 Apr 2017 09:19:18 GMT
http://hdl.handle.net/2117/103410
20170406T09:19:18Z
Hou, Xiaofei
Rigola Serrano, Joaquim
Lehmkuhl Barba, Oriol
Oliva Llena, Asensio
In the present paper, the Fractional Step method usually used in single fluid flow is here extended and applied for the twofluid model resolution using the finite volume discretization. The use of a projection method resolution instead of the usual pressurecorrection method for multifluid flow, successfully avoids iteration processes. On the other hand, the main weakness of the two fluid model used for simulations of free surface flows, which is the numerical diffusion of the interface, is also solved by means of the conservative Level Set method (interface sharpening) (Strubelj et al., 2009). Moreover, the use of the algorithm proposed has allowed presenting different freesurface cases with or without Level Set implementation even under coarse meshes under a wide range of density ratios. Thus, the numerical results presented, numerically verified, experimentally validated and converged under high density ratios, shows the capability and reliability of this resolution method for both mixed and unmixed flows.

On the evolution of flow topology in turbulent RayleighBénard convection
http://hdl.handle.net/2117/103099
On the evolution of flow topology in turbulent RayleighBénard convection
Dabbagh, Firas; Trias Miquel, Francesc Xavier; Gorobets, Andrei; Oliva Llena, Asensio
Smallscale dynamics is the spirit of turbulence physics. It implicates many attributes of flow topology evolution, coherent structures, hairpin vorticity dynamics, and mechanism of the kinetic energy cascade. In this work, several dynamical aspects of the smallscale motions have been numerically studied in a framework of RayleighBenard convection (RBC). To do so, direct numerical simulations have been carried out at two Rayleigh numbers Ra = 10(8) and 10(10), inside an airfilled rectangular cell of aspect ratio unity and pi spanwise openended distance. As a main feature, the average rate of the invariants of the velocity gradient tensor (Q(G), RG) has displayed the socalled
Copyright 2016 AIP Publishing. This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing.
Thu, 30 Mar 2017 10:44:50 GMT
http://hdl.handle.net/2117/103099
20170330T10:44:50Z
Dabbagh, Firas
Trias Miquel, Francesc Xavier
Gorobets, Andrei
Oliva Llena, Asensio
Smallscale dynamics is the spirit of turbulence physics. It implicates many attributes of flow topology evolution, coherent structures, hairpin vorticity dynamics, and mechanism of the kinetic energy cascade. In this work, several dynamical aspects of the smallscale motions have been numerically studied in a framework of RayleighBenard convection (RBC). To do so, direct numerical simulations have been carried out at two Rayleigh numbers Ra = 10(8) and 10(10), inside an airfilled rectangular cell of aspect ratio unity and pi spanwise openended distance. As a main feature, the average rate of the invariants of the velocity gradient tensor (Q(G), RG) has displayed the socalled

On the properties of discrete spatial filters for CFD
http://hdl.handle.net/2117/102792
On the properties of discrete spatial filters for CFD
Baez Vidal, Aleix; Lehmkuhl Barba, Oriol; Trias Miquel, Francesc Xavier; Pérez Segarra, Carlos David
The spatial filtering of variables in the context of Computational Fluid Dynamics (CFD) is a common practice. Most of the discrete filters used in CFD simulations are locally accurate models of continuous operators. However, when filters are adaptative, i.e. the filter width is not constant, or meshes are irregular, discrete filters sometimes break relevant global properties of the continuous models they are based on. For example, the principle of maxima and minima reduction or conservation are eventually infringed. In this paper, we analyze the properties of analytic continuous convolution filters and extract those we consider to define filtering. Then, we impose the accomplishment of these properties on explicit discrete filters by means of constraints. Three filters satisfying the derived conditions are deduced and compared to common differential discrete CFD filters on synthetic fields. Tests on the developed discrete filters show the fulfillment of the imposed properties. In particular, the problem of maxima and minima generation is resolved for physically relevant cases. The tests are conducted on the basis of the eigenvectors of graph Laplacian matrices of meshes. Thus, insight into the relations between filtering and oscillation growth on general meshes is provided. Further tests on singularity fields and on isentropic vortices have also been conducted to evaluate the performance of filters on basic CFD fields. Results confirm that imposing the proposed conditions makes discrete filters properties consistent with those of the continuous ones.
© 2016. This version is made available under the CCBYNCND 4.0 license http://creativecommons.org/licenses/byncnd/4.0/
Wed, 22 Mar 2017 12:10:06 GMT
http://hdl.handle.net/2117/102792
20170322T12:10:06Z
Baez Vidal, Aleix
Lehmkuhl Barba, Oriol
Trias Miquel, Francesc Xavier
Pérez Segarra, Carlos David
The spatial filtering of variables in the context of Computational Fluid Dynamics (CFD) is a common practice. Most of the discrete filters used in CFD simulations are locally accurate models of continuous operators. However, when filters are adaptative, i.e. the filter width is not constant, or meshes are irregular, discrete filters sometimes break relevant global properties of the continuous models they are based on. For example, the principle of maxima and minima reduction or conservation are eventually infringed. In this paper, we analyze the properties of analytic continuous convolution filters and extract those we consider to define filtering. Then, we impose the accomplishment of these properties on explicit discrete filters by means of constraints. Three filters satisfying the derived conditions are deduced and compared to common differential discrete CFD filters on synthetic fields. Tests on the developed discrete filters show the fulfillment of the imposed properties. In particular, the problem of maxima and minima generation is resolved for physically relevant cases. The tests are conducted on the basis of the eigenvectors of graph Laplacian matrices of meshes. Thus, insight into the relations between filtering and oscillation growth on general meshes is provided. Further tests on singularity fields and on isentropic vortices have also been conducted to evaluate the performance of filters on basic CFD fields. Results confirm that imposing the proposed conditions makes discrete filters properties consistent with those of the continuous ones.