https://sam.ensam.eu:443 The DSpace digital repository system captures, stores, indexes, preserves, and distributes digital research material. Mon, 08 Jun 2026 14:49:24 GMT 2026-06-08T14:49:24Z http://hdl.handle.net/10985/8488 One and two-fiber orientation kinetic theories of fiber suspensions GRMELA, Miroslav; AMMAR, Amine; CHINESTA SORIA, Francisco The morphology influencing rheological properties of suspensions of rigid spheres constitutes the flow induced collective ordering of the spheres characterized by two or more sphere distribution functions. When the rigid spheres are replaced by rigid fibers, the collective order in the position of the spheres is replaced by the flow induced orientation of the fibers that suffices to be characterized by one-fiber orientation distribution function. A flow induced collective ordering of fibers (both in position and orientation), that can only be characterized by two or more fiber distribution functions, can still however constitute an important part of the morphology. We show that two types of interaction among fibers, one being the Onsager-type topological interaction entering the free energy and the other the hydrodynamics interaction entering the dissipative part of the time evolution, give indeed rise to a collective order in the orientation influencing the rheology of fiber suspensions. http://dx.doi.org/10.1016/j.jnnfm.2012.10.009 Tue, 01 Jan 2013 00:00:00 GMT http://hdl.handle.net/10985/8488 2013-01-01T00:00:00Z GRMELA, Miroslav AMMAR, Amine CHINESTA SORIA, Francisco The morphology influencing rheological properties of suspensions of rigid spheres constitutes the flow induced collective ordering of the spheres characterized by two or more sphere distribution functions. When the rigid spheres are replaced by rigid fibers, the collective order in the position of the spheres is replaced by the flow induced orientation of the fibers that suffices to be characterized by one-fiber orientation distribution function. A flow induced collective ordering of fibers (both in position and orientation), that can only be characterized by two or more fiber distribution functions, can still however constitute an important part of the morphology. We show that two types of interaction among fibers, one being the Onsager-type topological interaction entering the free energy and the other the hydrodynamics interaction entering the dissipative part of the time evolution, give indeed rise to a collective order in the orientation influencing the rheology of fiber suspensions. http://hdl.handle.net/10985/8473 An overview of the proper generalized decomposition with applications in computational rheology CHINESTA SORIA, 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 GMT http://hdl.handle.net/10985/8473 2011-01-01T00:00:00Z CHINESTA SORIA, 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. http://hdl.handle.net/10985/10244 Parametric solutions involving geometry: A step towards efficient shape optimization AMMAR, Amine; HUERTA, Antonio; CHINESTA SORIA, Francisco; CUETO, Elias; LEYGUE, Adrien Optimization of manufacturing processes or structures involves the optimal choice of many parameters (process parameters, material parameters or geometrical parameters). Usual strategies proceed by defining a trial choice of those parameters and then solving the resulting model. Then, an appropriate cost function is evaluated and its optimality checked. While the optimum is not reached, the process parameters should be updated by using an appropriate optimization procedure, and then the model must be solved again for the updated process parameters. Thus, a direct numerical solution is needed for each choice of the process parameters, with the subsequent impact on the computing time. In this work we focus on shape optimization that involves the appropriate choice of some parameters defining the problem geometry. The main objective of this work is to describe an original approach for computing an off-line parametric solution. That is, a solution able to include information for different parameter values and also allowing to compute readily the sensitivities. The curse of dimensionality is circumvented by invoking the Proper Generalized Decomposition (PGD) introduced in former works, which is applied here to compute geometrically parametrized solutions. Wed, 01 Jan 2014 00:00:00 GMT http://hdl.handle.net/10985/10244 2014-01-01T00:00:00Z AMMAR, Amine HUERTA, Antonio CHINESTA SORIA, Francisco CUETO, Elias LEYGUE, Adrien Optimization of manufacturing processes or structures involves the optimal choice of many parameters (process parameters, material parameters or geometrical parameters). Usual strategies proceed by defining a trial choice of those parameters and then solving the resulting model. Then, an appropriate cost function is evaluated and its optimality checked. While the optimum is not reached, the process parameters should be updated by using an appropriate optimization procedure, and then the model must be solved again for the updated process parameters. Thus, a direct numerical solution is needed for each choice of the process parameters, with the subsequent impact on the computing time. In this work we focus on shape optimization that involves the appropriate choice of some parameters defining the problem geometry. The main objective of this work is to describe an original approach for computing an off-line parametric solution. That is, a solution able to include information for different parameter values and also allowing to compute readily the sensitivities. The curse of dimensionality is circumvented by invoking the Proper Generalized Decomposition (PGD) introduced in former works, which is applied here to compute geometrically parametrized solutions. http://hdl.handle.net/10985/6800 Shear-strain step response in linear regime of dilute suspensions of naturally bent carbon nanotubes CRUZ, Camilo; ILLOUL, Lounès; REGNIER, Gilles; CHINESTA SORIA, Francisco Impressive enhancements of the storage modulus have been documented when low volume fractions of single wall carbon nanotubes (SWNTs) are added to a Newtonian solvent for obtaining dilute suspensions. The intrinsic bending dynamics of carbon nanotubes (CNTs) has been proposed to explain such elasticity. CNTs contain topological defects inducing naturally bent structures in absence of external forces and, hence, a semiflexible filament with a bent configuration at minimal internal-bending-energy is used for mimicking the structure of SWNTs in suspension. Previous continuous model is discretized as a non-freely jointed bead-rod chain with a naturally bent configuration for simulating the rheological behaviour after a shear-strain step in linear regime of SWNT dilute suspension by using a Brownian dynamics (BD) approach. In general, bead-rod chains exhibit an instantaneous relaxation after a high shear-strain step. Bending rigidity and number of constitutive rods are found to be determinant parameters in the internal-energy relaxation behaviour of non-freely jointed bead-rod chains in dilute solution. Proper comparisons between the BD simulation results and the experimental data for treated SWNT dilute suspensions confirm the consistency of the physical model mimicking the structure of a SWNT. Sun, 01 Jan 2012 00:00:00 GMT http://hdl.handle.net/10985/6800 2012-01-01T00:00:00Z CRUZ, Camilo ILLOUL, Lounès REGNIER, Gilles CHINESTA SORIA, Francisco Impressive enhancements of the storage modulus have been documented when low volume fractions of single wall carbon nanotubes (SWNTs) are added to a Newtonian solvent for obtaining dilute suspensions. The intrinsic bending dynamics of carbon nanotubes (CNTs) has been proposed to explain such elasticity. CNTs contain topological defects inducing naturally bent structures in absence of external forces and, hence, a semiflexible filament with a bent configuration at minimal internal-bending-energy is used for mimicking the structure of SWNTs in suspension. Previous continuous model is discretized as a non-freely jointed bead-rod chain with a naturally bent configuration for simulating the rheological behaviour after a shear-strain step in linear regime of SWNT dilute suspension by using a Brownian dynamics (BD) approach. In general, bead-rod chains exhibit an instantaneous relaxation after a high shear-strain step. Bending rigidity and number of constitutive rods are found to be determinant parameters in the internal-energy relaxation behaviour of non-freely jointed bead-rod chains in dilute solution. Proper comparisons between the BD simulation results and the experimental data for treated SWNT dilute suspensions confirm the consistency of the physical model mimicking the structure of a SWNT. http://hdl.handle.net/10985/10242 Parametric solution of the Rayleigh-Benard convection model by using the PGD Application to nanofluids AGHIGHI, Mohammad Saeid; AMMAR, Amine; METIVIER, Christel; CHINESTA SORIA, Francisco Purpose – The purpose of this paper is to focus on the advanced solution of the parametric non-linear model related to the Rayleigh-Benard laminar flow involved in the modeling of natural thermal convection. This flow is fully determined by the dimensionless Prandtl and Rayleigh numbers. Thus, if one could precompute (off-line) the model solution for any possible choice of these two parameters the analysis of many possible scenarios could be performed on-line and in real time. Design/methodology/approach – In this paper both parameters are introduced as model extracoordinates, and then the resulting multidimensional problem solved thanks to the space-parameters separated representation involved in the proper generalized decomposition (PGD) that allows circumventing the curse of dimensionality. Thus the parametric solution will be available fast and easily. Findings – Such parametric solution could be viewed as a sort of abacus, but despite its inherent interest such calculation is at present unaffordable for nowadays computing availabilities because one must solve too many problems and of course store all the solutions related to each choice of both parameters. Originality/value – Parametric solution of coupled models by using the PGD. Model reduction of complex coupled flow models. Analysis of Rayleigh-Bernard flows involving nanofluids. Thu, 01 Jan 2015 00:00:00 GMT http://hdl.handle.net/10985/10242 2015-01-01T00:00:00Z AGHIGHI, Mohammad Saeid AMMAR, Amine METIVIER, Christel CHINESTA SORIA, Francisco Purpose – The purpose of this paper is to focus on the advanced solution of the parametric non-linear model related to the Rayleigh-Benard laminar flow involved in the modeling of natural thermal convection. This flow is fully determined by the dimensionless Prandtl and Rayleigh numbers. Thus, if one could precompute (off-line) the model solution for any possible choice of these two parameters the analysis of many possible scenarios could be performed on-line and in real time. Design/methodology/approach – In this paper both parameters are introduced as model extracoordinates, and then the resulting multidimensional problem solved thanks to the space-parameters separated representation involved in the proper generalized decomposition (PGD) that allows circumventing the curse of dimensionality. Thus the parametric solution will be available fast and easily. Findings – Such parametric solution could be viewed as a sort of abacus, but despite its inherent interest such calculation is at present unaffordable for nowadays computing availabilities because one must solve too many problems and of course store all the solutions related to each choice of both parameters. Originality/value – Parametric solution of coupled models by using the PGD. Model reduction of complex coupled flow models. Analysis of Rayleigh-Bernard flows involving nanofluids. http://hdl.handle.net/10985/8479 On the solution of the multidimensional Langer’s equation using the proper generalized decomposition method for modeling phase transitions LAMARI, Hajer; AMMAR, Amine; LEYGUE, Adrien; CHINESTA SORIA, Francisco The dynamics of phase transition in a binary mixture occurring during a quench is studied taking into account composition fluctuations by solving Langer’s equation in a domain composed of a certain number of micro-domains. The resulting Langer’s equation governing the evolution of the distribution function becomes multidimensional. Circumventing the curse of dimensionality the proper generalized decomposition is applied. The influence of the interaction parameter in the vicinity of the critical point is analyzed. First we address the case of a system composed of a single micro-domain in which phase transition occurs by a simple symmetry change. Next, we consider a system composed of two micro-domains in which phase transition occurs by phase separation, with special emphasis on the effect of the Landau free energy non-local term. Finally, some systems consisting of many micro-domains are considered. http://dx.doi.org/10.1088/0965-0393/20/1/015007 Sun, 01 Jan 2012 00:00:00 GMT http://hdl.handle.net/10985/8479 2012-01-01T00:00:00Z LAMARI, Hajer AMMAR, Amine LEYGUE, Adrien CHINESTA SORIA, Francisco The dynamics of phase transition in a binary mixture occurring during a quench is studied taking into account composition fluctuations by solving Langer’s equation in a domain composed of a certain number of micro-domains. The resulting Langer’s equation governing the evolution of the distribution function becomes multidimensional. Circumventing the curse of dimensionality the proper generalized decomposition is applied. The influence of the interaction parameter in the vicinity of the critical point is analyzed. First we address the case of a system composed of a single micro-domain in which phase transition occurs by a simple symmetry change. Next, we consider a system composed of two micro-domains in which phase transition occurs by phase separation, with special emphasis on the effect of the Landau free energy non-local term. Finally, some systems consisting of many micro-domains are considered. http://hdl.handle.net/10985/8472 Simulating microstructure evolution during passive mixing MAITREJEAN, Guillaume; AMMAR, Amine; CHINESTA SORIA, Francisco The prediction of microstructure evolution during passive mixing is of major interest in order to qualify and quantify mixing devices as well as to predict the final morphology of the resulting blend. Direct numerical simulation fails because of the different characteristic lengths of the microstructure and the process itself. Micro-macro approaches could be a valuable alternative but the computational cost remains tremendous. For this reason many authors proposed the introduction of some microstructural variables able to qualify and quantify the mixing process at a mesoscale level. Some proposals considered only the effects induced by the flow kinematics, other introduced also the effects of shape relaxation due to the surface tension and coalescence. The most advanced integrate also the break-up process. However, the derivation of the evolution equations governing the evolution of such microstructural variables needs the introduction of some closure relations whose impact on the computed solution should be evaluated before applying it for simulating complex mixing flows. In this work we consider the Lee and Park’s model that considers the flow kinematics, the surface tension, the coalescence and the break-up mechanisms in the evolution of the area tensor. The accuracy of both a quadratic closure and an orthotropic relations will be analyzed in the first part of this work, and then the resulting closed model by using a quadratic closure will be used for simulating complex mixing flows. Sun, 01 Jan 2012 00:00:00 GMT http://hdl.handle.net/10985/8472 2012-01-01T00:00:00Z MAITREJEAN, Guillaume AMMAR, Amine CHINESTA SORIA, Francisco The prediction of microstructure evolution during passive mixing is of major interest in order to qualify and quantify mixing devices as well as to predict the final morphology of the resulting blend. Direct numerical simulation fails because of the different characteristic lengths of the microstructure and the process itself. Micro-macro approaches could be a valuable alternative but the computational cost remains tremendous. For this reason many authors proposed the introduction of some microstructural variables able to qualify and quantify the mixing process at a mesoscale level. Some proposals considered only the effects induced by the flow kinematics, other introduced also the effects of shape relaxation due to the surface tension and coalescence. The most advanced integrate also the break-up process. However, the derivation of the evolution equations governing the evolution of such microstructural variables needs the introduction of some closure relations whose impact on the computed solution should be evaluated before applying it for simulating complex mixing flows. In this work we consider the Lee and Park’s model that considers the flow kinematics, the surface tension, the coalescence and the break-up mechanisms in the evolution of the area tensor. The accuracy of both a quadratic closure and an orthotropic relations will be analyzed in the first part of this work, and then the resulting closed model by using a quadratic closure will be used for simulating complex mixing flows. http://hdl.handle.net/10985/9991 Review on the Brownian Dynamics Simulation of Bead-Rod-Spring Models Encountered in Computational Rheology CRUZ, Camilo; REGNIER, Gilles; CHINESTA SORIA, Francisco Kinetic theory is a mathematical framework intended to relate directly the most relevant characteristics of the molecular structure to the rheological behavior of the bulk system. In other words, kinetic theory is a micro-to-macro approach for solving the flow of complex fluids that circumvents the use of closure relations and offers a better physical description of the phenomena involved in the flow processes. Cornerstone models in kinetic theory employ beads, rods and springs for mimicking the molecular structure of the complex fluid. The generalized bead-rod-spring chain includes the most basic models in kinetic theory: the freely jointed bead-spring chain and the freely-jointed bead-rod chain. Configuration of simple coarse-grained models can be represented by an equivalent Fokker-Planck (FP) diffusion equation, which describes the evolution of the configuration distribution function in the physical and configurational spaces. FP equation can be a complex mathematical object, given its multidimensionality, and solving it explicitly can become a difficult task. Even more, in some cases, obtaining an equivalent FP equation is not possible given the complexity of the coarse-grained molecular model. Brownian dynamics can be employed as an alternative extensive numerical method for approaching the configuration distribution function of a given kinetic-theory model that avoid obtaining and/or resolving explicitly an equivalent FP equation. The validity of this discrete approach is based on the mathematical equivalence between a continuous diffusion equation and a stochastic differential equation as demonstrated by Itô in the 1940s. This paper presents a review of the fundamental issues in the BD simulation of the linear viscoelastic behavior of bead-rod-spring coarse grained models in dilute solution. In the first part of this work, the BD numerical technique is introduced. An overview of the mathematical framework of the BD and a review of the scope of applications are presented. Subsequently, the links between the rheology of complex fluids, the kinetic theory and the BD technique are established at the light of the stochastic nature of the bead-rod-spring models. Finally, the pertinence of the present state-of-the-art review is explained in terms of the increasing interest for the stochastic micro-to-macro approaches for solving complex fluids problems. In the second part of this paper, a detailed description of the BD algorithm used for simulating a small-amplitude oscillatory deformation test is given. Dynamic properties are employed throughout this work to characterise the linear viscoelastic behavior of bead-rod-spring models in dilute solution. In the third and fourth part of this article, an extensive discussion about the main issues of a BD simulation in linear viscoelasticity of diluted suspensions is tackled at the light of the classical multi-bead-spring chain model and the multi-bead-rod chain model, respectively. Kinematic formulations, integration schemes and expressions to calculate the stress tensor are revised for several classical models: Rouse and Zimm theories in the case of multi-bead-spring chains, and Kramers chain and semi-flexible filaments in the case of multi-bead-rod chains. The implemented BD technique is, on the one hand, validated in front of the analytical or exact numerical solutions known of the equivalent FP equations for those classic kinetic theory models; and, on the other hand, is control-set thanks to the analysis of the main numerical issues involved in a BD simulation. Finally, the review paper is closed by some concluding remarks. Sun, 01 Jan 2012 00:00:00 GMT http://hdl.handle.net/10985/9991 2012-01-01T00:00:00Z CRUZ, Camilo REGNIER, Gilles CHINESTA SORIA, Francisco Kinetic theory is a mathematical framework intended to relate directly the most relevant characteristics of the molecular structure to the rheological behavior of the bulk system. In other words, kinetic theory is a micro-to-macro approach for solving the flow of complex fluids that circumvents the use of closure relations and offers a better physical description of the phenomena involved in the flow processes. Cornerstone models in kinetic theory employ beads, rods and springs for mimicking the molecular structure of the complex fluid. The generalized bead-rod-spring chain includes the most basic models in kinetic theory: the freely jointed bead-spring chain and the freely-jointed bead-rod chain. Configuration of simple coarse-grained models can be represented by an equivalent Fokker-Planck (FP) diffusion equation, which describes the evolution of the configuration distribution function in the physical and configurational spaces. FP equation can be a complex mathematical object, given its multidimensionality, and solving it explicitly can become a difficult task. Even more, in some cases, obtaining an equivalent FP equation is not possible given the complexity of the coarse-grained molecular model. Brownian dynamics can be employed as an alternative extensive numerical method for approaching the configuration distribution function of a given kinetic-theory model that avoid obtaining and/or resolving explicitly an equivalent FP equation. The validity of this discrete approach is based on the mathematical equivalence between a continuous diffusion equation and a stochastic differential equation as demonstrated by Itô in the 1940s. This paper presents a review of the fundamental issues in the BD simulation of the linear viscoelastic behavior of bead-rod-spring coarse grained models in dilute solution. In the first part of this work, the BD numerical technique is introduced. An overview of the mathematical framework of the BD and a review of the scope of applications are presented. Subsequently, the links between the rheology of complex fluids, the kinetic theory and the BD technique are established at the light of the stochastic nature of the bead-rod-spring models. Finally, the pertinence of the present state-of-the-art review is explained in terms of the increasing interest for the stochastic micro-to-macro approaches for solving complex fluids problems. In the second part of this paper, a detailed description of the BD algorithm used for simulating a small-amplitude oscillatory deformation test is given. Dynamic properties are employed throughout this work to characterise the linear viscoelastic behavior of bead-rod-spring models in dilute solution. In the third and fourth part of this article, an extensive discussion about the main issues of a BD simulation in linear viscoelasticity of diluted suspensions is tackled at the light of the classical multi-bead-spring chain model and the multi-bead-rod chain model, respectively. Kinematic formulations, integration schemes and expressions to calculate the stress tensor are revised for several classical models: Rouse and Zimm theories in the case of multi-bead-spring chains, and Kramers chain and semi-flexible filaments in the case of multi-bead-rod chains. The implemented BD technique is, on the one hand, validated in front of the analytical or exact numerical solutions known of the equivalent FP equations for those classic kinetic theory models; and, on the other hand, is control-set thanks to the analysis of the main numerical issues involved in a BD simulation. Finally, the review paper is closed by some concluding remarks. http://hdl.handle.net/10985/9962 A mesoscopic rheological model of moderately concentrated colloids GRMELA, Miroslav; AMMAR, Amine; CHINESTA SORIA, Francisco; MAITREJEAN, Guillaume We extend the Maffettone–Minale model by including non-elliptical shapes of dispersed particles, a new family of internal forces controlling particle deformations, and particle–particle interactions. The last extension is made by transposing the way the chain-chain interactions are mathematically expressed in the reptation theory to suspensions. The particle–particle interactions are regarded as a confinement to cages formed by surrounding particles and by introducing a new dissipative motion (an analog of the reptation motion) inside the cages. Nonlinear responses to imposed shear and elongational flows are found to be in qualitative agreement with available experimental data. Wed, 01 Jan 2014 00:00:00 GMT http://hdl.handle.net/10985/9962 2014-01-01T00:00:00Z GRMELA, Miroslav AMMAR, Amine CHINESTA SORIA, Francisco MAITREJEAN, Guillaume We extend the Maffettone–Minale model by including non-elliptical shapes of dispersed particles, a new family of internal forces controlling particle deformations, and particle–particle interactions. The last extension is made by transposing the way the chain-chain interactions are mathematically expressed in the reptation theory to suspensions. The particle–particle interactions are regarded as a confinement to cages formed by surrounding particles and by introducing a new dissipative motion (an analog of the reptation motion) inside the cages. Nonlinear responses to imposed shear and elongational flows are found to be in qualitative agreement with available experimental data. http://hdl.handle.net/10985/6476 Empirical Natural Closure Relation for Short Fiber Suspension Models PRULIERE, Etienne; AMMAR, Amine; CHINESTA SORIA, Francisco This work focuses on the resolution of the Fokker-Planck equation that governs the evolution of the fibers orientation distribution. To reduce the computing time, that equation is solved along some flow trajectories in order to extract the significant information of the solution from the application of the Karhunen-Loève decomposition. Now, from this information one could solve the Fokker-Planck equation everywhere in the flow domain or simply adjust a closure relation that becomes optimal for such flow, solving the evolution of some orientation moments which require a less amount of computation. This paper focuses on this last strategy. For this purpose we start introducing the Karhunen-Loève decomposition that is applied later to automatically extract the main solution characteristics for adjusting empirically a natural closure relation. L'auteur Francisco CHINESTA faisait parti en 2007 du Laboratoire de Mécanique des Systèmes et des Procédés (LMSP). Depuis 2010, le LMSP a fusionné avec deux autres unités de recherche en un seul laboratoire intitulé PIMM (Procédés et Ingénierie en Mécanique et Matériaux). Mon, 01 Jan 2007 00:00:00 GMT http://hdl.handle.net/10985/6476 2007-01-01T00:00:00Z PRULIERE, Etienne AMMAR, Amine CHINESTA SORIA, Francisco This work focuses on the resolution of the Fokker-Planck equation that governs the evolution of the fibers orientation distribution. To reduce the computing time, that equation is solved along some flow trajectories in order to extract the significant information of the solution from the application of the Karhunen-Loève decomposition. Now, from this information one could solve the Fokker-Planck equation everywhere in the flow domain or simply adjust a closure relation that becomes optimal for such flow, solving the evolution of some orientation moments which require a less amount of computation. This paper focuses on this last strategy. For this purpose we start introducing the Karhunen-Loève decomposition that is applied later to automatically extract the main solution characteristics for adjusting empirically a natural closure relation.