SAM
https://sam.ensam.eu:443
The DSpace digital repository system captures, stores, indexes, preserves, and distributes digital research material.Fri, 30 Oct 2020 21:16:36 GMT2020-10-30T21:16:36ZReview on the Brownian Dynamics Simulation of Bead-Rod-Spring Models Encountered in Computational Rheology
http://hdl.handle.net/10985/9991
Review on the Brownian Dynamics Simulation of Bead-Rod-Spring Models Encountered in Computational Rheology
CRUZ, Camilo; CHINESTA, Francisco; REGNIER, Gilles
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 GMThttp://hdl.handle.net/10985/99912012-01-01T00:00:00ZCRUZ, CamiloCHINESTA, FranciscoREGNIER, GillesKinetic 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.Shear-strain step response in linear regime of dilute suspensions of naturally bent carbon nanotubes
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; CHINESTA, Francisco; REGNIER, Gilles
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 GMThttp://hdl.handle.net/10985/68002012-01-01T00:00:00ZCRUZ, CamiloILLOUL, LounèsCHINESTA, FranciscoREGNIER, GillesImpressive 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.Frontal weld lines in injection-molded short fiber-reinforced PBT: Extensive microstructure characterization for mechanical performance evaluation
http://hdl.handle.net/10985/18406
Frontal weld lines in injection-molded short fiber-reinforced PBT: Extensive microstructure characterization for mechanical performance evaluation
BARADI, Mohamed Besher; CRUZ, Camilo; RIEDEL, Thomas; RÉGNIER, Gilles
Different factors contribute to the weakness of weld lines (WLs) induced by injection molding such as unsuitable fiber orientation (FO), incomplete polymer matrix diffusion, voids and V-notches. This study aims to characterize the contribution of each factor on the weakness of frontal WLs in a short glass fiber-reinforced polybutylene-terephthalate characterized by extensive X-ray computed tomography and mechanical tensile testing assisted with digital image correlation. A reduction of 50% of the stress at break and almost 40% of the strain at break is observed despite the complete matrix healing at the WL interface and the absence of V-notches. Frontal WLs induce a FO gradient starting 2 to 3 mm before the WL plane. The fibers in the WL region mainly orient in transverse-to-flow and thickness direction. This FO gradient localizes the deformations, which leads to failure at a strength near to the one of the unreinforced variant. Voids formation in frontal WLs seems to be driven by large gradients of FO and subsequent anisotropic shrinkage. In addition, this shrinkage behavior at the WL causes an increase of thickness. By applying higher packing pressures, the fibers orient more in flow direction at the core of the WL, leading to a higher tensile strength and a lower content of voids. Finally, we can conclude that the FO is the dominant factor controlling the mechanical performance in frontal WLs.
Tue, 01 Jan 2019 00:00:00 GMThttp://hdl.handle.net/10985/184062019-01-01T00:00:00ZBARADI, Mohamed BesherCRUZ, CamiloRIEDEL, ThomasRÉGNIER, GillesDifferent factors contribute to the weakness of weld lines (WLs) induced by injection molding such as unsuitable fiber orientation (FO), incomplete polymer matrix diffusion, voids and V-notches. This study aims to characterize the contribution of each factor on the weakness of frontal WLs in a short glass fiber-reinforced polybutylene-terephthalate characterized by extensive X-ray computed tomography and mechanical tensile testing assisted with digital image correlation. A reduction of 50% of the stress at break and almost 40% of the strain at break is observed despite the complete matrix healing at the WL interface and the absence of V-notches. Frontal WLs induce a FO gradient starting 2 to 3 mm before the WL plane. The fibers in the WL region mainly orient in transverse-to-flow and thickness direction. This FO gradient localizes the deformations, which leads to failure at a strength near to the one of the unreinforced variant. Voids formation in frontal WLs seems to be driven by large gradients of FO and subsequent anisotropic shrinkage. In addition, this shrinkage behavior at the WL causes an increase of thickness. By applying higher packing pressures, the fibers orient more in flow direction at the core of the WL, leading to a higher tensile strength and a lower content of voids. Finally, we can conclude that the FO is the dominant factor controlling the mechanical performance in frontal WLs.Mechanical and microstructural characterization of flowing weld lines in injection-molded short fiber-reinforced PBT
http://hdl.handle.net/10985/14636
Mechanical and microstructural characterization of flowing weld lines in injection-molded short fiber-reinforced PBT
BARADI, Mohamed Besher; CRUZ, Camilo; RIEDEL, Thomas; REGNIER, Gilles
The aim of this work is an extensive experimental mechanical and microstructural characterization of flowing weld lines (WLs) in injection-molded short glass fiber-reinforced polybutylenterephthalate (PBT) using X-ray computed tomography and digital image correlation (DIC). It was found that the induced fiber orientation (FO) in a flowing WL is similar to that induced by flow along a wall. In this test case, the impact of the flowing WL on the FO did not vanish after a flow length of 70 mm. The shape of the inserts, which originated the flowing WLs, only had a marginal effect on the induced FO gradient. By reducing part thickness, the erasing of the FO gradient induced by the WL is reached at shorter flow distances. At the WLs, there is a reduction of the fiber volume fraction in comparison to the regions far from the WL plane. DIC results show a pronounced strain localization at the WL, which can be explained by the FO gradient induced by the WL.
Tue, 01 Jan 2019 00:00:00 GMThttp://hdl.handle.net/10985/146362019-01-01T00:00:00ZBARADI, Mohamed BesherCRUZ, CamiloRIEDEL, ThomasREGNIER, GillesThe aim of this work is an extensive experimental mechanical and microstructural characterization of flowing weld lines (WLs) in injection-molded short glass fiber-reinforced polybutylenterephthalate (PBT) using X-ray computed tomography and digital image correlation (DIC). It was found that the induced fiber orientation (FO) in a flowing WL is similar to that induced by flow along a wall. In this test case, the impact of the flowing WL on the FO did not vanish after a flow length of 70 mm. The shape of the inserts, which originated the flowing WLs, only had a marginal effect on the induced FO gradient. By reducing part thickness, the erasing of the FO gradient induced by the WL is reached at shorter flow distances. At the WLs, there is a reduction of the fiber volume fraction in comparison to the regions far from the WL plane. DIC results show a pronounced strain localization at the WL, which can be explained by the FO gradient induced by the WL.Effects of a bent structure on the linear viscoelastic response of diluted carbon nanotube suspensions
http://hdl.handle.net/10985/17960
Effects of a bent structure on the linear viscoelastic response of diluted carbon nanotube suspensions
CRUZ, Camilo; ILLOUL, Lounès; CHINESTA, Francisco; RÉGNIER, Gilles
Commonly isolated carbon nanotubes in suspension have been modelled as a perfectly straight structure. Nevertheless, single-wall carbon nanotubes (SWNTs) contain naturally side-wall defects and, in consequence, natural bent configurations. Hence, a semi-flexile filament model with a natural bent configuration was proposed to represent physically the SWNT structure. This continuous model was discretized as a non-freely jointed multi-bead-rod system with a natural bent configuration. Using a Brownian dynamics algorithm the dynamical mechanical contribution to the linear viscoelastic response of naturally bent SWNTs in dilute suspension was simulated. The dynamics of such system shows the apparition of new relaxation processes at intermediate frequencies characterized mainly by the activation of a mild elasticity. Storage modulus evolution at those intermediate frequencies strongly depends on the flexibility of the system, given by the rigidity constant of the bending potential and the number of constitutive rods.
Fri, 01 Jan 2010 00:00:00 GMThttp://hdl.handle.net/10985/179602010-01-01T00:00:00ZCRUZ, CamiloILLOUL, LounèsCHINESTA, FranciscoRÉGNIER, GillesCommonly isolated carbon nanotubes in suspension have been modelled as a perfectly straight structure. Nevertheless, single-wall carbon nanotubes (SWNTs) contain naturally side-wall defects and, in consequence, natural bent configurations. Hence, a semi-flexile filament model with a natural bent configuration was proposed to represent physically the SWNT structure. This continuous model was discretized as a non-freely jointed multi-bead-rod system with a natural bent configuration. Using a Brownian dynamics algorithm the dynamical mechanical contribution to the linear viscoelastic response of naturally bent SWNTs in dilute suspension was simulated. The dynamics of such system shows the apparition of new relaxation processes at intermediate frequencies characterized mainly by the activation of a mild elasticity. Storage modulus evolution at those intermediate frequencies strongly depends on the flexibility of the system, given by the rigidity constant of the bending potential and the number of constitutive rods.