SAM
https://sam.ensam.eu:443
The DSpace digital repository system captures, stores, indexes, preserves, and distributes digital research material.Fri, 02 Dec 2022 08:05:49 GMT2022-12-02T08:05:49ZSuccessive bifurcations in a fully three-dimensional open cavity flow
http://hdl.handle.net/10985/13014
Successive bifurcations in a fully three-dimensional open cavity flow
PICELLA, Francesco; LOISEAU, Jean-Christophe; LUSSEYRAN, F; ROBINET, Jean-Christophe; CHERUBINI, Stefania; PASTUR, L
The transition to unsteadiness of a three-dimensional open cavity flow is investigated using the joint application of direct numerical simulations and fully three-dimensional linear stability analyses, providing a clear understanding of the first two bifurcations occurring in the flow. The first bifurcation is characterized by the emergence of Taylor–Görtler-like vortices resulting from a centrifugal instability of the primary vortex core. Further increasing the Reynolds number eventually triggers self-sustained periodic oscillations of the flow in the vicinity of the spanwise end walls of the cavity. This secondary instability causes the emergence of a new set of Taylor–Görtler vortices experiencing a spanwise drift directed toward the spanwise end walls of the cavity. While a two-dimensional stability analysis would fail to capture this secondary instability due to the neglect of the lateral walls, it is the first time to our knowledge that this drifting of the vortices can be entirely characterized by a three-dimensional linear stability analysis of the flow. Good agreements with experimental observations and measurements strongly support our claim that the initial stages of the transition to turbulence of three-dimensional open cavity flows are solely governed by modal instabilities.
Mon, 01 Jan 2018 00:00:00 GMThttp://hdl.handle.net/10985/130142018-01-01T00:00:00ZPICELLA, FrancescoLOISEAU, Jean-ChristopheLUSSEYRAN, FROBINET, Jean-ChristopheCHERUBINI, StefaniaPASTUR, LThe transition to unsteadiness of a three-dimensional open cavity flow is investigated using the joint application of direct numerical simulations and fully three-dimensional linear stability analyses, providing a clear understanding of the first two bifurcations occurring in the flow. The first bifurcation is characterized by the emergence of Taylor–Görtler-like vortices resulting from a centrifugal instability of the primary vortex core. Further increasing the Reynolds number eventually triggers self-sustained periodic oscillations of the flow in the vicinity of the spanwise end walls of the cavity. This secondary instability causes the emergence of a new set of Taylor–Görtler vortices experiencing a spanwise drift directed toward the spanwise end walls of the cavity. While a two-dimensional stability analysis would fail to capture this secondary instability due to the neglect of the lateral walls, it is the first time to our knowledge that this drifting of the vortices can be entirely characterized by a three-dimensional linear stability analysis of the flow. Good agreements with experimental observations and measurements strongly support our claim that the initial stages of the transition to turbulence of three-dimensional open cavity flows are solely governed by modal instabilities.Dense-gas effects on compressible boundary-layer stability
http://hdl.handle.net/10985/18556
Dense-gas effects on compressible boundary-layer stability
GLOERFELT, Xavier; ROBINET, Jean-Christophe; SCIACOVELLI, Luca; CINNELLA, Paola; GRASSO, Francesco
A study of dense-gas effects on the stability of compressible boundary-layer flows is conducted. From the laminar similarity solution, the temperature variations are small due to the high specific heat of dense gases, leading to velocity profiles close to the incompressible ones. Concurrently, the complex thermodynamic properties of dense gases can lead to unconventional compressibility effects. In the subsonic regime, the Tollmien–Schlichting viscous mode is attenuated by compressibility effects and becomes preferentially skewed in line with the results based on the ideal-gas assumption. However, the absence of a generalized inflection point precludes the sustainability of the first mode by inviscid mechanisms. On the contrary, the viscous mode can be completely stable at supersonic speeds. At very high speeds, we have found instances of radiating supersonic instabilities with substantial amplification rates, i.e. waves that travel supersonically relative to the free-stream velocity. This acoustic mode has qualitatively similar features for various thermodynamic conditions and for different working fluids. This shows that the leading parameters governing the boundary-layer behaviour for the dense gas are the constant-pressure specific heat and, to a minor extent, the density-dependent viscosity. A satisfactory scaling of the mode characteristics is found to be proportional to the height of the layer near the wall that acts as a waveguide where acoustic waves may become trapped. This means that the supersonic mode has the same nature as Mack’s modes, even if its frequency for maximal amplification is greater. Direct numerical simulation accurately reproduces the development of the supersonic mode and emphasizes the radiation of the instability waves.
Wed, 01 Jan 2020 00:00:00 GMThttp://hdl.handle.net/10985/185562020-01-01T00:00:00ZGLOERFELT, XavierROBINET, Jean-ChristopheSCIACOVELLI, LucaCINNELLA, PaolaGRASSO, FrancescoA study of dense-gas effects on the stability of compressible boundary-layer flows is conducted. From the laminar similarity solution, the temperature variations are small due to the high specific heat of dense gases, leading to velocity profiles close to the incompressible ones. Concurrently, the complex thermodynamic properties of dense gases can lead to unconventional compressibility effects. In the subsonic regime, the Tollmien–Schlichting viscous mode is attenuated by compressibility effects and becomes preferentially skewed in line with the results based on the ideal-gas assumption. However, the absence of a generalized inflection point precludes the sustainability of the first mode by inviscid mechanisms. On the contrary, the viscous mode can be completely stable at supersonic speeds. At very high speeds, we have found instances of radiating supersonic instabilities with substantial amplification rates, i.e. waves that travel supersonically relative to the free-stream velocity. This acoustic mode has qualitatively similar features for various thermodynamic conditions and for different working fluids. This shows that the leading parameters governing the boundary-layer behaviour for the dense gas are the constant-pressure specific heat and, to a minor extent, the density-dependent viscosity. A satisfactory scaling of the mode characteristics is found to be proportional to the height of the layer near the wall that acts as a waveguide where acoustic waves may become trapped. This means that the supersonic mode has the same nature as Mack’s modes, even if its frequency for maximal amplification is greater. Direct numerical simulation accurately reproduces the development of the supersonic mode and emphasizes the radiation of the instability waves.Global Stability Analyses Unraveling Roughness-induced Transition Mechanisms
http://hdl.handle.net/10985/17848
Global Stability Analyses Unraveling Roughness-induced Transition Mechanisms
LOISEAU, Jean-Christophe; ROBINET, Jean-Christophe; CHERUBINI, Stefania; LERICHE, Emmanuel
The linear global instability and resulting transition to turbulence induced by a cylindrical roughness element of heighth and diameter d=3h immersed within an incompressible boundary layer flow along a flat plate is investigated using the joint application of direct numerical simulations and three-dimensional stability analyses. The configuration investigated is the same as the one investigated experimentally by Fransson et al. Base flow computations show that the roughness element induces a wake composed of a central low-speed region surrounded by a three-dimensional shear layer and a pair of low- and high-speed streaks on each side. Results from the global stability analyses highlight the unstable nature of the central low-speed region and its crucial importance in the laminar-turbulent transition process. For the set of parameters considered, it is able to sustain a varicose global instability for which the predicted critical Reynolds number is only 6% larger than the one reported in Ref. 10. A kinetic energy budget and wavemaker analysis revealed that this mode finds its root in the reversed flow region right downstream the roughness element and extracts most of its energy from the central low-speed region and streaks further downstream. Direct numerical simulations of the flow past this roughness element puts in the limelight the ability for this linear instability to give birth to hairpin vortices and thus trigger transition to turbulence.
Thu, 01 Jan 2015 00:00:00 GMThttp://hdl.handle.net/10985/178482015-01-01T00:00:00ZLOISEAU, Jean-ChristopheROBINET, Jean-ChristopheCHERUBINI, StefaniaLERICHE, EmmanuelThe linear global instability and resulting transition to turbulence induced by a cylindrical roughness element of heighth and diameter d=3h immersed within an incompressible boundary layer flow along a flat plate is investigated using the joint application of direct numerical simulations and three-dimensional stability analyses. The configuration investigated is the same as the one investigated experimentally by Fransson et al. Base flow computations show that the roughness element induces a wake composed of a central low-speed region surrounded by a three-dimensional shear layer and a pair of low- and high-speed streaks on each side. Results from the global stability analyses highlight the unstable nature of the central low-speed region and its crucial importance in the laminar-turbulent transition process. For the set of parameters considered, it is able to sustain a varicose global instability for which the predicted critical Reynolds number is only 6% larger than the one reported in Ref. 10. A kinetic energy budget and wavemaker analysis revealed that this mode finds its root in the reversed flow region right downstream the roughness element and extracts most of its energy from the central low-speed region and streaks further downstream. Direct numerical simulations of the flow past this roughness element puts in the limelight the ability for this linear instability to give birth to hairpin vortices and thus trigger transition to turbulence.Influence of the Shape on the Roughness-Induced Transition
http://hdl.handle.net/10985/17803
Influence of the Shape on the Roughness-Induced Transition
LOISEAU, Jean-Christophe; CHERUBINI, Stefania; ROBINET, Jean-Christophe; LERICHE, Emmanuel
lobal instability analysis of the three-dimensional flow past two rough- ness elements of different shape, namely a cylinder and a bump, is presented. In both cases, the eigenspectrum is made of modes characterised by a varicose symmetry and localised mostly in the zones of large base flow shear. The primary instabil- ity exhibited is the same in both cases and consists in an isolated unstable mode closely related to streaks local instability. For the cylinder however, a whole branch of modes is in addition destabilised as the Reynolds number is further increased.
Thu, 01 Jan 2015 00:00:00 GMThttp://hdl.handle.net/10985/178032015-01-01T00:00:00ZLOISEAU, Jean-ChristopheCHERUBINI, StefaniaROBINET, Jean-ChristopheLERICHE, Emmanuellobal instability analysis of the three-dimensional flow past two rough- ness elements of different shape, namely a cylinder and a bump, is presented. In both cases, the eigenspectrum is made of modes characterised by a varicose symmetry and localised mostly in the zones of large base flow shear. The primary instabil- ity exhibited is the same in both cases and consists in an isolated unstable mode closely related to streaks local instability. For the cylinder however, a whole branch of modes is in addition destabilised as the Reynolds number is further increased.Optimal perturbation for two-dimensional vortex systems: route to non-axisymmetric state
http://hdl.handle.net/10985/14216
Optimal perturbation for two-dimensional vortex systems: route to non-axisymmetric state
NAVROSE; JOHNSON, Author H. G.; BRION, Véronique; JACQUIN, Laurent; ROBINET, Jean-Christophe
We investigate perturbations that maximize the gain of disturbance energy in a two-dimensional isolated vortex and a counter-rotating vortex pair. The optimization is carried out using the method of Lagrange multipliers. For low initial energy of the perturbation ( ), the nonlinear optimal perturbation/gain is found to be the same as the linear optimal perturbation/gain. Beyond a certain threshold , the optimal perturbation/gain obtained from linear and nonlinear computations are different. There exists a range of for which the nonlinear optimal gain is higher than the linear optimal gain. For an isolated vortex, the higher value of nonlinear optimal gain is attributed to interaction among different azimuthal components, which is otherwise absent in a linearized system. Spiral dislocations are found in the nonlinear optimal perturbation at the radial location where the most dominant wavenumber changes. Long-time nonlinear evolution of linear and nonlinear optimal perturbations is studied. The evolution shows that, after the initial increment of perturbation energy, the vortex attains a quasi-steady state where the mean perturbation energy decreases on a slow time scale. The quasi-steady vortex state is non-axisymmetric and its shape depends on the initial perturbation. It is observed that the lifetime of a quasi-steady vortex state obtained using the nonlinear optimal perturbation is longer than that obtained using the linear optimal perturbation. For a counter-rotating vortex pair, the mechanism that maximizes the energy gain is found to be similar to that of the isolated vortex. Within the linear framework, the optimal perturbation for a vortex pair can be either symmetric or antisymmetric, whereas the structure of the nonlinear optimal perturbation, beyond the threshold, is always asymmetric. No quasi-steady state for a counter-rotating vortex pair is observed.
Mon, 01 Jan 2018 00:00:00 GMThttp://hdl.handle.net/10985/142162018-01-01T00:00:00ZNAVROSEJOHNSON, Author H. G.BRION, VéroniqueJACQUIN, LaurentROBINET, Jean-ChristopheWe investigate perturbations that maximize the gain of disturbance energy in a two-dimensional isolated vortex and a counter-rotating vortex pair. The optimization is carried out using the method of Lagrange multipliers. For low initial energy of the perturbation ( ), the nonlinear optimal perturbation/gain is found to be the same as the linear optimal perturbation/gain. Beyond a certain threshold , the optimal perturbation/gain obtained from linear and nonlinear computations are different. There exists a range of for which the nonlinear optimal gain is higher than the linear optimal gain. For an isolated vortex, the higher value of nonlinear optimal gain is attributed to interaction among different azimuthal components, which is otherwise absent in a linearized system. Spiral dislocations are found in the nonlinear optimal perturbation at the radial location where the most dominant wavenumber changes. Long-time nonlinear evolution of linear and nonlinear optimal perturbations is studied. The evolution shows that, after the initial increment of perturbation energy, the vortex attains a quasi-steady state where the mean perturbation energy decreases on a slow time scale. The quasi-steady vortex state is non-axisymmetric and its shape depends on the initial perturbation. It is observed that the lifetime of a quasi-steady vortex state obtained using the nonlinear optimal perturbation is longer than that obtained using the linear optimal perturbation. For a counter-rotating vortex pair, the mechanism that maximizes the energy gain is found to be similar to that of the isolated vortex. Within the linear framework, the optimal perturbation for a vortex pair can be either symmetric or antisymmetric, whereas the structure of the nonlinear optimal perturbation, beyond the threshold, is always asymmetric. No quasi-steady state for a counter-rotating vortex pair is observed.Multiple-correction hybrid k -exact schemes for high-order compressible RANS-LES simulations on fully unstructured grids
http://hdl.handle.net/10985/15515
Multiple-correction hybrid k -exact schemes for high-order compressible RANS-LES simulations on fully unstructured grids; Multiple-correction hybrid k -exact schemes for high-order compressible RANS-LES simulations on fully unstructured grids
PONT, Grégoire; PONT, Grégoire; BRENNER, Pierre; BRENNER, Pierre; CINNELLA, Paola; CINNELLA, Paola; MAUGARS, Bruno; MAUGARS, Bruno; ROBINET, Jean-Christophe; ROBINET, Jean-Christophe
A Godunov's type unstructured finite volume method suitable for highly compressible turbulent scale-resolving simulations around complex geometries is constructed by using a successive correction technique. First, a family of k-exact Godunov schemes is developed by recursively correcting the truncation error of the piecewise polynomial representation of the primitive variables. The keystone of the proposed approach is a quasi-Green gradient operator which ensures consistency on general meshes. In addition, a high-order single-point quadrature formula, based on high-order approximations of the successive derivatives of the solution, is developed for flux integration along cell faces. The proposed family of schemes is compact in the algorithmic sense, since it only involves communications between direct neighbors of the mesh cells. The numerical properties of the schemes up to fifth-order are investigated, with focus on their resolvability in terms of number of mesh points required to resolve a given wavelength accurately. Afterwards, in the aim of achieving the best possible trade-off between accuracy, computational cost and robustness in view of industrial flow computations, we focus more specifically on the third-order accurate scheme of the family, and modify locally its numerical flux in order to reduce the amount of numerical dissipation in vortex-dominated regions. This is achieved by switching from the upwind scheme, mostly applied in highly compressible regions, to a fourth-order centered one in vortex-dominated regions. An analytical switch function based on the local grid Reynolds number is adopted in order to warrant numerical stability of the recentering process. Numerical applications demonstrate the accuracy and robustness of the proposed methodology for compressible scale-resolving computations. In particular, supersonic RANS/LES computations of the flow over a cavity are presented to show the capability of the scheme to predict flows with shocks, vortical structures and complex geometries.; A Godunov's type unstructured finite volume method suitable for highly compressible turbulent scale-resolving simulations around complex geometries is constructed by using a successive correction technique. First, a family of k-exact Godunov schemes is developed by recursively correcting the truncation error of the piecewise polynomial representation of the primitive variables. The keystone of the proposed approach is a quasi-Green gradient operator which ensures consistency on general meshes. In addition, a high-order single-point quadrature formula, based on high-order approximations of the successive derivatives of the solution, is developed for flux integration along cell faces. The proposed family of schemes is compact in the algorithmic sense, since it only involves communications between direct neighbors of the mesh cells. The numerical properties of the schemes up to fifth-order are investigated, with focus on their resolvability in terms of number of mesh points required to resolve a given wavelength accurately. Afterwards, in the aim of achieving the best possible trade-off between accuracy, computational cost and robustness in view of industrial flow computations, we focus more specifically on the third-order accurate scheme of the family, and modify locally its numerical flux in order to reduce the amount of numerical dissipation in vortex-dominated regions. This is achieved by switching from the upwind scheme, mostly applied in highly compressible regions, to a fourth-order centered one in vortex-dominated regions. An analytical switch function based on the local grid Reynolds number is adopted in order to warrant numerical stability of the recentering process. Numerical applications demonstrate the accuracy and robustness of the proposed methodology for compressible scale-resolving computations. In particular, supersonic RANS/LES computations of the flow over a cavity are presented to show the capability of the scheme to predict flows with shocks, vortical structures and complex geometries.
Sun, 01 Jan 2017 00:00:00 GMThttp://hdl.handle.net/10985/155152017-01-01T00:00:00ZPONT, GrégoirePONT, GrégoireBRENNER, PierreBRENNER, PierreCINNELLA, PaolaCINNELLA, PaolaMAUGARS, BrunoMAUGARS, BrunoROBINET, Jean-ChristopheROBINET, Jean-ChristopheA Godunov's type unstructured finite volume method suitable for highly compressible turbulent scale-resolving simulations around complex geometries is constructed by using a successive correction technique. First, a family of k-exact Godunov schemes is developed by recursively correcting the truncation error of the piecewise polynomial representation of the primitive variables. The keystone of the proposed approach is a quasi-Green gradient operator which ensures consistency on general meshes. In addition, a high-order single-point quadrature formula, based on high-order approximations of the successive derivatives of the solution, is developed for flux integration along cell faces. The proposed family of schemes is compact in the algorithmic sense, since it only involves communications between direct neighbors of the mesh cells. The numerical properties of the schemes up to fifth-order are investigated, with focus on their resolvability in terms of number of mesh points required to resolve a given wavelength accurately. Afterwards, in the aim of achieving the best possible trade-off between accuracy, computational cost and robustness in view of industrial flow computations, we focus more specifically on the third-order accurate scheme of the family, and modify locally its numerical flux in order to reduce the amount of numerical dissipation in vortex-dominated regions. This is achieved by switching from the upwind scheme, mostly applied in highly compressible regions, to a fourth-order centered one in vortex-dominated regions. An analytical switch function based on the local grid Reynolds number is adopted in order to warrant numerical stability of the recentering process. Numerical applications demonstrate the accuracy and robustness of the proposed methodology for compressible scale-resolving computations. In particular, supersonic RANS/LES computations of the flow over a cavity are presented to show the capability of the scheme to predict flows with shocks, vortical structures and complex geometries.
A Godunov's type unstructured finite volume method suitable for highly compressible turbulent scale-resolving simulations around complex geometries is constructed by using a successive correction technique. First, a family of k-exact Godunov schemes is developed by recursively correcting the truncation error of the piecewise polynomial representation of the primitive variables. The keystone of the proposed approach is a quasi-Green gradient operator which ensures consistency on general meshes. In addition, a high-order single-point quadrature formula, based on high-order approximations of the successive derivatives of the solution, is developed for flux integration along cell faces. The proposed family of schemes is compact in the algorithmic sense, since it only involves communications between direct neighbors of the mesh cells. The numerical properties of the schemes up to fifth-order are investigated, with focus on their resolvability in terms of number of mesh points required to resolve a given wavelength accurately. Afterwards, in the aim of achieving the best possible trade-off between accuracy, computational cost and robustness in view of industrial flow computations, we focus more specifically on the third-order accurate scheme of the family, and modify locally its numerical flux in order to reduce the amount of numerical dissipation in vortex-dominated regions. This is achieved by switching from the upwind scheme, mostly applied in highly compressible regions, to a fourth-order centered one in vortex-dominated regions. An analytical switch function based on the local grid Reynolds number is adopted in order to warrant numerical stability of the recentering process. Numerical applications demonstrate the accuracy and robustness of the proposed methodology for compressible scale-resolving computations. In particular, supersonic RANS/LES computations of the flow over a cavity are presented to show the capability of the scheme to predict flows with shocks, vortical structures and complex geometries.Optimal perturbation for two-dimensional vortex systems: route to non-axisymmetric state
http://hdl.handle.net/10985/17877
Optimal perturbation for two-dimensional vortex systems: route to non-axisymmetric state
NAVROSE; JOHNSON, Author H. G.; BRION, Véronique; JACQUIN, Laurent; ROBINET, Jean-Christophe
We investigate perturbations that maximize the gain of disturbance energy in a two-dimensional isolated vortex and a counter-rotating vortex pair. The optimization is carried out using the method of Lagrange multipliers. For low initial energy of the perturbation (E.0/), the nonlinear optimal perturbation/gain is found to be the same as the linear optimal perturbation/gain. Beyond a certain threshold E.0/, the optimal perturbation/gain obtained from linear and nonlinear computations are different. There exists a range of E.0/ for which the nonlinear optimal gain is higher than the linear optimal gain. For an isolated vortex, the higher value of nonlinear optimal gain is attributed to interaction among different azimuthal components, which is otherwise absent in a linearized system. Spiral dislocations are found in the nonlinear optimal perturbation at the radial location where the most dominant wavenumber changes. Long-time nonlinear evolution of linear and nonlinear optimal perturbations is studied. The evolution shows that, after the initial increment of perturbation energy, the vortex attains a quasi-steady state where the mean perturbation energy decreases on a slow time scale. The quasi-steady vortex state is non-axisymmetric and its shape depends on the initial perturbation. It is observed that the lifetime of a quasi-steady vortex state obtained using the nonlinear optimal perturbation is longer than that obtained using the linear optimal perturbation. For a counter-rotating vortex pair, the mechanism that maximizes the energy gain is found to be similar to that of the isolated vortex. Within the linear framework, the optimal perturbation for a vortex pair can be either symmetric or antisymmetric, whereas the structure of the nonlinear optimal perturbation, beyond the threshold E.0/, is always asymmetric. No quasi-steady state for a counter-rotating vortex pair is observed.
Mon, 01 Jan 2018 00:00:00 GMThttp://hdl.handle.net/10985/178772018-01-01T00:00:00ZNAVROSEJOHNSON, Author H. G.BRION, VéroniqueJACQUIN, LaurentROBINET, Jean-ChristopheWe investigate perturbations that maximize the gain of disturbance energy in a two-dimensional isolated vortex and a counter-rotating vortex pair. The optimization is carried out using the method of Lagrange multipliers. For low initial energy of the perturbation (E.0/), the nonlinear optimal perturbation/gain is found to be the same as the linear optimal perturbation/gain. Beyond a certain threshold E.0/, the optimal perturbation/gain obtained from linear and nonlinear computations are different. There exists a range of E.0/ for which the nonlinear optimal gain is higher than the linear optimal gain. For an isolated vortex, the higher value of nonlinear optimal gain is attributed to interaction among different azimuthal components, which is otherwise absent in a linearized system. Spiral dislocations are found in the nonlinear optimal perturbation at the radial location where the most dominant wavenumber changes. Long-time nonlinear evolution of linear and nonlinear optimal perturbations is studied. The evolution shows that, after the initial increment of perturbation energy, the vortex attains a quasi-steady state where the mean perturbation energy decreases on a slow time scale. The quasi-steady vortex state is non-axisymmetric and its shape depends on the initial perturbation. It is observed that the lifetime of a quasi-steady vortex state obtained using the nonlinear optimal perturbation is longer than that obtained using the linear optimal perturbation. For a counter-rotating vortex pair, the mechanism that maximizes the energy gain is found to be similar to that of the isolated vortex. Within the linear framework, the optimal perturbation for a vortex pair can be either symmetric or antisymmetric, whereas the structure of the nonlinear optimal perturbation, beyond the threshold E.0/, is always asymmetric. No quasi-steady state for a counter-rotating vortex pair is observed.Bifurcation analysis and frequency prediction in shear-driven cavity flow
http://hdl.handle.net/10985/18039
Bifurcation analysis and frequency prediction in shear-driven cavity flow
BENGANA, Y.; LOISEAU, Jean-Christophe; ROBINET, Jean-Christophe; TUCKERMAN, L. S.
A comprehensive study of the two-dimensional incompressible shear-driven flow in an open square cavity is carried out. Two successive bifurcations lead to two limit cycles with different frequencies and different numbers of structures which propagate along the top of the cavity and circulate in its interior. A branch of quasi-periodic states produced by secondary Hopf bifurcations transfers the stability from one limit cycle to the other. A full analysis of this scenario is obtained by means of nonlinear simulations, linear stability analysis and Floquet analysis. We characterize the temporal behaviour of the limit cycles and quasi-periodic state via Fourier transforms and their spatial behaviour via the Hilbert transform. We address the relevance of linearization about the mean flow. Although here the nonlinear frequencies are not very far from those obtained by linearization about the base flow, the difference is substantially reduced when eigenvalues are obtained instead from linearization about the mean and in addition, the corresponding growth rate is small, a combination of properties called RZIF (real zero imaginary frequency). Moreover growth rates obtained by linearization about the mean of one limit cycle are correlated with relative stability to the other limit cycle. Finally, we show that the frequencies of the successive modes are separated by a constant increment.
Tue, 01 Jan 2019 00:00:00 GMThttp://hdl.handle.net/10985/180392019-01-01T00:00:00ZBENGANA, Y.LOISEAU, Jean-ChristopheROBINET, Jean-ChristopheTUCKERMAN, L. S.A comprehensive study of the two-dimensional incompressible shear-driven flow in an open square cavity is carried out. Two successive bifurcations lead to two limit cycles with different frequencies and different numbers of structures which propagate along the top of the cavity and circulate in its interior. A branch of quasi-periodic states produced by secondary Hopf bifurcations transfers the stability from one limit cycle to the other. A full analysis of this scenario is obtained by means of nonlinear simulations, linear stability analysis and Floquet analysis. We characterize the temporal behaviour of the limit cycles and quasi-periodic state via Fourier transforms and their spatial behaviour via the Hilbert transform. We address the relevance of linearization about the mean flow. Although here the nonlinear frequencies are not very far from those obtained by linearization about the base flow, the difference is substantially reduced when eigenvalues are obtained instead from linearization about the mean and in addition, the corresponding growth rate is small, a combination of properties called RZIF (real zero imaginary frequency). Moreover growth rates obtained by linearization about the mean of one limit cycle are correlated with relative stability to the other limit cycle. Finally, we show that the frequencies of the successive modes are separated by a constant increment.Large scale dynamics of a high Reynolds number axisymmetric separating/reattaching flow
http://hdl.handle.net/10985/17990
Large scale dynamics of a high Reynolds number axisymmetric separating/reattaching flow
PAIN, R; WEISS, P-E; DECK, S; ROBINET, Jean-Christophe
A numerical study is conducted to unveil the large scale dynamics of a high Reynolds number axisymmetric separating/reattaching flow at M∞ = 0.7. The numerical simulation allows us to acquire a high rate sampled unsteady volumetric dataset. This huge amount of spatial and temporal information is exploited in the Fourier space to visualize for the first time in physical space and at such a high Reynolds number (ReD = 1.2 × 106) the statistical signature of the helical structure related to the antisymmetric mode (m = 1) at StD = 0.18. The main hydrodynamic mechanisms are identified through the spatial distribution of the most energetic frequencies, i.e., StD = 0.18 and StD ≥ 3.0 corresponding to the vortex-shedding and Kelvin-Helmholtz instability phenomena, respectively. In particular, the dynamics related to the dimensionless shedding frequency is shown to become dominant for 0.35 ≤ x/D ≤ 0.75 in the whole radial direction as it passes through the shear layer. The spatial distribution of the coherence function for the most significant modes as well as a three-dimensional Fourier decomposition suggests the global features of the flow mechanisms. More specifically, the novelty of this study lies in the evidence of the flow dynamics through the use of cross-correlation maps plotted with a frequency selection guided by the characteristic Strouhal number formerly identified in a local manner in the flow field or at the wall. Moreover and for the first time, the understanding of the scales at stake is supported both by a Fourier analysis and a dynamic mode decomposition in the complete three-dimensional space surrounding the afterbody zone.
Tue, 01 Jan 2019 00:00:00 GMThttp://hdl.handle.net/10985/179902019-01-01T00:00:00ZPAIN, RWEISS, P-EDECK, SROBINET, Jean-ChristopheA numerical study is conducted to unveil the large scale dynamics of a high Reynolds number axisymmetric separating/reattaching flow at M∞ = 0.7. The numerical simulation allows us to acquire a high rate sampled unsteady volumetric dataset. This huge amount of spatial and temporal information is exploited in the Fourier space to visualize for the first time in physical space and at such a high Reynolds number (ReD = 1.2 × 106) the statistical signature of the helical structure related to the antisymmetric mode (m = 1) at StD = 0.18. The main hydrodynamic mechanisms are identified through the spatial distribution of the most energetic frequencies, i.e., StD = 0.18 and StD ≥ 3.0 corresponding to the vortex-shedding and Kelvin-Helmholtz instability phenomena, respectively. In particular, the dynamics related to the dimensionless shedding frequency is shown to become dominant for 0.35 ≤ x/D ≤ 0.75 in the whole radial direction as it passes through the shear layer. The spatial distribution of the coherence function for the most significant modes as well as a three-dimensional Fourier decomposition suggests the global features of the flow mechanisms. More specifically, the novelty of this study lies in the evidence of the flow dynamics through the use of cross-correlation maps plotted with a frequency selection guided by the characteristic Strouhal number formerly identified in a local manner in the flow field or at the wall. Moreover and for the first time, the understanding of the scales at stake is supported both by a Fourier analysis and a dynamic mode decomposition in the complete three-dimensional space surrounding the afterbody zone.Numerical investigation of three-dimensional partial cavitation in a Venturi geometry
http://hdl.handle.net/10985/20502
Numerical investigation of three-dimensional partial cavitation in a Venturi geometry
GOUIN, Camille; JUNQUEIRA-JUNIOR, Carlos; GONCALVES DA SILVA, Eric; ROBINET, Jean-Christophe
Sheet cavitation appears in many hydraulic applications and can lead to technical issues. Some fundamental outcomes, such as, the complex topology of 3-dimensional cavitation pockets and their associated dynamics need to be carefully visited. In the paper, the dynamics of partial cavitation developing in a 3D Venturi geometry and the interaction with sidewalls are numerically investigated. The simulations are performed using a one-fluid compressible Reynolds-averaged Navier–Stokes solver associated with a nonlinear turbulence model and a void ratio transport equation model. A detailed analysis of this cavitating flow is carried out using innovative tools, such as, spectral proper orthogonal decompositions. Particular attention is paid in the study of 3D effects by comparing the numerical results obtained with sidewalls and periodic conditions. A three-dimensional dynamics of the sheet cavitation, unrelated to the presence of sidewalls, is identified and discussed.
Fri, 01 Jan 2021 00:00:00 GMThttp://hdl.handle.net/10985/205022021-01-01T00:00:00ZGOUIN, CamilleJUNQUEIRA-JUNIOR, CarlosGONCALVES DA SILVA, EricROBINET, Jean-ChristopheSheet cavitation appears in many hydraulic applications and can lead to technical issues. Some fundamental outcomes, such as, the complex topology of 3-dimensional cavitation pockets and their associated dynamics need to be carefully visited. In the paper, the dynamics of partial cavitation developing in a 3D Venturi geometry and the interaction with sidewalls are numerically investigated. The simulations are performed using a one-fluid compressible Reynolds-averaged Navier–Stokes solver associated with a nonlinear turbulence model and a void ratio transport equation model. A detailed analysis of this cavitating flow is carried out using innovative tools, such as, spectral proper orthogonal decompositions. Particular attention is paid in the study of 3D effects by comparing the numerical results obtained with sidewalls and periodic conditions. A three-dimensional dynamics of the sheet cavitation, unrelated to the presence of sidewalls, is identified and discussed.