SAM
https://sam.ensam.eu:443
The DSpace digital repository system captures, stores, indexes, preserves, and distributes digital research material.Wed, 28 Feb 2024 15:28:31 GMT2024-02-28T15:28:31ZAn overview of the proper generalized decomposition with applications in computational rheology
http://hdl.handle.net/10985/8473
An overview of the proper generalized decomposition with applications in computational rheology
CHINESTA, Francisco; AMMAR, Amine; LEYGUE, Adrien; KEUNINGS, Roland
We review the foundations and applications of the proper generalized decomposition (PGD), a powerful model reduction technique that computes a priori by means of successive enrichment a separated representation of the unknown field. The computational complexity of the PGD scales linearly with the dimension of the space wherein the model is defined, which is in marked contrast with the exponential scaling of standard grid-based methods. First introduced in the context of computational rheology by Ammar et al. [3] and [4], the PGD has since been further developed and applied in a variety of applications ranging from the solution of the Schrödinger equation of quantum mechanics to the analysis of laminate composites. In this paper, we illustrate the use of the PGD in four problem categories related to computational rheology: (i) the direct solution of the Fokker-Planck equation for complex fluids in configuration spaces of high dimension, (ii) the development of very efficient non-incremental algorithms for transient problems, (iii) the fully three-dimensional solution of problems defined in degenerate plate or shell-like domains often encountered in polymer processing or composites manufacturing, and finally (iv) the solution of multidimensional parametric models obtained by introducing various sources of problem variability as additional coordinates.
Sat, 01 Jan 2011 00:00:00 GMThttp://hdl.handle.net/10985/84732011-01-01T00:00:00ZCHINESTA, FranciscoAMMAR, AmineLEYGUE, AdrienKEUNINGS, RolandWe review the foundations and applications of the proper generalized decomposition (PGD), a powerful model reduction technique that computes a priori by means of successive enrichment a separated representation of the unknown field. The computational complexity of the PGD scales linearly with the dimension of the space wherein the model is defined, which is in marked contrast with the exponential scaling of standard grid-based methods. First introduced in the context of computational rheology by Ammar et al. [3] and [4], the PGD has since been further developed and applied in a variety of applications ranging from the solution of the Schrödinger equation of quantum mechanics to the analysis of laminate composites. In this paper, we illustrate the use of the PGD in four problem categories related to computational rheology: (i) the direct solution of the Fokker-Planck equation for complex fluids in configuration spaces of high dimension, (ii) the development of very efficient non-incremental algorithms for transient problems, (iii) the fully three-dimensional solution of problems defined in degenerate plate or shell-like domains often encountered in polymer processing or composites manufacturing, and finally (iv) the solution of multidimensional parametric models obtained by introducing various sources of problem variability as additional coordinates.On the multi-scale description of micro-structured fluids composed of aggregating rods
http://hdl.handle.net/10985/17967
On the multi-scale description of micro-structured fluids composed of aggregating rods
PEREZ, Marta; SCHEUER, Adrien; ABISSET-CHAVANNE, Emmanuelle; AMMAR, Amine; CHINESTA, Francisco; KEUNINGS, Roland
When addressing the flow of concentrated suspensions composed of rods, dense clusters are observed. Thus, the adequate modelling and simulation of such a flow requires addressing the kinematics of these dense clusters and their impact on the flow in which they are immersed. In a former work, we addressed a first modelling framework of these clusters, assumed so dense that they were considered rigid and their kinematics (flow-induced rotation) were totally defined by a symmetric tensor c with unit trace representing the cluster conformation. Then, the rigid nature of the clusters was relaxed, assuming them deformable, and a model giving the evolution of both the cluster shape and its microstructural orientation descriptor (the so-called shape and orientation tensors) was proposed. This paper compares the predictions coming from those models with finer-scale discrete simulations inspired from molecular dynamics modelling.
Tue, 01 Jan 2019 00:00:00 GMThttp://hdl.handle.net/10985/179672019-01-01T00:00:00ZPEREZ, MartaSCHEUER, AdrienABISSET-CHAVANNE, EmmanuelleAMMAR, AmineCHINESTA, FranciscoKEUNINGS, RolandWhen addressing the flow of concentrated suspensions composed of rods, dense clusters are observed. Thus, the adequate modelling and simulation of such a flow requires addressing the kinematics of these dense clusters and their impact on the flow in which they are immersed. In a former work, we addressed a first modelling framework of these clusters, assumed so dense that they were considered rigid and their kinematics (flow-induced rotation) were totally defined by a symmetric tensor c with unit trace representing the cluster conformation. Then, the rigid nature of the clusters was relaxed, assuming them deformable, and a model giving the evolution of both the cluster shape and its microstructural orientation descriptor (the so-called shape and orientation tensors) was proposed. This paper compares the predictions coming from those models with finer-scale discrete simulations inspired from molecular dynamics modelling.Incremental dynamic mode decomposition: A reduced-model learner operating at the low-data limit
http://hdl.handle.net/10985/18539
Incremental dynamic mode decomposition: A reduced-model learner operating at the low-data limit
REILLE, Agathe; HASCOET, Nicolas; GHNATIOS, Chady; AMMAR, Amine; CUETO, Elias; DUVAL, Jean-Louis; CHINESTA, Francisco; KEUNINGS, Roland
The present work aims at proposing a new methodology for learning reduced models from a small amount of data. It is based on the fact that discrete models, or their transfer function counterparts, have a low rank and then they can be expressed very efficiently using few terms of a tensor decomposition. An efficient procedure is proposed as well as a way for extending it to nonlinear settings while keeping limited the impact of data noise. The proposed methodology is then validated by considering a nonlinear elastic problem and constructing the model relating tractions and displacements at the observation points.
Tue, 01 Jan 2019 00:00:00 GMThttp://hdl.handle.net/10985/185392019-01-01T00:00:00ZREILLE, AgatheHASCOET, NicolasGHNATIOS, ChadyAMMAR, AmineCUETO, EliasDUVAL, Jean-LouisCHINESTA, FranciscoKEUNINGS, RolandThe present work aims at proposing a new methodology for learning reduced models from a small amount of data. It is based on the fact that discrete models, or their transfer function counterparts, have a low rank and then they can be expressed very efficiently using few terms of a tensor decomposition. An efficient procedure is proposed as well as a way for extending it to nonlinear settings while keeping limited the impact of data noise. The proposed methodology is then validated by considering a nonlinear elastic problem and constructing the model relating tractions and displacements at the observation points.Modelling the effect of particle inertia on the orientation kinematics of fibres and spheroids immersed in a simple shear flow
http://hdl.handle.net/10985/18388
Modelling the effect of particle inertia on the orientation kinematics of fibres and spheroids immersed in a simple shear flow
SCHEUER, Adrien; GRÉGOIRE, Guillaume; ABISSET-CHAVANNE, Emmanuelle; CHINESTA, Francisco; KEUNINGS, Roland
Simulations of flows containing non-spherical particles (fibres or ellipsoids) rely on the knowledge of the equation governing the particle motion in the flow. Most models used nowadays are based on the pioneering work of Jeffery (1922), who obtained an equation for the motion of an ellipsoidal particle immersed in a Newtonian fluid, despite the fact that this model relies on strong assumptions: negligible inertia, unconfined flow, dilute regime, flow unperturbed by the presence of the suspended particle, etc. In this work, we propose a dumbbell-based model aimed to describe the motion of an inertial fibre or ellipsoid suspended in a Newtonian fluid. We then use this model to study the orientation kinematics of such particle in a linear shear flow and compare it to the inertialess case. In the case of fibres, we observe the appearance of periodic orbits (whereas inertialess fibres just align in the flow field). For spheroids, our model predicts an orbit drift towards the flow-gradient plane, either gradually (slight inertia) or by first rotating around a moving oblique axis (heavy particles). Multi-Particle Collision Dynamics (MPCD) simulations were carried out to assess the model predictions in the case of inertial fibres and revealed similar behaviours.
Wed, 01 Jan 2020 00:00:00 GMThttp://hdl.handle.net/10985/183882020-01-01T00:00:00ZSCHEUER, AdrienGRÉGOIRE, GuillaumeABISSET-CHAVANNE, EmmanuelleCHINESTA, FranciscoKEUNINGS, RolandSimulations of flows containing non-spherical particles (fibres or ellipsoids) rely on the knowledge of the equation governing the particle motion in the flow. Most models used nowadays are based on the pioneering work of Jeffery (1922), who obtained an equation for the motion of an ellipsoidal particle immersed in a Newtonian fluid, despite the fact that this model relies on strong assumptions: negligible inertia, unconfined flow, dilute regime, flow unperturbed by the presence of the suspended particle, etc. In this work, we propose a dumbbell-based model aimed to describe the motion of an inertial fibre or ellipsoid suspended in a Newtonian fluid. We then use this model to study the orientation kinematics of such particle in a linear shear flow and compare it to the inertialess case. In the case of fibres, we observe the appearance of periodic orbits (whereas inertialess fibres just align in the flow field). For spheroids, our model predicts an orbit drift towards the flow-gradient plane, either gradually (slight inertia) or by first rotating around a moving oblique axis (heavy particles). Multi-Particle Collision Dynamics (MPCD) simulations were carried out to assess the model predictions in the case of inertial fibres and revealed similar behaviours.From dilute to entangled fibre suspensions involved in the flow of reinforced polymers: A unified framework
http://hdl.handle.net/10985/12434
From dilute to entangled fibre suspensions involved in the flow of reinforced polymers: A unified framework
PEREZ, Marta; GUEVELOU, S; ABISSET-CHAVANNE, Emmanuelle; CHINESTA, Francisco; KEUNINGS, Roland
Most suspension descriptions nowadays employed are based on Jeffery model and some of its phenomenological adaptations that do not take into account the possible existence of a relative velocity between the fibres and the suspending fluid when the fibre interactions increase. It is expected that at very low density of contacts, as predicted by standard suspension models, fibres move with the suspending fluid velocity. When the density of fibre interactions becomes extremely high and a percolated network of fibre contacts is established within the suspension, fibres cannot move anymore and then the fluid flows throughout the rigid or moderately deformable entangled fibre skeleton, like a fluid flowing through a porous medium. In between these two limit cases, one could expect that fibres move but with a velocity lower than the one of the suspending fluid. Thus, two contributions are expected, one coming from standard suspension theory in which fibres and fluid move with the same velocity, and the other resulting in a Darcy contribution consisting of the relative fibre/fluid velocity. In this paper, we elaborate a general model able to adapt continuously to all these flow regimes.
Sun, 01 Jan 2017 00:00:00 GMThttp://hdl.handle.net/10985/124342017-01-01T00:00:00ZPEREZ, MartaGUEVELOU, SABISSET-CHAVANNE, EmmanuelleCHINESTA, FranciscoKEUNINGS, RolandMost suspension descriptions nowadays employed are based on Jeffery model and some of its phenomenological adaptations that do not take into account the possible existence of a relative velocity between the fibres and the suspending fluid when the fibre interactions increase. It is expected that at very low density of contacts, as predicted by standard suspension models, fibres move with the suspending fluid velocity. When the density of fibre interactions becomes extremely high and a percolated network of fibre contacts is established within the suspension, fibres cannot move anymore and then the fluid flows throughout the rigid or moderately deformable entangled fibre skeleton, like a fluid flowing through a porous medium. In between these two limit cases, one could expect that fibres move but with a velocity lower than the one of the suspending fluid. Thus, two contributions are expected, one coming from standard suspension theory in which fibres and fluid move with the same velocity, and the other resulting in a Darcy contribution consisting of the relative fibre/fluid velocity. In this paper, we elaborate a general model able to adapt continuously to all these flow regimes.Data-driven upscaling of orientation kinematics in suspensions of rigid fibres
http://hdl.handle.net/10985/15419
Data-driven upscaling of orientation kinematics in suspensions of rigid fibres
SCHEUER, Adrien; AMMAR, Amine; ABISSET-CHAVANNE, Emmanuelle; CUETO, Elías; CHINESTA, Francisco; KEUNINGS, Roland; ADVANI, Suresh G.
Describing the orientation state of the particles is often critical in fibre suspension applications. Macroscopic descriptors, the so-called second-order orientation tensor (or moment) leading the way, are often preferred due to their low computational cost. Closure problems however arise when evolution equations for the moments are derived from the orientation distribution functions and the impact of the chosen closure is often unpredictable. In this work, our aim is to provide macroscopic simulations of orientation that are cheap, accurate and closure-free. To this end, we propose an innovative data-based approach to the upscaling of orientation kinematics in the context of fibre suspensions. Since the physics at the microscopic scale can be modelled reasonably enough, the idea is to conduct accurate offline direct numerical simulations at that scale and to extract the corresponding macroscopic descriptors in order to build a database of scenarios. During the online stage, the macroscopic descriptors can then be updated quickly by combining adequately the items from the database instead of relying on an imprecise macroscopic model. This methodology is presented in the well-known case of dilute fibre suspensions (where it can be compared against closure-based macroscopic models) and in the case of suspensions of confined or electrically-charged fibres, for which state-of-the-art closures proved to be inadequate or simply do not exist.
Mon, 01 Jan 2018 00:00:00 GMThttp://hdl.handle.net/10985/154192018-01-01T00:00:00ZSCHEUER, AdrienAMMAR, AmineABISSET-CHAVANNE, EmmanuelleCUETO, ElíasCHINESTA, FranciscoKEUNINGS, RolandADVANI, Suresh G.Describing the orientation state of the particles is often critical in fibre suspension applications. Macroscopic descriptors, the so-called second-order orientation tensor (or moment) leading the way, are often preferred due to their low computational cost. Closure problems however arise when evolution equations for the moments are derived from the orientation distribution functions and the impact of the chosen closure is often unpredictable. In this work, our aim is to provide macroscopic simulations of orientation that are cheap, accurate and closure-free. To this end, we propose an innovative data-based approach to the upscaling of orientation kinematics in the context of fibre suspensions. Since the physics at the microscopic scale can be modelled reasonably enough, the idea is to conduct accurate offline direct numerical simulations at that scale and to extract the corresponding macroscopic descriptors in order to build a database of scenarios. During the online stage, the macroscopic descriptors can then be updated quickly by combining adequately the items from the database instead of relying on an imprecise macroscopic model. This methodology is presented in the well-known case of dilute fibre suspensions (where it can be compared against closure-based macroscopic models) and in the case of suspensions of confined or electrically-charged fibres, for which state-of-the-art closures proved to be inadequate or simply do not exist.Microscopic modelling of orientation kinematics of non-spherical particles suspended in confined flows using unilateral mechanics
http://hdl.handle.net/10985/13304
Microscopic modelling of orientation kinematics of non-spherical particles suspended in confined flows using unilateral mechanics
SCHEUER, Adrien; ABISSET-CHAVANNE, Emmanuelle; CHINESTA, Francisco; KEUNINGS, Roland
The properties of reinforced polymers strongly depend on the microstructural state, that is, the orientation state of the fibres suspended in the polymeric matrix, induced by the forming process. Understanding flow-induced anisotropy is thus a key element to optimize both materials and process. Despite the important progresses accomplished in the modelling and simulation of suspensions, few works addressed the fact that usual processing flows evolve in confined configurations, where particles characteristic lengths may be greater than the thickness of the narrow gaps in which the flow takes place. In those circumstances, orientation kinematics models proposed for unconfined flows must be extended to the confined case. In this short communication, we propose an alternative modelling framework based on the use of unilateral mechanics, consequently exhibiting a clear analogy with plasticity and contact mechanics. This framework allows us to revisit the motion of confined particles in Newtonian and non-Newtonian matrices. We also prove that the confined kinematics provided by this model are identical to those derived from microstructural approaches
Mon, 01 Jan 2018 00:00:00 GMThttp://hdl.handle.net/10985/133042018-01-01T00:00:00ZSCHEUER, AdrienABISSET-CHAVANNE, EmmanuelleCHINESTA, FranciscoKEUNINGS, RolandThe properties of reinforced polymers strongly depend on the microstructural state, that is, the orientation state of the fibres suspended in the polymeric matrix, induced by the forming process. Understanding flow-induced anisotropy is thus a key element to optimize both materials and process. Despite the important progresses accomplished in the modelling and simulation of suspensions, few works addressed the fact that usual processing flows evolve in confined configurations, where particles characteristic lengths may be greater than the thickness of the narrow gaps in which the flow takes place. In those circumstances, orientation kinematics models proposed for unconfined flows must be extended to the confined case. In this short communication, we propose an alternative modelling framework based on the use of unilateral mechanics, consequently exhibiting a clear analogy with plasticity and contact mechanics. This framework allows us to revisit the motion of confined particles in Newtonian and non-Newtonian matrices. We also prove that the confined kinematics provided by this model are identical to those derived from microstructural approachesOn the multi‑scale description of electrical conducting suspensions involving perfectly dispersed rods
http://hdl.handle.net/10985/10253
On the multi‑scale description of electrical conducting suspensions involving perfectly dispersed rods
PEREZ, Marta; ABISSET-CHAVANNE, Emmanuelle; BARASINSKI, Anais; CHINESTA, Francisco; AMMAR, Amine; KEUNINGS, Roland
Nanocomposites allow for a significant enhancement of functional properties, in particular electrical conduction. In order to optimize materials and parts, predictive models are required to evaluate particle distribution and orientation. Both are key parameters in order to evaluate percolation and the resulting electrical networks. Many forming processes involve flowing suspensions for which the final particle orientation could be controlled by means of the flow and the electric field. In view of the multiscale character of the problem, detailed descriptions are defined at the microscopic scale and then coarsened to be applied efficiently in process simulation at the macroscopic scale. The first part of this work revisits the different modeling approaches throughout the different description scales. Then, modeling of particle contacts is addressed as they determine the final functional properties, in particular electrical conduction. Different descriptors of rod contacts are proposed and analyzed. Numerical results are discussed, in particular to evaluate the impact of closure approximations needed to derive a macroscopic description.
Thu, 01 Jan 2015 00:00:00 GMThttp://hdl.handle.net/10985/102532015-01-01T00:00:00ZPEREZ, MartaABISSET-CHAVANNE, EmmanuelleBARASINSKI, AnaisCHINESTA, FranciscoAMMAR, AmineKEUNINGS, RolandNanocomposites allow for a significant enhancement of functional properties, in particular electrical conduction. In order to optimize materials and parts, predictive models are required to evaluate particle distribution and orientation. Both are key parameters in order to evaluate percolation and the resulting electrical networks. Many forming processes involve flowing suspensions for which the final particle orientation could be controlled by means of the flow and the electric field. In view of the multiscale character of the problem, detailed descriptions are defined at the microscopic scale and then coarsened to be applied efficiently in process simulation at the macroscopic scale. The first part of this work revisits the different modeling approaches throughout the different description scales. Then, modeling of particle contacts is addressed as they determine the final functional properties, in particular electrical conduction. Different descriptors of rod contacts are proposed and analyzed. Numerical results are discussed, in particular to evaluate the impact of closure approximations needed to derive a macroscopic description.Flow modelling of quasi-Newtonian fluids in two-scale fibrous fabrics: Advanced simulations
http://hdl.handle.net/10985/11390
Flow modelling of quasi-Newtonian fluids in two-scale fibrous fabrics: Advanced simulations
AMMAR, Amine; ABISSET-CHAVANNE, Emmanuelle; CHINESTA, Francisco; KEUNINGS, Roland
Permeability is the fundamental macroscopic material property needed to quantify the flow in a fibrous medium viewed as a porous medium. Composite processing models require the permeability as input data to predict flow patterns and pressure fields. In a previous work, the expressions of macroscopic permeability were derived in a double-scale porosity medium for both Newtonian and generalized Newtonian (shear-thinning) resins. In the linear case, only a microscopic calculation on a representative volume is required, implying as many microscopic calculations as there are representative microscopic volumes in the whole fibrous structure. In the non-linear case, and even when the porous microstructure can be described by a unique representative volume, a large number of microscopic calculations must be carried out as the microscale resin viscosity depends on the macroscopic velocity, which in turn depends on the permeability that results from a microscopic calculation. An original and efficient offline-online procedure was proposed for the solution of non-linear flow problems related to generalized Newtonian fluids in porous media. In this paper, this procedure is generalized to quasi-Newtonian fluids in order to evaluate the effect of extensional viscosity on the resulting upscaled permeability. This work constitutes a natural step forward in the definition of equivalent saturated permeabilities for linear and non-linear fluids.
Fri, 01 Jan 2016 00:00:00 GMThttp://hdl.handle.net/10985/113902016-01-01T00:00:00ZAMMAR, AmineABISSET-CHAVANNE, EmmanuelleCHINESTA, FranciscoKEUNINGS, RolandPermeability is the fundamental macroscopic material property needed to quantify the flow in a fibrous medium viewed as a porous medium. Composite processing models require the permeability as input data to predict flow patterns and pressure fields. In a previous work, the expressions of macroscopic permeability were derived in a double-scale porosity medium for both Newtonian and generalized Newtonian (shear-thinning) resins. In the linear case, only a microscopic calculation on a representative volume is required, implying as many microscopic calculations as there are representative microscopic volumes in the whole fibrous structure. In the non-linear case, and even when the porous microstructure can be described by a unique representative volume, a large number of microscopic calculations must be carried out as the microscale resin viscosity depends on the macroscopic velocity, which in turn depends on the permeability that results from a microscopic calculation. An original and efficient offline-online procedure was proposed for the solution of non-linear flow problems related to generalized Newtonian fluids in porous media. In this paper, this procedure is generalized to quasi-Newtonian fluids in order to evaluate the effect of extensional viscosity on the resulting upscaled permeability. This work constitutes a natural step forward in the definition of equivalent saturated permeabilities for linear and non-linear fluids.A simple microstructural viscoelastic model for flowing foams
http://hdl.handle.net/10985/16835
A simple microstructural viscoelastic model for flowing foams
IBÁÑEZ, Rubén; SCHEUER, Adrien; ABISSET-CHAVANNE, Emmanuelle; CHINESTA, Francisco; HUERTA, Antonio; KEUNINGS, Roland
The numerical modelling of forming processes involving the flow of foams requires taking into account the different problem scales. Thus, in industrial applications a macroscopic approach is suitable, whereas the macroscopic flow parameters depend on the cellular structure: cell size, shape, orientation, etc. Moreover, the shape and orientation of the cells are induced by the flow. A fully microscopic description remains useful to understand the foam behaviour and the topological changes induced by the cell elongation or distortion, however, from an industrial point of view, microscopic simulations remain challenging to address practical applications involving flows in complex 3D geometries. In this paper, we propose a viscoelastic flow model where the foam microstructure is represented from suitable microstructure descriptors whose evolution is governed by the macroscopic flow kinematics.
Mon, 01 Jan 2018 00:00:00 GMThttp://hdl.handle.net/10985/168352018-01-01T00:00:00ZIBÁÑEZ, RubénSCHEUER, AdrienABISSET-CHAVANNE, EmmanuelleCHINESTA, FranciscoHUERTA, AntonioKEUNINGS, RolandThe numerical modelling of forming processes involving the flow of foams requires taking into account the different problem scales. Thus, in industrial applications a macroscopic approach is suitable, whereas the macroscopic flow parameters depend on the cellular structure: cell size, shape, orientation, etc. Moreover, the shape and orientation of the cells are induced by the flow. A fully microscopic description remains useful to understand the foam behaviour and the topological changes induced by the cell elongation or distortion, however, from an industrial point of view, microscopic simulations remain challenging to address practical applications involving flows in complex 3D geometries. In this paper, we propose a viscoelastic flow model where the foam microstructure is represented from suitable microstructure descriptors whose evolution is governed by the macroscopic flow kinematics.