SAM
https://sam.ensam.eu:443
The DSpace digital repository system captures, stores, indexes, preserves, and distributes digital research material.Tue, 27 Feb 2024 09:23:38 GMT2024-02-27T09:23:38ZAn a posteriori-implicit turbulent model with automatic dissipation adjustment for Large Eddy Simulation of compressible flows
http://hdl.handle.net/10985/18002
An a posteriori-implicit turbulent model with automatic dissipation adjustment for Large Eddy Simulation of compressible flows
NOGUEIRA, Xesús; RAMÍREZ, Luis; FERNÁNDEZ-FIDALGO, Javier; DELIGANT, Michael; KHELLADI, Sofiane; CHASSAING, Jean-Camille; NAVARRINA, Fermín
In this work we present an a posteriori high-order finite volume scheme for the computation of compressible turbulent flows. An automatic dissipation adjustment (ADA) method is combined with the a posteriori paradigm, in order to obtain an implicit subgrid scale model and preserve the stability of the numerical method. Thus, the numerical scheme is designed to increase the dissipation in the control volumes where the flow is under-resolved, and to decrease the dissipation in those cells where there is excessive dissipation. This is achieved by adding a multiplicative factor to the dissipative part of the numerical flux. In order to keep the stability of the numerical scheme, the a posteriori approach is used. It allows to increase the dissipation quickly in cells near shocks if required, ensuring the stability of the scheme. Some numerical tests are performed to highlight the accuracy and robustness of the proposed numerical scheme.
Wed, 01 Jan 2020 00:00:00 GMThttp://hdl.handle.net/10985/180022020-01-01T00:00:00ZNOGUEIRA, XesúsRAMÍREZ, LuisFERNÁNDEZ-FIDALGO, JavierDELIGANT, MichaelKHELLADI, SofianeCHASSAING, Jean-CamilleNAVARRINA, FermínIn this work we present an a posteriori high-order finite volume scheme for the computation of compressible turbulent flows. An automatic dissipation adjustment (ADA) method is combined with the a posteriori paradigm, in order to obtain an implicit subgrid scale model and preserve the stability of the numerical method. Thus, the numerical scheme is designed to increase the dissipation in the control volumes where the flow is under-resolved, and to decrease the dissipation in those cells where there is excessive dissipation. This is achieved by adding a multiplicative factor to the dissipative part of the numerical flux. In order to keep the stability of the numerical scheme, the a posteriori approach is used. It allows to increase the dissipation quickly in cells near shocks if required, ensuring the stability of the scheme. Some numerical tests are performed to highlight the accuracy and robustness of the proposed numerical scheme.New high-resolution-preserving sliding mesh techniques for higher-order finite volume schemes
http://hdl.handle.net/10985/17818
New high-resolution-preserving sliding mesh techniques for higher-order finite volume schemes
RAMÍREZ, Luis; FOULQUIÉ, Charles; NOGUEIRA, Xesús; KHELLADI, Sofiane; CHASSAING, Jean-Camille; COLOMINAS, Ignasi
This paper presents a new sliding mesh technique for the computation of unsteady viscous flows in the presence of rotating bodies. The compressible Euler and incompressible Navier–Stokes equations are solved using a higher-order (>2) finite volume method on unstructured grids. A sliding mesh approach is employed at the interface between computational grids in relative motion. In order to prevent loss of accuracy, two distinct families of higher-order sliding mesh interfaces are developed. These approaches fit naturally in a high-order finite volume framework. To this end, Moving Least Squares (MLS) approximants are used for the transmission of the information from one grid to another. A particular attention is paid for the study of the accuracy and conservation properties of the numerical scheme for static and rotating grids. The capabilities of the present solver to compute complex unsteady vortical flow motions created by rotating geometries are illustrated on a cross-flow configuration.
Thu, 01 Jan 2015 00:00:00 GMThttp://hdl.handle.net/10985/178182015-01-01T00:00:00ZRAMÍREZ, LuisFOULQUIÉ, CharlesNOGUEIRA, XesúsKHELLADI, SofianeCHASSAING, Jean-CamilleCOLOMINAS, IgnasiThis paper presents a new sliding mesh technique for the computation of unsteady viscous flows in the presence of rotating bodies. The compressible Euler and incompressible Navier–Stokes equations are solved using a higher-order (>2) finite volume method on unstructured grids. A sliding mesh approach is employed at the interface between computational grids in relative motion. In order to prevent loss of accuracy, two distinct families of higher-order sliding mesh interfaces are developed. These approaches fit naturally in a high-order finite volume framework. To this end, Moving Least Squares (MLS) approximants are used for the transmission of the information from one grid to another. A particular attention is paid for the study of the accuracy and conservation properties of the numerical scheme for static and rotating grids. The capabilities of the present solver to compute complex unsteady vortical flow motions created by rotating geometries are illustrated on a cross-flow configuration.A high-order density-based finite volume method for the computation of all-speed flows
http://hdl.handle.net/10985/17850
A high-order density-based finite volume method for the computation of all-speed flows
NOGUEIRA, Xesús; RAMÍREZ, Luis; KHELLADI, Sofiane; CHASSAING, Jean-Camille; COLOMINAS, Ignasi
In this paper we present a high-order density-based finite-volume framework for all-speed flows. The formulation is based on high-order variable reconstructions performed using Moving Least Squares approximations. In particular, we show that combining high-order discretization schemes with low-Mach fixes, it is possible to remove the grid dependency problem at low Mach numbers on both structured and unstructured grids. In order to maintain the accuracy and the robustness of the numerical method at transonic conditions, different procedures are proposed, based on the use of a selective limiting.
Fri, 01 Jan 2016 00:00:00 GMThttp://hdl.handle.net/10985/178502016-01-01T00:00:00ZNOGUEIRA, XesúsRAMÍREZ, LuisKHELLADI, SofianeCHASSAING, Jean-CamilleCOLOMINAS, IgnasiIn this paper we present a high-order density-based finite-volume framework for all-speed flows. The formulation is based on high-order variable reconstructions performed using Moving Least Squares approximations. In particular, we show that combining high-order discretization schemes with low-Mach fixes, it is possible to remove the grid dependency problem at low Mach numbers on both structured and unstructured grids. In order to maintain the accuracy and the robustness of the numerical method at transonic conditions, different procedures are proposed, based on the use of a selective limiting.Model-form and predictive uncertainty quantification in linear aeroelasticity
http://hdl.handle.net/10985/15562
Model-form and predictive uncertainty quantification in linear aeroelasticity
NITSCHKE, Christian T.; CINNELLA, Paola; LUCOR, Didier; CHASSAING, Jean-Camille
In this work, Bayesian techniques are employed to quantify model-form and predictive uncertainty in the linear behavior of an elastically mounted airfoil undergoing pitching and plunging motions. The Bayesian model averaging approach is used to construct an adjusted stochastic model from different model classes for time-harmonic incompressible flows. From a set of deterministic function approximations, we construct different stochastic models, whose uncertain coefficients are calibrated using Bayesian inference with regard to the critical flutter velocity. Results show substantial reductions in the predictive uncertainties of the critical flutter speed compared to non-calibrated stochastic simulations. In particular, it is shown that an efficient adjusted model can be derived by considering a possible bias in the random error term on the posterior predictive distributions of the flutter index.
Sun, 01 Jan 2017 00:00:00 GMThttp://hdl.handle.net/10985/155622017-01-01T00:00:00ZNITSCHKE, Christian T.CINNELLA, PaolaLUCOR, DidierCHASSAING, Jean-CamilleIn this work, Bayesian techniques are employed to quantify model-form and predictive uncertainty in the linear behavior of an elastically mounted airfoil undergoing pitching and plunging motions. The Bayesian model averaging approach is used to construct an adjusted stochastic model from different model classes for time-harmonic incompressible flows. From a set of deterministic function approximations, we construct different stochastic models, whose uncertain coefficients are calibrated using Bayesian inference with regard to the critical flutter velocity. Results show substantial reductions in the predictive uncertainties of the critical flutter speed compared to non-calibrated stochastic simulations. In particular, it is shown that an efficient adjusted model can be derived by considering a possible bias in the random error term on the posterior predictive distributions of the flutter index.3D global optimal forcing and response of the supersonic boundary layer
http://hdl.handle.net/10985/17986
3D global optimal forcing and response of the supersonic boundary layer
CHASSAING, Jean-Camille; SAGAUT, P.; BUGEAT, Benjamin; ROBINET, Jean-Christophe
3D optimal forcing and response of a 2D supersonic boundary layer are obtained by computing the largest singular value and the associated singular vectors of the global resolvent matrix. This approach allows to take into account both convective-type and component-type non-normalities responsible for the non-modal growth of perturbations in noise selective amplifier flows. It is moreover a fully non-parallel approach that does not require any particular assumptions on the baseflow. The numerical method is based on the explicit calculation of the Jacobian matrix proposed by Mettot et al. [1] for 2D perturbations. This strategy uses the numerical residual of the compressible Navier-Stokes equations imported from a finite-volume solver that is then linearised employing a finite difference method. Extension to 3D perturbations, which are expanded into modes of wave number, is here proposed by decomposing the Jacobian matrix according to the direction of the derivatives contained in its coefficients. Validation is performed on a Blasius boundary layer and a supersonic boundary layer, in comparison respectively to global and local results. Application of the method to a boundary layer at M = 4.5 recovers three regions of receptivity in the frequency-transverse wave number space. Finally, the energy growth of each optimal response is studied and discussed.
Tue, 01 Jan 2019 00:00:00 GMThttp://hdl.handle.net/10985/179862019-01-01T00:00:00ZCHASSAING, Jean-CamilleSAGAUT, P.BUGEAT, BenjaminROBINET, Jean-Christophe3D optimal forcing and response of a 2D supersonic boundary layer are obtained by computing the largest singular value and the associated singular vectors of the global resolvent matrix. This approach allows to take into account both convective-type and component-type non-normalities responsible for the non-modal growth of perturbations in noise selective amplifier flows. It is moreover a fully non-parallel approach that does not require any particular assumptions on the baseflow. The numerical method is based on the explicit calculation of the Jacobian matrix proposed by Mettot et al. [1] for 2D perturbations. This strategy uses the numerical residual of the compressible Navier-Stokes equations imported from a finite-volume solver that is then linearised employing a finite difference method. Extension to 3D perturbations, which are expanded into modes of wave number, is here proposed by decomposing the Jacobian matrix according to the direction of the derivatives contained in its coefficients. Validation is performed on a Blasius boundary layer and a supersonic boundary layer, in comparison respectively to global and local results. Application of the method to a boundary layer at M = 4.5 recovers three regions of receptivity in the frequency-transverse wave number space. Finally, the energy growth of each optimal response is studied and discussed.Low-frequency resolvent analysis of the laminar oblique shock wave/boundary layer interaction
http://hdl.handle.net/10985/22015
Low-frequency resolvent analysis of the laminar oblique shock wave/boundary layer interaction
BUGEAT, Benjamin; CHASSAING, Jean-Camille; SAGAUT, Pierre; ROBINET, Jean-Christophe
Resolvent analysis is used to study the low-frequency behaviour of the laminar oblique shock wave/boundary layer interaction (SWBLI). It is shown that the computed optimal gain, which can be seen as a transfer function of the system, follows a first-order low-pass filter equation, recovering the results of Touber & Sandham (J. Fluid Mech., vol. 671, 2011, pp. 417–465). This behaviour is understood as proceeding from the excitation of a single stable, steady global mode whose damping rate sets the time scale of the filter. Different Mach and Reynolds numbers are studied, covering different recirculation lengths L. This damping rate is found to scale as 1/L, leading to a constant Strouhal number StL as observed in the literature. It is associated with a breathing motion of the recirculation bubble. This analysis furthermore supports the idea that the low-frequency dynamics of the SWBLI is a forced dynamics, in which background perturbations continuously excite the flow. The investigation is then carried out for three-dimensional perturbations for which two regimes are identified. At low wavenumbers of the order of L, a modal mechanism
similar to that of two-dimensional perturbations is found and exhibits larger values of the optimal gain. At larger wavenumbers, of the order of the boundary layer thickness, the growth of streaks, which results from a non-modal mechanism, is detected. No interaction with the recirculation region is observed. Based on these results, the potential prevalence of three-dimensional effects in the low-frequency dynamics of the SWBLI is discussed.
Fri, 01 Apr 2022 00:00:00 GMThttp://hdl.handle.net/10985/220152022-04-01T00:00:00ZBUGEAT, BenjaminCHASSAING, Jean-CamilleSAGAUT, PierreROBINET, Jean-ChristopheResolvent analysis is used to study the low-frequency behaviour of the laminar oblique shock wave/boundary layer interaction (SWBLI). It is shown that the computed optimal gain, which can be seen as a transfer function of the system, follows a first-order low-pass filter equation, recovering the results of Touber & Sandham (J. Fluid Mech., vol. 671, 2011, pp. 417–465). This behaviour is understood as proceeding from the excitation of a single stable, steady global mode whose damping rate sets the time scale of the filter. Different Mach and Reynolds numbers are studied, covering different recirculation lengths L. This damping rate is found to scale as 1/L, leading to a constant Strouhal number StL as observed in the literature. It is associated with a breathing motion of the recirculation bubble. This analysis furthermore supports the idea that the low-frequency dynamics of the SWBLI is a forced dynamics, in which background perturbations continuously excite the flow. The investigation is then carried out for three-dimensional perturbations for which two regimes are identified. At low wavenumbers of the order of L, a modal mechanism
similar to that of two-dimensional perturbations is found and exhibits larger values of the optimal gain. At larger wavenumbers, of the order of the boundary layer thickness, the growth of streaks, which results from a non-modal mechanism, is detected. No interaction with the recirculation region is observed. Based on these results, the potential prevalence of three-dimensional effects in the low-frequency dynamics of the SWBLI is discussed.